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ABSTRACT 


An  experimental  study  of  laminar  flow  heat  transfer  of  an  ethylene  glycoi/water 
mixture  in  an  electrically  heated  horizontal  tube  using  wire  mesh  (HEATEX)  and 
twisted  tape  inserts  was  investigated.  Twelve  thermocouples,  inserted  in  the  tube  wall 
at  four  longitudinal  locations,  enabled  a  mean  inside  experimental  heat-transfer 
coefficient  to  be  accurately  measured.  A  constant  wall  heat  flux  boundary  condition 
was  placed  on  the  wall  by  wrapping  six  200  W  flexible  heater  tapes  tightly  around  the 
tube.  The  ethylene  glycol/water  mixture  provided  a  coolant  Reynolds  number 
between  200-5000  and  a  Prandtl  number  between  30-140.  Two  smooth  inside 
diameters  and  a  roped  tube  profile  were  tested  with  and  without  the  inserts.  Heat- 
transfer  correlations  for  tubes  without  inserts  were  developed  and  compared  with 
theory  for  both  thermally  and  hydrodynamically  developing  flow.  Correlations  were 
also  developed  for  the  two  types  of  inserts.  Nusselt  numbers  for  fully  developed  flow 
were  found  to  be  a  function  of  Reynolds  and  Prandtl  numbers  for  the  wire  mesh 
insert  and  a  function  of  tape  twist  ratio,  Reynolds  and  Prandtl  numbers  for  the 
twisted  tape  insert  Heat  transfer  enhancements  of  over  7  for  the  wire  mesh  insert 
and  over  4  for  the  twisted  tape  insert  at  high  Reynolds  numbers  were  obtained  over 
the  empty  tube.  By  using  these  correlations  in  conjunction  with  earlier  obtained 
refrigerant  condensation  data  (using  the  same  tubes,  inserts,  and  coolant),  more 
accurate  values  of  the  outside  condensation  heat-transfer  coefficients  were  obtained. 
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NOMENCLATURE 


A 

area,  m2 

c/» 

specific  heat,  J/kg  K 

D 

diameter,  m 

Gz 

Graetz  number 

h 

heat  transfer  coefficient,  W/m^ 

H 

pitch  for  180°  rotation  of  tape 

I 

current,  amps 

k 

thermal  conductivity,  W/m  K 

L 

length,  m 

LMTD 

log  mean  temperature  difference 

m 

mass  flow  rate,  kg/s 

Nu 

Nusselt  number 

P 

perimeter,  m  or  Power,  W 

Pr 

Prandtl  number 

q 

heat  transfer  rate,  W 

heat  flux,  W/nt2 

r 

radius,  m 

R 

electrical  resistance,  ohms 

Rw 

wall  resistance,  K/W 

Xlll 


Re  Reynolds  number 

Re,  Reynolds  number  for  swirl  flow 

T  Temperature,  K 

t  fin  thickness,  m 

U  overall  heat  transfer  coefficient.  W/n?  K 

u  velocity  component  in  x  direction,  m/s 

v  velocity  component  in  y  direction,  m/s 

V  velocity,  m/s  or  Voltage,  volts 

x*  dimensionless  axial  coordinate,  (x/D)/(Re  Pr) 

y  twist  ratio,  H/D 

Greek  Symbols 

5  tape  thickness,  m  or  boundary  layer  thickness 

e  heat  transfer  enhancement  ratio 

li  dynamic  viscosity,  kg/m  s 

v  kinematic  viscosity,  m2/s 

p  density,  kg/m3 

Subscripts 

c  coolant  or  cross  section 

D  tube  diameter 

fd  fully  developed 

f  fin 

h  hydrodynamically 


xiv 


surface 


saturation 

thermally 

thermocouple 


I.  INTRODUCTION 


A.  BACKGROUND 

With  the  increasing  complexity  of  shipboard  weapons  and  combat  systems  used 
aboard  United  States  Naval  vessels,  the  need  for  more  energy  efficient,  lightweight, 
and  high  capacity  refrigeration,  air  conditioning,  and  chilled  water  systems  has  grown 
tremendously.  One  way  to  meet  this  need  is  to  redesign  the  refrigeration  condenser. 
In  redesigning  these  condensers,  changes  in  condenser  tube  material  can  be 
combined  with  various  heat  transfer  enhancement  methods  to  achieve  a  lighter 
weight,  more  efficient  condenser  tube. 

Emphasis  on  refrigerant  condensation  research  reached  a  much  higher  level  of 
importance  due  to  the  pressing  environmental  issues  of  ozone  depletion  and  global 
warming.  In  September  1987,  an  international  conference  was  held  in  Montreal, 
Canada,  in  which  24  nations  representing  the  United  Nation’s  Environment  Program 
(UNEP),  signed  the  Montreal  Protocol  on  Substances  that  Deplete  the  Ozone  Layer. 
This  agreement  called  for  a  near  term  freeze  on  the  production  and  consumption  of 
chlorofluorocarbons  (CFCs),  •  the  major  contributor  to  ozone  depletion 
[Ref.  1].  In  August  1988,  the  Environmental  Protection  Agency  (EPA) 
adopted  the  regulations  of  the  Montreal  Protocol  under  the  Clean  Air  Act  [Ref.  2]. 
These  regulations  required  a  freeze  in  both  production  and  consumption  of  CFC-11, 
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12,  113,  114  and  115  at  1986  levels.  The  U.S.  Navy  currently  uses  CFC- 12  and  CFC- 
114  in  its  shipboard  air  conditioning  and  refrigeration  systems.  In  June  1990  at  a 
progress  meeting  held  in  London,  England,  UNEP  delegates  agreed  to  even  larger 
reductions,  calling  for  a  50%  reduction  by  1995,  85%  reduction  by  1997,  and  for 
CFCs  to  be  totally  phased  out  by  the  year  2000  [Kef.  3].  In  the  spring  of  1992, 
President  Bush  moved  up  the  CFC  phase  out  deadline  to  1995. 

Because  of  the  Navy’s  commitment  to  the  CFC  phase  out  policy  and  the 
increasing  need  for  lighter  weight,  more  efficient  refrigeration  systems,  refrigerant 
condensation  research  has  surfaced  to  the  forefront  of  Naval  Research. 

B.  REFRIGERANT  CONDENSATION  RESEARCH  AT  NPS 

In  order  to  examine  heat  transfer  characteristics  of  various  condenser  tube 
designs  and  alternative  refrigerants,  a  condenser/evaporator  test  platform  was 
constructed  at  the  Naval  Postgraduate  School.  Construction  was  begun  by  Zebrowski 
[Ref.  4]  and  completed  by  Mabrey  [Ref.  5]  in  1988.  The  test  platform,  shown  in 
Figure  1.1,  consisted  of  an  evaporator  and  a  condenser  section  with  an  associated 
ethylene  glycol/water  mixture  coolant  system. 

The  evaporator  section  was  designed  so  that  in  addition  to  supplying  the 
refrigerant  vapor  to  the  condenser,  it  could  also  be  utilized  to  conduct  bundle  boiling 
experiments.  The  condenser  section  consisted  mainly  of  four  instrumented  condenser 
tubes  and  a  series  of  auxiliary  condenser  coils  used  to  control  pressure  in  the 
apparatus.  Details  of  design,  construction  and  operation  of  the  condenser/evaporator 
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test  platform  are  provided  by  Zebrowski  [Ref.  4],  Mabrey  [Ref.  5]  and  Mazzone 
[Ref.  6]. 

Unlike  a  typical  vapor  compression  refrigeration  system  used  aboard  ships,  the 
test  platform  at  NPS  lacked  a  compressor  between  the  evaporator  and  the  condenser. 
This  compressor  is  normally  required  to  increase  the  pressure  (and  hence 
temperature)  of  the  refrigerant  vapor  sufficiently  enough  so  that  it  can  give  up  its 
latent  heat  to  the  ambient  sea  water,  the  primary  coolant  Because  there  was  no 
compressor  in  the  NFS  test  platform,  the  condenser  must  operate  at  the  saturation 
conditions  of  the  evaporator.  Therefore,  a  much  colder  coolant  a  refrigerated 
ethylene  glycol/water  mixture,  was  needed  to  condense  the  vapor,  leading  to  laminar 
flow  conditions  within  the  condenser  tubes.  This  in  turn  leads  to  poor  heat  transfer 
due  to  a  large  inside  thermal  resistance. 

The  primary  objective  of  the  condensation  experiments  of  Mabrey  [Ref.  5]  and 
Mazzone  [Ref.  6]  was  to  accurately  determine  outside  heat  transfer  coefficients. 
Briefly,  this  was  accomplished  by  first  determining  the  heat  transfer  rate  to  the 
coolant  by: 

4  -  rh  cp  AT  '  (1.1) 

where  m  is  the  coolant  mass  flow  rate  and  AT  the  coolant  temperature  rise. 
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THERMOCOUPLE  **  1  THERMOCOUPLE  **  2 


Figure  1.1  NPS  Evaporator/Condenser  Refrigeration  Apparatus. 
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The  overall  heat  transfer  coefficient,  could  then  be  determined  by: 


where  A,  is  the  outside  surface  area  of  the  tube  and  LMTD  is  the  log  mean 
temperature  difference,  defined  as: 


LMTD 


(1.3) 


The  overall  heat  transfer  coefficient  can  also  be  thought  of  in  terms  of  a  sum  of 
thermal  resistances.  The  overall  thermal  resistance  is  given  by: 


The  overall  thermal  resistance  can  be  measured  experimentally,  the  wall  resistance, 
Rw,  is  known,  and  the  inside  resistance  can  be  determined  from  a  chosen  correlation. 


Therefore,  the  outside  heat  transfer  coefficient,  h0,  can  be  calculated  by: 


There  is  one  major  drawback  to  this  method.  As  mentioned  earlier,  when 
laminar  flow  heat  transfer  conditions  exist  within  the  condenser  tube,  it  leads  to  poor 
heat  transfer  performance  and  therefore  represents  the  dominant  thermal  resistance 
within  the  condenser.  This  inside  resistance  then  becomes  the  controlling  factor  in 
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the  determination  of  the  outside  heat  transfer  coefficient  Indeed,  if  the  inside 
resistance  is  a  large  proportion  of  the  overall  resistance,  then  the  accuracy  of  the 
calculated  outside  heat  transfer  coefficient  becomes  uncertain.  Therefore,  the 
accuracy  of  the  outside  heat  transfer  coefficient  relies  heavily  on  an  accurate 
correlation  to  determine  h;. 

An  alternative  is  to  reduce  the  inside  controlling  resistance  by  augmenting  the 
inside  heat-transfer  coefficient  Augmentation  techniques  for  lowering  the  inside 
thermal  resistance  primarily  focus  on  repeatedly  bringing  fresh  fluid  from  the  coolant 
bulk  stream  to  the  heated  inside  surface.  This  results  in  breaking  up  the  thermal 
boundary  layer  and  thereby  reducing  the  thermal  resistance.  The  technique  utilized 
in  theNPS  condensation  experiments  to  break  up  this  thermal  boundaiy  layer  made 
use  of  two  types  of  inserts,  a  twisted  tape  insert  and  a  wire  mesh  (HEATEX)  insert 
These  inserts  and  their  use  will  be  discussed  in  more  detail  in  later  chapters. 

The  results  of  condensation  experiments  conducted  by  Mazzone  [Ref.  6]  show 
signs  of  inaccuracies  related  to  the  choice  and  accuracy  of  the  inside  heat  transfer 
coefficient  correlation.  For  example,  in  attempting  to  determine  the  outside  heat 
transfer  coefficient  for  a  smooth  tube  with  no  insert,  Mazzone  found  that  the 
correlation  used  to  determine  the  inside  resistance  yielded  values  that  were  greater 
than  the  actual  measured  overall  resistance!  Obviously,  no  sensible  values  of  h0  could 
be  determined.  Similar  problems  existed  (but  to  a  lesser  extent,  due  to  the  greatly 
reduced  inside  thermal  resistance)  when  using  various  types  of  inserts.  Mazzone 
concluded  that  "...  laminar  flow  correlations  for  twisted  tape,  HEATEX,  and  no 


insert  conditions  must  be  developed  which  have  better  accuracy  than  those  used  ... 
allowing  for  more  accurate  determination  of  outside  heat  transfer  coefficients." 
Therefore,  the  determination  of  accurate  inside  heat  transfer  coefficients  is 
paramount  in  the  final  outcome  of  the  refrigeration  condensation  project. 

C  OBJECTIVES 

The  main  focus  of  this  thesis  was  therefore  to  accurately  determine  inside  heat 
transfer  coefficients  in  support  of  refrigerant  condensation  experiments  being 
conducted  at  NFS. 

The  specific  objectives  of  this  study  were  to: 

1.  Design  and  build  a  test  apparatus  to  accurately  determine  inside  heat  transfer 
coefficients  for  the  same  tubes  and  inserts  used  in  the  condensation 
experiments. 

2.  Develop  correlations  for  inside  heat  transfer  coefficients  under  laminar  flow 
conditions  both  with  and  without  inserts. 

3.  Using  the  correlations  mentioned  above,  reprocess  existing  condensation  data 
to  determine  more  accurate  values  of  the  outside  heat  transfer  coefficients. 
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II.  LITERATURE  SURVEY 


A.  INTRODUCTION 

T  aminar  fully  developed  and  developing  forced  convection  flows  have  been 
analyzed  in  great  detail  for  various  boundary  conditions.  The  study  of  heat  transfer 
in  laminar  flow  through  a  closed  duct  was  first  made  by  Graetz  [Ref.  7]  in  1883  and 
later  by  Nusselt  [Ref.  8]  in  1910.  An  excellent  and  thorough  study  of  the  theoretical 
development  (with  both  analytical  and  experimental  solutions)  for  various  laminar 
forced  convection  internal  flow  boundary  conditions  is  given  by  Shah  and  London 
[Ref.  9].  - 

This  survey  will  be  concerned  with  laminar  flow  heat  transfer  in  fully  developed 
and  developing  flows  subjected  to  both  a  constant  wall  temperature  and  constant 
heat  flux  boundary  condition.  It  will  also  summarize  some  of  the  augmentation 
techniques  used  to  enhance  tubeside  laminar  flow  heat  transfer  and  discuss  the 
details  involved  when  using  twisted  tape  and  wire  mesh  inserts  to  enhance  inside  heat 
transfer. 
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B.  LAMINAR  FLOW  IN  CIRCULAR  DUCTS 


1.  Boundary  Layers  and  Entrance  Lengths 

* 

a.  Hydrodynamic  Boundary  Layer 

For  laminar  flow  in  a  circular  tube,  fluid  enters  the  tube  with  a 
uniform  velocity  as  -hown  in  Figure  2.1;  but  when  the  fluid  particles  make  contact 
with  the  tube  wall,  viscous  effects  cause  them  to  assume  a  zero  velocity.  These 
viscous  effects,  associated  with  shear  stresses  acting  in  planes  parallel  to  the  tube 
wall,  tend  to  retard  the  motion  of  particles  in  adjoining  fluid  layers,  causing  a 
momentum  boundary  layer  to  develop.  The  momentum  boundary  layer  thickness,  S, 
increases  until  its  growth  is  stopped  by  symmetry  at  the  tube  centerline.  The  viscous 
effects  then  extend  over  the  entire  cross  section  of  the  tube.  The  flow  is  then  said  to 
be  hydrodynamically  fully  developed.  At  this  point,  the  fluid  velocity  distribution  at 
a  given  cross  section  is  independent  of  the  axial  distance,  x  (i.e.,  du/dx=0)  [Refs.  9, 
10,  11,  12]. 

The  hydrodynamic  entry  length  for  laminar  flow  (ReD  £  2300)  is 
defined  as  that  region  where  the  hydrodynamic  boundary  layer  is  developing  and  may 
be  obtained  from  the  approximation  [Ref.  10]: 

-  0.05 ReD  (21) 

n  u 
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This  expression  shows  that  for  ReD=2300,  the  entrance  length  is  115  times  the  tube 
diameter,  a  considerable  length  for  most  experimental  facilities. 


Figure  2. 1  Laminar  Hydrodynamic  Boundary  Layer  Development  in  a  Circular  Tube 
(Courtesy  [Ref.  10]) 

b.  Thermal  Boundary  Layer 

When  a  fluid  with  a  uniform  temperature  distribution  enters  a  circular 
duct  with  a  wall  temperature  that  is  greater  than  the  uniform  fluid  temperature, 
temperature  gradients  in  the  fluid  develop  and  convective  heat  transfer  occurs,  as  shown 
in  Figure  2.2.  The  region  of  fluid  where  these  temperature  gradients  exist  is  called  the 
thermal  boundary  layer.  The  thermal  boundary  layer  thickness,  <5„  grows  until  a 
dimensionless  fluid  temperature  distribution  is  independent  of  x.  The  thermal  boundary 
layer  is  then  considered  fully  developed.  This  point  is  usually  expressed  in 
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Surftc*  condition 

T  T/m  Al  ** 


Figure  2.2  Thermal  Boundary  Layer  Development  in  a  Heated  Circular  Tut  ■.* 
(Courtesy  of  [Ref.  10]) 

terms  of  a  non-dimensional  temperature  profile  given  by  [R  1>  j: 

J_  ,  0  (2.2) 

*  [T,<x)-Tm(x) J 

where  T,  is  the  tube  surface  temperature,  T  is  the  local  fluid  temperature,  and  T„ 
is  the  mean  temperature  of  the  fluid  over  the  cross  section  of  the  tube.  The  exact 
shape  of  the  fully  developed  temperature  profile  is  dependent  upon  the  boundary 
conditions  imposed,  either  a  uniform  surface  temperature  or  constant  heat  flux. 

For  laminar  flow,  the  thermal  entry  length  may  be  expressed  as  the 
following  [Ref.  10]: 


-  0.05 ReDPr  (2'3) 

D 

It  should  be  noted  that  the  rate  of  development  of  the  velocity  and 
temperature  profiles  in  the  entry  region  depend  upon  the  fluid  Prandtl  number.  For 
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Pi  =  1,  both  velocity  and  temperature  profiles  develop  at  the  same  rate.  However, 
for  Pr  >  1,  the  hydrodynamic  boundary  layer  develops  more  rapidly  than  the  thermal 
boundary  layer  (Xf4h  <  For  very  large  Prandtl  number  (Pr  >  100),  is  so 
much  smaller  than  X^,  that  it  is  reasonable  to  assume  that  a  fully  developed  velocity 
profile  exists  throughout  the  thermal  entry  region.  [Refs.  9,  10,  11,  12] 

2.  Fully  Developed  Flow 

Fully  developed  laminar  flow  of  an  incompressible  fluid  through  a  circular 
duct  exists  when  both  the  thermal  and  velocity  profiles  are  fully  developed.  Two 
boundary  conditions  must  now  be  addressed,  that  of  a  uniform  wall  temperature  and 
constant  wall  heat  flux. 

For  these  conditions,  the  general  form  of  the  energy  equation  in  cylindrical 
coordinates  for  an  incompressible  fluid  having  constant  properties  (with  no  heat 
generation  and  neglecting  body  forces,  pressure  gradients  and  viscous  dissipation), 
is  given  by  [Ref.  11]: 


ff  ar  ^  w  st  _  a  d  (rn 

dx  dr  r  dr  dr  , 


(2.4) 


It  can  be  shown  that  for  these  two  cases,  the  energy  equation  reduces  to 
the  following  differential  equations, 


Constant  Wall  Heat  Flux: 


1 

r  dr 


(,*D 

Srj 


u 

a 


dh. 

dx 


\ 


(2.5) 


72 


(2.6) 


Uniform  Wall  Temperature: 


-  u 

'  TrT  ' 

(  \ 

<nm 

a  1 

T-T 

1  m) 

k  , 

For  the  constant  wall  heat  flux  boundary  condition,  Kays  [Ref.  11]  solves 
the  differential  equation  (equation  2.5)  for  a  fully  developed  velocity  and 
temperature  profile  and  expresses  the  Nusselt  number  as: 

Nu  -  —  »  4.364  (2.7) 

k 

Note  that  for  this  specified  case,  the  heat  transfer  coefficient,  h,  depends 
only  on  the  thermal  conductivity,  k,  and  the  tube  diameter,  D. 

A  solution  to  the  uniform  wall  temperature  case  is  more  complex,  and 
involves  evaluating  the  infinite  series  solution  to  the  differential  equation  (equation 
2.6)  for  the  hydrodynamically  developed,  thermally  developing  flow  situation.  In 
evaluating  the  series  solution,  the  resulting  asymptotic  Nusselt  number  is  given  by 
[Ref.  11]: 

Nu  *  3.657  (2-8) 

This  result  is  16%  less  than  the  solution  for  the  constant  wall  heat  flux  boundary 
condition.  The  asymptotic  values  in  equations  2.7  and  2.8  are  shown  in  Figure  2.3. 
Note  that  the  inverse  of  the  Graetz  number  is  equivalent  to  the  dimensionless  axial 
distance  x*. 
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3.  Entry  Length  Problem 


The  solution  to  the  energy  equation  for  the  entry  region  is  far  more 
difficult  to  obtain  because  velocity  and  temperature  vary  both  axially  and  radially. 
There  are  two  entry  length  problems  to  be  considered,  a  combined  entry  length, 
where  both  temperature  and  velocity  profiles  develop  simultaneously,  and  the 
thermal  entry  length  problem,  where  the  thermal  conditions  develop  in  the  presence 
of  a  fully  developed  velocity  profile.  This  latter  case  would  exist  if  the  place  at  which 
heat  transfer  begins  is  preceded  by  an  unheated  starting  length,  a  situation  that  exists 
in  most  experimental  facilities.  It  should  be  noted  that  for  Pr  <  <  1  (i.e.,  liquid 
metals),  the  thermal  boundary  layer  would  develop  faster  than  the  hydrodynamic 
boundary  layer.  However,  this  case  has  little  practical  application  and  is  therefore 
rarely  considered  in  the  literature. 

Solutions  have  been  obtained  for  both  entry  length  conditions  mentioned 
above  [Refs.  9,  10,  11,  12]  and  these  results  are  also  presented  in  Figure  2.3.  The 
Nusselt  number  is  plotted  against  a  dimensionless  axial  coordinate  x*  (the  inverse  of 
the  Graetz  number)  given  by: 


x ID 
Refr 


(2.9) 


Fully  developed  conditions  are  reached  at  x*  —  0.05  where  the  solutions  for  the 
thermal  entry  length  converge  with  those  for  the  combined  entry  length.  For  x*  > 
0.05,  the  two  types  of  entry  length  asymptotically  reach  the  fully  developed  solutions 
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of  4.36  and  3.66  for  the  constant  wall  heat  flux  and  uniform  wall  temperature 
boundary  conditions  respectively. 


4.36 

3£6 


Figure  2.3  Local  Nusselt  Number  Obtained  from  Entry  Length  Solutions  for 
Laminar  Flow  in  a  Circular  Tube  (Courtesy  [Ref.  10]) 

a.  Thermal  Entry  Region  -  Uniform  Wall  Temperature 

The  thermal  entry  length  problem,  first  studied  by  Graetz  in  1883  and 
later  by  Nusselt  in  1910,  is  more  commonly  known  as  the  Graetz  Problem.  They 
considered  an  incompressible  fluid  with  constant  properties  flowing  through  a 
circular  tube  having  a  fully  developed  laminar  velocity  profile  and  a  developing 
laminar  temperature  profile  subjected  to  a  uniform  wall  temperature.  Shah  and 
London  [Ref.  9]  present  a  numerically  derived  infinite  series  solution  to  the  Graetz 
Problem. 
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A  more  useful  correlation  presented  by  Incropera  and  Dewitt  [Ref. 
10]  for  the  uniform  wall  temperature  thermal  entry  length  problem  was  developed 
by  Hausen  and  is  of  the  form: 


„  „  _  0.0668  (D/L)ReDPr 

Nu  *  3.66  + - 

1  +Q.Q4[(D/L)Re0Pr]vi 


(2.10) 


This  correlation  was  used  by  Ma2zone  [Ref.  6]  in  his  refrigerant  condensation 
experiments  to  predict  the  average  inside  heat  transfer  coefficient  for  a  smooth  tube 
without  any  type  of  insert 

Shah  and  London  [Ref.  9]  present  a  more  complete  set  of  correlations 
derived  from  experimental  data  for  heat  exchanger  tubes  in  laminar  flow  for  local 
(Nu,)  and  mean  (Nu,,)  Nusselt  numbers  based  on  various  x*  ranges.  The  correlations 
for  local  Nusselt  number  are  given  by: 


Nux  *  1.077(x*)*w-0.7  forx'S  0.01 

Nux  =  3.657  +  6.874(103x*)'a488e  'S7J*'  for  x*  >  0.01 

and  for  mean  Nusselt  numbers  by: 

iVu(lt  =  1.615(x,)-‘/,-0.7  for  x*  s  0.005 

Num  =  X.6l5(xym -02  for  0.005<x*  <  0.03 

Num  =  3.657  *  H499  for  x  *  >  0.03 

x* 


(2.11) 


(2.12) 


The  reader  is  reminded  that  both  the  Hausen  correlation  and  the 
Shah  and  London  correlations  presented  above  are  for  a  thermally  developing  flow 
with  a  fully  developed  velocity  profile  subjected  to  a  uniform  wall  temperature. 
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b.  Thermal  Entry  Region  -  Constant  Wall  Heat  Flux 

For  the  thermal  entry  length  problem  under  a  constant  wall  heat  flux 
boundary  condition,  Shah  and  London  [Ref.  9]  present  several  numerical  solutions 
but  recommend  the  following  approximate  equations  for  the  local  Nusselt  number: 

Nux  =  lJ02(x*)'1/3-l  forx'  jS  5X10'3 

Nux  -  l-302(;c*)‘‘/3-0-5  for  5xl0-3  £  x’  S  1.5  xlO'3  (2.13) 

Nux  =  4.364 +8.68(103r*)‘<X5“e  ‘4U*  for  x*  S  1.5 x  10° 

Shah  and  London  also  recommend  the  following  approximate  equations  for  the  mean 
Nusselt  number: 


Num  =  1.953(x*)"w  for  xm  s  0.03 

Nu  -  4.364  >  0  0722  forx'  >  0.03 

xm 


(2.14) 


The  above  two  sets  of  equations  are  for  a  thermally  developing, 
hydrodynamically  developed  flow  subjected  to  a  constant  wall  heat  flux. 


c.  Combined  Entry  Length 

As  discussed  earlier,  the  rate  at  which  the  temperature  and  velocity 
profiles  develop  is  dependent  upon  the  Prandtl  number.  At  medium  to  high  Prandtl 
numbers  (Pr  >  100),  the  velocity  profile  develops  much  more  rapidly  than  the 
temperature  profile  and  a  fully  developed  velocity  profile  assumption  is  valid;  for 
such  conditions  the  equations  given  above  should  be  used. 

For  1  <  Pr  <  100,  however,  a  fully  developed  velocity  profile  is  no 
longer  a  valid  assumption  and  a  combined  entry  length  solution  must  be  developed. 
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Again,  Shah  and  London  [Ref.  9]  present  a  comprehensive  study  of  various  solutions 
to  the  combined  entry  length  problem,  including  two  solutions  for  local  Nusselt 
numbers  developed  by  Churchill  and  Ozoe  [Ref  13],  the  first,  for  a  uniform  wall 
temperature: 


Nil, +1.7 


1+ 


5.357(1  +[(388/w)x  •]-*9p8 
ir/(284x*) 

[  {l+(Fr/0.0468)w}w(l+[(388/ir)x*]-M)3/4  J 


V3 


v» 


(2.15) 


and  the  second,  for  a  constant  wall  heat  flux* 


Afo,+1 


1  + 


5.364(1  +((220/tr)x*]'uv9)3/10 

ir/(U5.2x  *) 

[  (1  +(/V/0.0207)M)w(l  +[(220/ir)r  •]-1<w}3'3 


5/3 


3/10 


(2.16) 


From  the  preceding  paragraphs,  it  can  be  seen  that  the  subject  of 
laminar  flow  forced  convection  in  circular  ducts  has  been  studied  in  great  detail 
yielding  both  analytical  and  experimental  results.  The  reader  is  referred  to  Shah  and 
London  [Ref.  9]  and  Kays  and  Crawford  [Ref.  11]  for  a  more  complete  and  thorough 
review  of  this  subject 
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C.  AUGMENTATION  TECHNIQUES  TO  ENHANCE  TUBESIDE  HEAT 
TRANSFER 

1.  General  Introduction 

As  discussed  earlier,  when  laminar  flow  heat  transfer  occurs  within  tubes, 
it  usually  represents  the  dominant  thermal  resistance  in  a  heat  exchanger.  These 
large  resistances  result  in  low  inside  heat  transfer  coefficients.  In  recent  years,  the 
requirement  for  more  efficient  heat  exchangers  has  stimulated  considerable  interest 
in  techniques  to  augment  or  enhance  tubeside  heat  transfer. 

The  various  methods  to  enhance  tubeside  heat  transfer  can  be  classified 
into  two  categories,  passive  and  active  techniques.  Passive  techniques  involve  no 
external  energy  (except  increased  pump  work)  to  produce  the  necessary 
enhancement  Passive  techniques  are  divided  into  the  following  methods: 


1.  Internally  finned  tubes  (longitudinal  or  spiral  fins)  designed  to  increase  the 
inside  heat  transfer  surface  area. 

2.  Surface  roughness  techniques  aimed  at  agitating  flow  rather  than  increasing 
heat  transfer  surface  area. 

3.  Swirl  flow  techniques  which  involve  the  use  of  inlet  vortex  generators, 
periodically  spread  propellers  or  twisted  tape  inserts  designed  to  cause 
enhancement  by  creating  a  rotating  and/or  secondary  flow. 

4.  Displaced  promoters  which  alter  the  fluid  flow  near  the  surface  by  the  use  of 
inserts  such  as  wire  mesh,  static  mixer  elements,  rings  or  disks. 

5.  Coolant  additives  designed  to  increase  the  conductivity  of  the  coolant 

6.  Compound  techniques  that  involve  more  than  one  of  the  above  mentioned 
methods. 
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Active  techniques  utilize  an  external  energy  source  to  promote  the 
required  augmentation.  These  techniques  include:  (1)  mechanical  aids,  (2)  heated 
surface  vibration,  (3)  fluid  pulsation,  (4)  electrostatic  fields,  and  (5)  suction  or 
injection. 

All  of  the  above  techniques  are  aimed  at  creating  turbulence  in  the 
relatively  slow  moving  or  stagnant  boundary  layer  fluid  by  displacing  this  fluid  from 
the  heateu  surface  and  mixing  it  with  the  bulk  fluid.  Essential  in  all  of  these 
techniques  is  a  balance  between  increased  heat  transfer  rate  and  increased  pumping 
power  requirements/costs.  Bergles  and  Joshi  [Ref.  14]  present  an  excellent  review 
of  augmentation  techniques  and  compare  heat  transfer  data  taken  at  low  Reynolds 
number  flow  subjected  to  both  uniforarwafl- temperature  and  constant  wall  heat  flux. 

Of  particular  interest  to  this  thesis  are  two  passive  techniques,  the  swirl 
flow  technique  using  twisted  tape  inserts  and  the  displaced  promoter  technique  using 
wire  mesh  inserts.  These  two  methods  of  enhancement  of  tubeside  heat  transfer  are 
discussed  in  more  detail  in  the  following  sections. 

2.  Heat  Transfer  Enhancement  with  Twisted  Tape  Elements 

Among  the  most  common  swirl  flow  augmentation  techniques  is  the  use 
of  twisted  tape  inserts.  The  tapes  usually  consist  of  a  thin  stainless  steel,  brass  or 
copper  strip  which  is  the  width  of  the  inside  diameter  of  the  tube.  The  tapes  are 
twisted  by  clamping  one  end  to  an  overhead,  attaching  weights  to  the  other  end  and 
then  twisting  to  the  desired  specifications  [Ref.  15].  Enhancement  occurs  not  only 
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due  to  the  increased  path  length  of  flow,  but  also  due  to  secondary  flow  effects  and 
fin  effects  caused  by  tape  contact  with  the  inside  tube  wall. 

Analytical  studies  of  twisted  tape  inserts  by  Date  and  Singham  [Ref.  16] 
suggest  that  swirl  flows  improve  heat  transfer  by  as  much  as  a  factor  of  70;  however, 
no  experimental  data  have  confirmed  their  unrealistically  large  predictions.  The  first 
experimental  work  of  swirl  flow  using  twisted  tapes  was  conducted  by  Hong  and 
Bergles  [Ref.  15].  Their  primary  objective  was  to  develop  experimentally  based 
correlations  for  predicting  heat  transfer  coefficients  for  laminar  flow  of  water  and 
ethylene  glycol  in  an  electrically  heated  tube  (heated  length  of  1.22  m,  inside 
diameter  of  10.2  mm)  with  two  sizes  of  twisted  tape  inserts.  They  showed  that  the 
Nusselt  number  was  independent  of  axial  location,  .suggesting  that  fully  developed 
conditions  existed.  Therefore,  the  mean  Nusselt  number  should  take  the  form: 


Num  -  m,Pr.y). 

Re,  is  the  Reynolds  number  for  the  swirl  flow  given  by: 


Re. 


4m 

irp(D -45) 


and  5  is  the  tape  thickness.  The  twist  ratio,  y,  is  defined  as: 


(2.17) 


(2.18) 


(2.19) 
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where  H  is  the  pitch  for  a  180°  twist  of  the  tape.  Hong  and  Bergles  then  combined 
Re,  and  y  as  a  single  parameter  and  came  up  with  the  following  correlation  for 
values  of  Re^y  >  10: 

Num  =  0.383 Pr^iReJyf632  (2-20) 

For  values  of  Re,/y  <  10,  the  following  correlation  was  obtained: 

Num  *  5.172(1  5.484  x  IQ-' Pr^Re/y)1*)0-5  (2-21) 

It  should  be  noted  that  all  fluid  properties  in  equations  2.20  and  2.21 
should  be  evaluated  at  the  bulk  fluid  temperature.  One  additional  important 
observation  was  made  by  Hong  and  Bergles  regarding  tube  wall  temperature.  They 
found  a  significant  variation  of  temperature  aronemrial  location  with  two  different 
tape  orientations  at  that  position  of  90°  and  120°,  suggesting  that  the  temperature 
profile  was  related  to  the  tape  orientation.  Hong  and  Bergles  concluded  that  for  a 
fully  developed  flow,  inside  heat  transfer  coefficients  depend  on  twist  ratio  and  could 
be  improved  by  a  factor  of  two  or  three  over  an  empty  tube.  These  conclusions  were 
confirmed  by  later  studies  of  Saha  eL  al.  [Ref.  17]  and  Marner  and  Bergles  [Ref. 
18]. 

3.  Heat  Transfer  Enhancement  with  Wire  Mesh  (HEATEX)  Elements 

The  use  of  wire  mesh  inserts  to  promote  inside  heat  transfer  falls  into  the 
displaced  promoter  technique  mentioned  earlier.  One  such  wire  mesh  insert  was 
commercially  developed  by  CAL-GAVIN  Limited  (UK)  and  called  "HEATEX." 
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This  HEATEX  insert  consists  of  a  central  wire  core  on  which  a  series  of 


wire  loops  or  petals  are  attached.  The  petals  are  attached  to  the  core  at  an  inclined 
angle  that  face  the  oncoming  fluid  flow.  Ideally,  each  petal  makes  contact  with  the 
tube  wall  in  the  form  of  an  arc,  thereby  maximizing  disturbance  of  fluid  along  the 
tube  wall.  In  addition  to  the  boundary  layer  disturbance,  radial  mixing  occurs  caused 
by  repeated  redistribution  and  mixing  of  fluid  through  the  mesh  of  loops  near  the 
tube  axis  [Ref.  19].  Various  wire  loop  densities  and  nominal  diameters  are 
available,  depending  on  the  application. 

Initial  studies  by  Gough  et  ai.  [Ref.  19]  with  HEATEX  inserts  suggest 
enhancements  of  a  factor  greater  than  ten  times  that  of  a  similar  tube  with  no  insert 
Additional  experiments  were  conducted  by  Oliver  and  Aldington  [Ref.  20]  in  which 
a  correlation  was  determined: 

Nu  -  0.232ReOJ4/Vw6  (2-22) 

It  must  be  remembered,  however,  that  this  correlation  is  for  a  specific  loop  density 
and  nominal  diameter  HEATEX  element 

The  two  types  of  HEATEX  mixing  elements  supplied  by  CAL-GAVIN 
(nominal  diameters  of  13.26  and  10.16  mm)  did  not  match  the  loop  densities  of 
previous  studies;  therefore  new  correlations  needed  to  be  developed.  Mazzone  [Ref. 
6]  developed  correlations  for  these  inserts  from  limited  data  supplied  by  CAL- 
GAVIN.  For  the  13.26  mm  and  10.16  mm  elements,  these  correlations  were 
respectively: 
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Nu  =  0.22 6Re0xaPr0M  (13.26  mm)  i2-23) 

and 

Nu  =  0.063 Re0nPr0M  (10.16mm)  (Z24) 

Since  these  correlations  were  developed  from  a  very  limited  amount  of 
data  taken  under  conditions  that  did  not  exactly  match  that  of  Mazzone’s 
condensation  experiments,  they  were  always  subject  to  inaccuracy.  Therefore,  the 
need  for  more  accurate  correlations  under  the  exact  conditions  used  in  the 
condensation  experiments  are  required. 
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III.  EXPERIMENTAL  APPARATUS 


A.  SYSTEM  OVERVIEW 

The  test  apparatus  utilized  is  shown  in  Figures  3.1  and  3.2.  The  apparatus  is  a 
closed  loop,  low  pressure  system  utilizing  an  ethylene  glycol/water  mixture  as  the 
coolant  The  coolant  is  circulated  by  a  centrifugal  pump  from  a  supply  tank,  through 
a  reducing  manifold  and  flowmeter  into  the  test  section.  Coolant  is  then  returned  to 
the  supply  tank  via  a  return  header. 

The  test  section  consisted  of  an  inlet  and  outlet  thermocouple  well/mixing 
chamber  to  measure  inlet  and  outlet  coolant  temperatures  and  the  test  tube  itself. 
The  tube  was  instrumented  with  twelve  thermocouples  (embedded  in  the  tube  wall) 
to  determine  inside  wall  temperatures.  A  constant  heat  flux  boundary  condition  was 
achieved  using  a  series  of  flexible  electrical  heating  elements  wrapped  around  the 
tube.  Monitoring  the  system  as  well  as  data  acquisition  and  data  processing  were 
provided  by  a  Hewlett-Packard  microcomputer  and  Data  Acquisition  System. 
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EXPERIMENTAL  APPARATUS 


Figure  3.1  Primary  and  Secondary  Coolant  Systems. 
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CHILLER 


DATA  ACQUISITION  &  CONTROL  SYSTEM 


Figure  3.2  Data  Acquisition  and  Control  Systems. 


27 


INLET  T/C 


B.  TUBES  AND  INSERTS  TESTED 
1.  Tubes  Tested 

Three  types  of  tubes  were  utilized  in  the  experiments,  a  copper  smooth 
tube,  a  copper/nickel  externally  finned  tube,  and  a  copper/nickel  corrugated  or  roped 
tube.  These  tubes  correspond  to  the  tubes  utilized  in  the  refrigerant  condensation 
experiments  of  Mabrey  [Ref.  5]  and  Mazzone  [Ref.  6]. 

The  smooth  copper  tube  provided  a  baseline  set  of  experimental  data 
which  were  then  compared  with  previous  experimental  and  theoretical  work  for  the 
inside  heat  transfer  coefficient  conducted  by  Shah  and  London  [Ref.  9]  for  no  insert, 
Oliver  and  Aldington  [Ref.  20]  for  the  HEATEX  insert,  and  Hong  and  Bergles  [Ref. 
15]  for  the  twisted  tape  insert 

The  Cu/Ni  finned  tube  provided  a  smaller  inside  diameter  to  be  used  in 
these  comparisons.  Note  that  the  fins  play  little  part  in  the  experiments.  The  wall 
thermocouples  were  buried  in  the  tube  wall  and  were  used  to  measure  the  inside  wall 
temperature,  which  had  a  smooth  bore.  Consequently,  for  the  purposes  of  these 
experiments,  the  tube  behaved  like  a  smooth  copper  tube  with  a  different  internal 
diameter. 

The  final  tube  tested  was  a  Cu/Ni  corrugated  (or  roped)  tube, 
commercially  referred  to  as  a  "Korodense"  tube.  The  tube  was  manufactured  and 
provided  by  Wolverine  Tube  Company  of  Decatur,  AJ.  All  tube  dimensions  and 
specifications  were  taken  from  the  Wolverine  Technical  Bulletin  No.  4020  [Ref. 
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21].  Figure  3.3  is  a  schematic  of  the  Korodense  tube.  Tube  specifications  for  all 
three  tubes  are  listed  in  Table  E.1  of  Appendix  E. 


Figure  3.3  Wolverine  Korodense  Tube  (Courtesy  [Ref.  21]) 

2.  Insert  Elements  Tested 

Two  types  of  insert  mixing  elements  were  utilized  in  this  investigation,  a 
twisted  tape  insert  and  a  wire  mesh  insert  commercially  referred  to  as  a  "HEATEX" 
tube  insert  In  both  cases,  two  sizes  were  required  due  to  the  two  different  tube 
inside  diameters  used. 

The  twisted  tape  inserts  were  made  of  a  strip  of  brass  0.559  mm  thick  and 
either  10.16  mm  or  13.26  mm  wide.  The  tapes  were  manufactured  in  the  NPS 
Mechanical  Engineering  machine  shop  by  clamping  them  to  the  overhead,  attaching 
a  weight  to  the  opposite  end,  and  twisting  them  to  the  desired  twist  ratio.  A  twist 
ratio  of  3  and  4  corresponding  to  the  two  respective  tube  inside  diameters  were 
manufactured.  Once  installed  in  the  tube,  the  edges  of  the  tape  should  ideally  have 
been  in  contact  with  the  tube  wall. 
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The  second  type  of  insert  used  was  a  commercially  available  product 
referred  to  as  HEATEX  radial  mixing  elements.  These  elements  were  supplied  by 
Cal  Gavin  of  Birmingham,  U.K.  They  are  manufactured  from  stainless  steel  and 
consist  of  a  central  wire  core  onto  which  a  series  of  wire  loops  or  petals  are  attached. 
Each  petal  was  inclined  at  an  angle  facing  the  oncoming  flow.  The  outside  diameter 
of  the  elements  were  slightly  oversized  for  the  tube,  ensuring  that  close  contact  with 
the  tube  wall  in  the  form  of  an  arc,  rather  than  a  point,  was  achieved.  As  with  the 
twisted  tape  elements,  different  diameter  HEATEX  elements  were  needed,  13.26 
mm  used  in  the  smooth  and  Korodense  tubes  and  the  other  10.15  mm  for  the  Cu/Ni 
finned  tube.  Figure  3.4  shows  both  the  twisted  tape  and  HEATEX  inserts. 


Figure  3.4  Twisted  Tape  and  HEATEX  Inserts 


30 


C  COOLANT  SYSTEM 


1.  Primary  Coolant  System 

An  ethylene  glycol/water  mixture  (54/46%  by  weight)  was  stored  in  a  1.5 
m3  insulated  supply  tank.  The  mixture  was  circulated  through  the  primary  system  by 
a  1  Vi  hp,  constant  speed,  centrifugal  pump.  The  coolant  was  discharged  from  the 
pump  via  a  76  mm  diameter  PVC  pipe  to  a  manifold  where  it  was  reduced  to  a  15.9 
mm  diameter  Tygon  flexible  tube.  A  rotameter  type  flowmeter  with  a  ball  valve  at 
its  entrance  was  used  to  control  the  flow  rate  through  the  system.  Flowmeter 
calibration  results  are  discussed  in  Appendix  A. 

Flow  entered  the  test  section  through  a  90°  bend  which  housed  the 
thermocouple  well  and  exited  through  a  mixing  chamber/thermocouple  well 
combination  (see  Figure  3.1).  Prior  to  entering  the  exit  mixing  chamber,  the  coolant 
passed  through  a  10  cm  length  of  flexible  tubing  containing  a  high  density  wire  mesh. 
This  wire  mesh,  in  conjunction  with  the  mixing  chamber,  served  to  minimize  radial 
temperature  variations  and  provide  a  uniform  outlet  coolant  temperature.  The  flow 
then  exited  the  test  section  through  more  flexible  tubing  to  a  manifold  and  return 
PVC  header,  suspended  above  the  apparatus,  and  back  to  the  coolant  supply  tank. 

2.  Secondary  Coolant  System 

In  order  to  achieve  a  wide  range  of  Reynolds  number,  coolant  inlet 
temperatures  were  varied  from  -10°C  to  +20°C.  This  was  accomplished  by  a 
secondary  coolant  system  comprised  of  a  28  kW  (8  ton)  external  refrigeration  system. 
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A  30  gpm,  .75  hp  centrifugal  pump  continuously  circulated  coolant  between  the 
supply  tank  and  a  chiller  barrel  where  it  was  cooled  to  the  desired  temperature. 
Temperature  was  maintained  by  a  thermostatically  controlled  solenoid  valve  at  the 
discharge  side  of  the  chiller  barrel.  This  system  allowed  for  a  range  of  coolant 
temperatures  between  -20°C  and  ambient  to  be  selected  while  maintaining  the 
coolant  inlet  temperature  approximately  constant  (±2°C)  during  each  run. 

D.  INSTRUMENTATION  AND  TEST  SECTION 
1.  Instrumentation 

All  temperature  measurements  were  made  using  copper/constantan  type-T 
teflon  coated  thermocouples  and  were  previously  calibrated  by  Mabrey  [Ref.  5].  The 
inlet  coolant  temperature  was  measured  by  a  single  thermocouple  while  the  outlet 
coolant  temperature  was  measured  by  two  thermocouples  located  immediately  after 
the  mixing  chamber. 

Each  test  tube  was  fitted  with  twelve  wall  thermocouples  to  determine  the 
average  inside  wall  temperature.  Thermocouples  were  placed  at  four  different 
longitudinal  positions  along  the  axis  of  the  tube.  Three  thermocouples  at  each 
location  were  spread  evenly  around  the  tube  at  120°C  intervals.  Figure  3.5  shows  a 
typical  tube  with  the  exact  longitudinal  and  nodal  thermocouple  locations.  Actual 
longitudinal  locations  varied  slightly  for  each  tube  and  are  given  in  Table  E.l  as  x/D; 
where  x  is  the  length  from  the  start  of  the  heated  length  to  the  thermocouple 
location  and  D;  in  the  inside  diameter  of  the  tube.  The  two  extreme  longitudinal 
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locations  of  the  thermocouples  corresponded  to  the  beginning  and  end  of  the 
"effective"  condensing  length  of  the  tubes  being  used  in  the  condensation  experiments 
of  Mabrey  [Ref.  5]  and  Mazzone  [Ref.  6]. 


INSTRUMENTED  TUBE 


inlet  oo 


POS  i 
I5*2cm 


7 


POS  2 
55.9cm 


POS  3 
96.5cm 


POS  4 
137.2cm 


I*  1.52m 
1.42m 

L/D  =  105.6 


Do  =  15. 88mm 
Di  *  13.26mm 
Dt/c-14.36mm 


Figure  3.5  Instrumented  Tube 


Di  OUTLET 
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2.  Thermocouple  Attachment  Method 


In  order  to  attach  the  thermocouples  to  the  wall  of  the  tube,  short 
longitudinal  grooves  were  cut  at  each  location.  These  longitudinal  grooves  were  cut 
to  a  depth  of  Vi  the  tube  wall  thickness.  It  was  decided  that  a  thermocouple 
embedded  within  the  tube  wall  vice  being  attached  to  the  surface  would  lead  to  a 
more  accurate  tube  wall  temperature  and  minimize  the  effect  of  direct  contact 
heating  from  the  heating  elements.  Thermocouples  were  then  placed  in  these  grooves 
and  either  soldered  in  place  or  tack  welded  and  covered  with  a  liquid  metal  epoxy. 
Different  methods  for  attaching  and  covering  the  thermocouples  were  used  for  each 
tube  based  on  the  tube  material  and  configuration  (discussed  below). 

In  the  case  of  the  smooth  copper  tube,  the  thermocouples  were  silver 
soldered  in  place  and  the  remainder  of  the  groove  then  covered  with  silver  solder. 
Since  the  tubes  were  copper,  this  was  relatively  easy  to  accomplish.  Care  was  taken 
to  ensure  that  the  thermocouple  remained  in  contact  with  the  tube  wall  and  that  no 
air  pockets  were  created  (which  would  cause  an  additional  large  resistance  to  heat 
flow).  Also,  great  care  was  taken  to  prevent  any  melting  of  the  teflon  insulation  past 
the  thermocouple  joint  This  condition  could  have  caused  an  electrical  contact 
between  the  two  wires  of  the  thermocouple  prior  to  the  tube  wall,  leading  to  an 
unrealistically  high  observed  wall  temperature. 

A  different  method  was  required  to  attach  the  thermocouples  to  the 
copper/nickel  finned  tube  due  not  only  to  the  external  fins  but  also  the  10%  nickel 
content  Initially,  a  small  diameter  hole  was  machined  into  the  tube  to  minimize  fin 
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damage  (rather  than  cutting  a  longitudinal  groove).  Since  the  thermocouple  was  to 
be  embedded  half  way  into  the  tube  wall,  this  hole  proved  to  be  relatively  deep. 
Also,  because  of  the  nickel  content  of  these  tubes,  the  preheat  temperature  required 
to  ensure  that  the  silver  solder  would  adhere  to  the  tube  wall  was  much  higher. 
These  higher  required  temperatures  and  the  hole  depth  caused  a  significant  amount 
of  the  thermocouple  teflon  insulation  to  be  melted,  allowing  for  potential  contact 
between  the  thermocouple  wires  at  some  point  other  than  the  tube  wall.  This 
technique  was  therefore  abandoned  in  favor  of  the  longitudinal  groove  method. 
However,  thermocouple  insulation  was  again  being  destroyed  due  to  the  higher 
required  preheat  temperatures.  An  alternative  solder,  soft  solder,  with  lower  preheat 
and  melting  temperatures,  was  therefore  used.  The  final  combination  of  longitudinal 
grooves  and  soft  solder  was  successful. 

For  the  copper/nickel  Korodense  tube,  the  groove  and  soft  solder  method 
was  again  tried.  Initial  tests  in  the  apparatus  showed  wall  temperatures  at  one 
longitudinal  position  varying  by  as  much  as  15°C.  On  closer  examination  of  the 
thermocouple  connections  (using  a  magnifying  glass),  it  was  revealed  that  insulation 
damage  was  evident  on  over  half  of  the  thermocouples.  Several  more  attempts  at  the 
groove  and  soft  solder  method  were  tried,  all  leading  to  the  same  results. 

A  third  method  utilizing  a  tack  welding  machine  (acquired  and  adapted 
for  our  use)  was  tried.  Thermocouples  were  tack  welded  into  the  bottom  of  the 
groove  resulting  in  no  damage  to  the  insulation  and  a  liquid  steel  epoxy  was  used  to 
fill  the  remainder  of  the  groove,  providing  additional  strength  and  protection  to  the 
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connection.  This  final  method  proved  to  be  the  most  successful  and  is  recommended 
for  use  in  any  follow-up  experiments. 

3.  Test  Section 

In  order  to  achieve  a  constant  heat  flux,  each  tube  was  tightly  wrapped 
with  six  silicone  rubber  flexible  electrical  heating  elements  rated  at  240  V,  200  W 
each.  Each  strip  was  2.5  cm  wide  and  50.6  cm  long.  The  heating  elements  were 
connected  in  parallel  to  a  240  Vac  power  supply  via  a  plug  panel  and  power 
switchboard.  Power  to  the  six  heaters  was  controlled  by  a  STACO,  240  Vac,  23.5  kVa 
rheostat  (see  Figure  3.2). 

To  minimize  heat  loss  from  the  heater  to  ambient,  the  tube  was  wrapped 
in  a  20  mm  thick  neoprene  type  insulation.  This  type  of  insulation  proved  to  be 
extremely  effective  giving  a  negligible  heat  loss.  It  also  provided  insulation  against 
heat  influx  from  ambient  at  the  colder  coolant  inlet  temperatures.  Actual  heat  loss 
and  heat  influx  experiments  and  their  results  are  discussed  in  more  detail  in 
Appendices  B  and  C.  All  other  tubing  and  PVC  piping  was  wrapped  in  similar  type 
insulation. 

The  whole  test  section  was  then  placed  in  a  restraining  stand  which 
minimized  outside  vibration  effects  and  ensured  that  each  tube  remained  horizontal 
and  in  the  same  relative  position  (to  the  remainder  of  the  apparatus)  for  each 
experiment 


E.  DATA  ACQUISITION  AND  CONTROL  SYSTEM 

A  Hewlett-Packard  9125  microcomputer  (MC)  was  used  to  control  a  Hewlett- 
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Packard  3497A  Automatic  Data  Acquisition  Unit  (DAU)  (see  Figure  3.2).  The  DAU 
read  the  output  from  the  twelve  tube  wall  thermocouples  and  three  coolant 
thermocouples  (one  inlet  and  two  outlet).  Additional  channels  were  assigned  to  sense 
heater  input  voltage  and  current  at  the  plug  panel.  The  voltage  was  converted  to 
heater  power  using  the  heater  resistance  and  the  calibration  discussed  in  Appendix 
D. 

Thermocouple  readings  were  made  in  microvolts  which  were  converted  to 
temperatures  (along  with  the  heat  flux  calculations)  using  the  data  reduction 
program.  The  data  reduction  program,  DRPSING,  will  be  discussed  in  more  detail 
in  a  later  chapter.  All  temperatures,  power,  and  heat  flux  measurements  were 
monitored  using  the  MC  and  DAU.  DAU  channel  assignments  are  listed  in  Table 
3.1. 
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TABLE  3.1  DATA  ACQUISITION  UNIT  CHANNEL  ASSIGNMENTS 


CHANNEL  NO. 

MEASURED  QUANTITY 

1 

Thermocouple  Pos  1 

2 

Thermocouple  Pos  2 

3 

Thermocouple  Pos  3 

4 

Thermocouple  Pos  4 

5 

Thermocouple  Pos  5 

6 

Thermocouple  Pos  6 

7 

Thermocouple  Pos  7 

8 

Thermocouple  Pos  8 

9 

Thermocouple  Pos  9 

10 

Thermocouple  Pos  10 

11 

Thermocouple  Pos  11 

12 

Thermocouple  Pos  12 

13 

Tube  Inlet 

14,  15 

Tube  Outlet 

25 

Input  Current 

29 

Input  Voltage 
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rv.  EXPERIMENTAL  PROCEDURES 


A.  CALIBRATION  EXPERIMENTS 

Several  calibration  experiments  were  required  prior  to  performing  the  main 
experiments.  Figure  4.1  depicts  a  typical  tube  being  tested.  The  primary  heat  source 
is  from  the  heating  elements  (q^).  However,  to  complete  the  heat  balance  on  the 
system,  several  other  heat  sources  and  heat  sinks  needed  to  be  accounted  for.  The 
first  experiment  conducted  was  to  determine  the  amount  of  heat  being  lost  from  the 
heating  elements  to  the  atmosphere  vice  going  to  the  coolant  Heat  loss  experiments 
are  discussed  in  Appendix  B. 

The  second  experiment  was  to  account  for  two  additional  potential  heat 
sources:  heat  leakage  into  the  coolant  from  ambient  (q^J  and  frictional  temperature 
rise  (qtnc).  Approximately  8  cm  at  either  end  of  the  tube,  along  with  the  flexible 
tubing  section  prior  to  the  thermocouple  wells,  were  not  covered  by  the  heating 
elements.  Since  coolant  inlet  temperatures  were  primarily  below  ambient 
temperature,  a  potential  for  heat  leakage  to  the  coolant  from  ambient  existed.  Also, 
because  of  these  cold  inlet  temperatures,  frictional  heating  could  exist,  especially 
when  using  the  various  inserts.  Heat  leakage  into  the  system  and  frictional  effects 
were  considered  together  as  one  heat  influx.  These  experiments  were  conducted  with 
the  heater  tapes  turned  off;  any  temperature  rise  was  then  monitored  using  the  inlet 
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and  outlet  thermocouples.  Appendix  C  discusses  the  experiments  conducted  to 
determine  this  part  of  the  heat  balance  equation.  As  mentioned  earlier,  the  test 
section  and  all  other  tubing  were  insulated  to  minimize  these  effects. 


HEAT  BALANCE 
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INLET  / 
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Figure  4.1  Heat  Balance. 

B.  PROCEDURES  FOR  DETERMINATION  OF  INSIDE  HEAT  TRANSFER 
RATE 

1.  General  Procedures 

Each  tube  utilized  in  this  study  was  tested  with  no  insert  and  the 
appropriate  twisted  tape  and  HEATEX  inserts  corresponding  to  the  tube  inside 
diameter.  For  each  insert  condition,  four  nominal  coolant  inlet  temperatures,  -10°C, 
0°C,  +  10°C  and  +20°C  were  chosen  to  increase  the  range  of  experimental  Reynolds 
and  Prandlt  numbers  that  could  be  attained.  To  check  that  heat  flux  variations  did 
not  affect  the  results,  data  were  taken  at  heat  fluxes  of  10  kW/m:  and  15  kW/m2  for 
each  inlet  temperature.  For  each  insert,  coolant  inlet  temperature  and  heat  flux 
setting,  coolant  velocity  was  varied  between  0.2  m/s  and  1.4  m/s.  Since  the  inside 
diameter  of  the  Cu/Ni  finned  tube  was  smaller  than  the  other  two  tubes,  higher 
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coolant  velocities  of  up  to  2.4  m/s  could  be  attained.  A  matrix  of  the  individual  data 
runs  and  data  tile  names  is  shown  in  Table  4.1. 

With  the  tube  installed  and  all  calibration  experiments  completed,  the 
experiments  to  determine  the  inside  heat  transfer  coefficient  were  conducted.  Since 
flowmeter  settings  for  specific  coolant  velocities  were  based  upon  the  inlet  coolant 
temperature,  the  system  was  initially  operated  at  a  100%  flowmeter  setting  with  no 
power  to  the  heating  elements  for  15  minutes.  This  allowed  the  inlet  thermocouple 
well  to  stabilize.  Once  the  desired  inlet  coolant  temperature  was  reached  and 
maintained,  power  to  the  heating  elements  was  set  to  the  predetermined  heat  flux 
setting  for  that  data  run  (10  or  15  kW/m2).  The  flowmeter  was  then  adjusted  to 
achieve  the  desired  coolant  velocity  based  upon  flowmeter  calibration  results 
(discussed  in  Appendix  A).  Flow  was  maintained  for  approximately  10  minutes  to 
ensure  steady  state  equilibrium  conditions  prior  to  taking  data. 


TABLE  4.1  DATA  FILE  NAMES 


TEMP(C) / 

Cl 

Smooth  Tt 

be 

Cu/Nl  Finned 

Tub* 

Cu/N: 

Korodense  Tube  ! 

q  •<•»/■*> 

Hfii  LIjlJ 

K  J 

MBTBM 

HEATEX 

PS! 

■ 

■RTTH 

HEATEX 

IBM 

■oWItI 

nm 

HEATEX 

20/10 

SUN  1 04 

• 

3HHX04 

CFNI04 

C13) 

ESI 

■SRH 

9 

■Bill 

KQMI01 

X0TT01 

K0HX01 

20/13. 

SMNI04 

• 

SUHX04 

,3521 

■coys 

1 

mmnm 

warn 

KONIOI 

KDTT01 

KDHX01 

HMfl 

E5H31 

SNTT01 

hi 

mm 

wsm 

IIS-SSI 

10/15 

SUN  1 03 

SMTT01 

SMHX01 

. 

CFTT02 

K0NI02 

K0TT02 

K0HX02 

0/10 

3MNI02 

34TTT02 

CFNI03 

CFTT03 

CFHX03 

K0NIO3 

K0TT03 

KOHX03 

mm 

SMNI02 

SMHX02 

CFNI00 

CFTT03 

CFHX03 

KDNI03 

KDTT03 

K0HX03 

■10/10 

SUNI01 

SJ4TT03 

SUHX03 

CFNI02 

CFTT04 

CFHX04 

KDNI04 

K0TT04 

KDHX04 

-10/15 

SMNI01 

3MTT03 

3MHX03 

CFNI02 

CFTT04 

CFHX04 

KDNI04 

KDTT04 

^<0HX04 
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Velocities  were  increased  from  0.2  m/s  to  1.4  m/s  in  steps  of  0.2  m/s  and 
then  decreased  from  1.3  m/s  to  0.3  m/s  in  similar  steps.  This  method  of  increasing 
on  even  values  and  decreasing  on  odd  values  was  done  not  only  to  determine  if  any 
hysteresis  effects  existed  (i.e.,  a  difference  in  the  increasing  and  decreasing  data),  but 
also  to  ultimately  increase  the  number  of  data  points  acquired.  Also,  since  data  were 
immediately  plotted  (as  Nua  versus  Re),  the  method  allowed  for  an  easy 
determination  as  to  whether  a  data  point  had  reached  a  steady  state  condition. 

Inlet  and  outlet  coolant  temperatures,  tube  wall  temperatures,  heat  flux 
and  heating  element  power  settings  were  all  continuously  monitored  using  the 
computer  program  DRPSING  in  conjunction  with  the  Data  Acquisition  System 
discussed  earlier  (DRPSING  is  listed  in  Appendix  G).  Individual  data  runs  for  each 
tube  and  insert  condition  were  conducted  in  a  similar  manner. 

2.  Specific  Experiments 
a.  Copper  Smooth  Tube 

The  first  tube  tested  was  the  copper  smooth  tube.  This  tube  was 
chosen  to  provide  a  baseline  set  of  data  needed  for  validation  of  system  performance 
and  for  comparison  with  data  acquired  using  other  tubes.  It  should  also  be  noted  that 
two  complete  independent  sets  of  data  covering  all  insert  conditions  and  inlet 
temperatures  were  taken  over  a  period  of  two  weeks.  Data  from  both  sets  were 
compared  and  found  to  be  in  veiy  good  agreement  (±2%);  this  validated  the 
repeatability  of  data  obtained  from  the  system  under  the  same  conditions. 
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The  initial  intention  for  all  tubes  was  to  start  with  a  given  insert  at 
the  highest  coolant  inlet  temperature,  take  data  at  each  of  the  heat  flux  settings  and 
then  proceed  to  the  next  coolant  temperature.  Once  one  insert  was  completed,  the 
system  could  then  be  shut  down  and  the  insert  changed.  This  procedure  would  then 
be  repeated  for  all  three  insert  conditions.  However,  the  coolant  supply  tank  and 
coolant  system  were  used  concurrently  for  other  experiments  which  required  the 
coolant  temperature  to  be  -20°C.  Since  it  took  approximately  36-48  hours  foi  the 
coolant  temperature  to  increase  from  -10°C  to  +20°C  (the  sump  was  not  fitted  with 
a  heater),  only  the  copper  smooth  tube  was  done  in  this  way.  The  method  for  the 
other  two  tubes  is  discussed  below. 

b.  Cu/Ni  Externally  Finned  Tube 

The  second  tube  tested  was  the  Cu/Ni  finned  tube.  For  a  given  inlet 
temperature  and  tube  insert,  the  same  procedure  as  with  the  copper  smooth  tube  was 
used.  However,  rather  than  decreasing  coolant  temperature  with  the  same  insert,  the 
coolant  temperature  was  maintained  and  the  insert  changed  until  all  three  of  the 
insert  Conditions  were  completed  for  a  given  inlet  temperature.  Once  all  the  data  for 
a  given  inlet  temperature  had  been  taken,  the  supply  tank  temperature  was 
decreased  to  the  next  inlet  temperature  condition.  It  took  only  about  one  hour  for 
the  secondary  coolant  system  to  decrease  the  temperature  by  10°C.  Since  it  was  a 
relatively  short  task  to  change  inserts,  the  36  to  48  hour  waiting  period  required  to 
raise  coolant  temperatures  back  to  +20°C  between  insert  changes  was  eliminated. 
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It  should  be  noted  that  for  the  Cu/Ni  finned  tube,  the  highest  nominal 
inlet  temperature  attained  was  only  +  15°C  vice  +20°C.  Colder  than  normal  ambient 
temperatures  existed  preventing  supply  tank  temperatures  from  reaching  +20°C. 
Also,  due  to  the  smaller  inside  diameter,  coolant  velocities  could  be  varied  between 
0.3  m/s  and  2.4  m/s. 

c.  Cu/Ni  Korodense  Tub* 

The  same  procedure  used  with  the  finned  tube  was  also  used  with  the 
Korodense  tube.  In  addition,  another  test  was  conducted  with  this  tube.  It  was  noted 
that  larger  than  normal  wall  temperature  variations  existed  at  one  longitudinal 
location  when  using  the  two  inserts.  It  was  suspected  that  the  point  of  contact 
between  the  insert  and  tube  wall  wasLdirectly_below  a  thermocouple  position  causing 
a  higher  wall  temperature  at  that  point  To  verify  this,  the  system  was  shut  down  by 
securing  the  coolant  pump  and  the  heating  elements  so  that  the  flowmeter  and 
power  settings  remained  exactly  the  same.  The  outlet  end  of  the  tube  was  opened 
and  the  insert  moved  downstream  by  approximately  2  cm.  The  system  was  then 
closed  and  the  coolant  flow  and  power  to  the  heating  elements  restored.  Wall 
temperature  variations  still  existed  but  were  significantly  less.  This  procedure  was 
repeated  several  times  with  similar  results. 
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V.  DATA  REDUCTION 


A.  AVERAGE  INSIDE  HEAT  TRANSFER  COEFFICIENT 

The  main  objective  of  the  data  reduction  scheme  was  to  determine  the  average 
inside  heat  transfer  coefficient  (h;).  This  in  turn  was  used  to  determine  a  mean 
Nusselt  number  (for  the  whole  heated  length  of  the  tube)  to  be  used  as  the  basis  for 
establishing  the  desired  correlations. 

To  accomplish  this  task,  the  Data  Acquisition  Unit  discussed  in  Chapter  III 
made  use  of  the  computer  program  DRPSING.  A  listing  of  the  program  is  given  in 
Appendix  E.  The  program  allowed  the  user  t®  monitor  information  from  the 
instrumentation  attached  to  the  apparatus,  acquire  and  then  process  the  data  to 
determine  the  desired  output  In  doing  this,  the  user  would  select  the  tube  and  insert 
condition  as  well  as  the  nominal  coolant  inlet  temperature  for  that  particular  data 
run.  The  desired  coolant  velocity  was  then  selected.  Based  on  the  actual  measured 
coolant  inlet  temperature,  coolant  mixture  properties  and  flowmeter  calibration 
results  (see  Appendix  A),  the  appropriate  flowmeter  setting  (as  a  %)  was  calculated 
and  set  on  the  flowmeter. 

The  program  then  monitored  coolant  inlet  and  outlet  temperatures,  the  twelve 
wall  thermocouple  temperatures  and  the  actual  measured  heat  flux  and  power 
settings.  Once  all  readings  reached  a  steady  state  condition  (approximately  10  min.). 


45 


the  user  could  terminate  the  display  session  and  acquire  data  for  that  coolant  velocity 
setting. 

The  DAU  measures  emf  from  the  thermocouples  and  then  converts  these 
values  to  temperatures  using  the  temperature  calibration  equation  developed  by 
Mabrey  [Ref.  5].  Heater  power  was  determined  by  measuring  the  voltage  at  the  plug 
panel  and  calculating  the  power  based  on  heater  resistance  and  the  power  calibration 
contained  in  Appendix  D. 

The  twelve  wall  thermocouples  (placed  axially  and  radially)  were  used  to 
minimize  the  uncertainty  in  determining  the  average  inside  wall  temperature.  These 
twelve  wall  temperatures  were  averaged  and  corrected  (to  account  for  the  depth  of 
burial  using  a  simple  conduction  calculation)  to  determine  the  average  inside  wall 
temperature,  T^. 

The  physical  and  thermodynamic  properties  of  the  ethylene  glycol/water 
mixture  (density,  dynamic  viscosity,  thermal  conductivity,  heat  capacity)  were 
determined  using  equations  and  figures  from  Cragoe  [Ref.  22]  and  Gallant  [Ref. 
23]  and  correlated  by  Mabr**y  [Ref.  5].  From  these  quantities,  the  coolant  mass 
flow  rate  was  calculated  by: 

m  =  VcPcA,  (51) 

where  A*  is  the  inside  cross  sectional  area  of  the  tube  and  Vc  is  the  set  coolant 
velocity.  The  coolant  Reynolds  number,  ReD  was  calculated  by: 
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ReD  = 


K  Dfic 

Pc 


(5.2) 


where  D,  is  the  inside  diameter  of  the  tube.  The  dimensionless  distance  (x‘)  in  the 
flow  direction  is  calculated  as: 


jc*  - 


xJD. 

ReDPr 


(5.3) 


where  x  is  the  axial  distance  from  the  beginning  of  the  heated  length. 

To  determine  the  heat  transfer  rate  to  the  coolant,  several  factors  were 
considered.  The  basic  equation  to  determine  the  heat  transfer  rate  is  given  by: 


d  »  me  AT  (5-4) 

*  P 

* 

where  AT  is  the  measured  coolant  temperature  rise  along  the  tube.  The  temperature 
rise  required  is  that  for  the  actual  heated  length  of  the  tube  that  is  subjected  to  the 
constant  heat  flux  condition,  ATHL.  To  determine  ATHL  from  AT,  corrections  have  to 
be  made  to  account  for  heat  loss,  heat  leakage  and  heat  transferred  to  the  unheated 
ends.  Figure  5.1  depicts  a  typical  tube  and  the  correction  terms  considered  in  the 
heat  balance. 

The  overall  AT  is  first  corrected  for  heat  loss  to  the  ambient,  (q^)  and  heat 
leakage  into  the  system,  (q^).  The  heat  leakage  term  consisted  of  any  temperature 
rise  due  to  the  influx  of  heat  from  ambient  as  well  as  that  caused  by  frictional 
effects.  Appendix  D  discusses  the  experiments  conducted  to  determine  both  this  heat 
loss  and  heat  leakage  correction. 
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TEMPERATURE  RISE  CORRECTIONS 


Figure  5.1  Temperature  Rise  Corrections. 

To  take  account  of  the  heat  transferred  to  the  coolant  from  the  unheated  ends, 
the  heat  transfer  rate  to  the  coolant  over  the  entire  length  of  the  tube  (heated  and 
unheated)  is  first  calculated.  The  average  heat  transfer  coefficient  can  then  be 
determined  by: 


K. 

< 


Qhl 

AtLMTD 


(5.5) 


where  A,  is  the  inside  surface  area  based  on  the  heated  length  (LH)  and  LMTD  is 
the  log  mean  temperature  difference  defined  as: 

At  this  point  in  the  calculation,  however,  the  heat  transfer  rate  to  the  coolant 
from  the  heated  length  alone  (qHL)  is  not  known.  Since  the  heating  elements  did  not 
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LMTD 


(5.6) 
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cover  the  entire  length  of  the  tube,  account  must  be  taken  of  heat  added  to  the 
coolant  from  the  unheated  lengths  at  the  inlet  and  outlet  of  the  tube.  This  was 
accomplished  by  considering  these  ends  as  annular  fins.  An  adiabatic  fin  tip  and 
outer  wall  boundary  condition  was  applied  and  the  heat  transfer  from  the  fins,  %  (as 
given  by  Incopera  and  Dewitt  [Ref.  10]),  determined  from: 


«/ 


,n/pkA, 


QjanhmLf 


where  L,  *  fin  length 
P  =  fin  perimeter  =  irDi 
Ag  =  fin  cross  sectional  area 


0b  =  T^-T*  (for  inlet  fin) 
=  Twl-Tout  (for  outlet  fin) 


The  fin  calculation  requires  a  known  average  inside  heat  transfer  coefficient  (h;). 
Therefore,  an  initial  estimate  of  hj  is  calculated  by: 


A.LMTD 


(5.8) 


where  qtub.  is  known.  The  amount  of  heat  that  passes  through  the  unheated  ends  is 
then  calculated  using  the  fin  correction  above.  The  heat  transfer  rate  to  the  coolant 
from  the  heated  length  (4hl)  &  then  calculated  by: 


4m.  =  4  **-4/  (5-9) 

A  new  average  inside  heat  transfer  coefficient  (hj  can  then  be  calculated  using 
equation  5.5  and  an  iteration  process  carried  out  until  convergence  to  a  final  average 
inside  heat  transfer  coefficient  is  reached. 

An  additional  correction  for  the  Cu/Ni  finned  tube  is  needed  to  account  for  the 
change  in  inside  diameter  of  the  smooth  ends  (needed  to  fit  through  the  tube  sheet 
in  the  condenser).  The  heat  transfer  coefficient  based  on  the  smaller  inside  diameter 
corresponding  to  the  finned  portion  of  the  tube  is  what  ultimately  is  desired.  As 
shown  in  Figure  5.2,  the  heated  length  extended  beyond  the  smaller  diameter  section. 
To  correct  for  this,  the  annular  fin  correction  discussed  earlier  was  applied  to 
determine  the  temperatures  at  the  inlet  and  outlet  of  the  heated  length.  The  coolant 
temperature  rise  along  the  heated  length  was  assumed  to  be  linear,  and  therefore  the 
temperature  rise  along  the  smooth  larger  inside  diameter  lengths  (Lt  and  L2)  could 
be  determined  by  the  simple  ratio: 
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(5.10) 


AT,  -  LfISu  and  A Tt 


L£Tk 


HL 


-HL  HL 

where  Lhl  is  the  heated  length  of  the  tube  and  AThl,  the  temperature  rise  along  the 
heated  length.  These  corrections  were  then  added  or  subtracted  from  the  corrected 
finned  inlet  and  outlet  temperatures. 


Cu/Ni  EXTERNALLY  FINNED  TUBE 
TEMPERATURE  RISE  CORRECTIONS 


inlet  —  AT.—  •—  AT  —  outlet 
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Figure  5.2  Cu/Ni  Finned  Tube  Temperature  Rise  Corrections. 

With  the  average  inside  heat  transfer  coefficient  known,  the  mean  Nusselt 
number  can  be  determined  by: 


Num 


(5.11) 


where  k  is  the  thermal  conductivity  of  the  ethylene  glycol/water  mixture. 
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The  experimentally  determined  mean  Nusselt  number  was  then  compared  to 
various  theoretical  and  experimental  correlations,  based  on  the  insert  condition. 


These  correlations  were  discussed  earlier  but  are  briefly  mentioned  here. 

For  the  no  insert  condition,  the  data  were  compared  to  two  theoretical 
correlations,  the  Hausen  correlation  [Ref.  10]  given  by: 


Su  s  3.66+ 

m  U0.V.[(D/L)ReDPr]z'1 


(5.12) 


and  the  Shah  and  London  correlation  [Ref.  9]  given  by: 

Num  =  1.93  (x*)'03  (513> 

The  reader  is  reminded  that  the  Hausen  correlation  is  for  a  uniform  wall 
temperature  boundary  condition  and  was  used  in  Mazzone’s  condensation 
experiments  [Ref.  6].  The  Shah  and  London  correlation  is  for  a  constant  wall  heat 
flux  condition.  It  should  also  be  noted  that  equation  5.13  is  used  since  x*  for  the  test 
section  was  always  less  than  0.03. 

For  the  twisted  tape  insert  the  experimental  mean  Nusselt  number  was 
compared  to  the  correlation  developed  by  Hong  and  Bergles  [Ref.  15]: 

Num  =  0.383Pr0JS(Re  Jy)0'632  (5-14) 

Equation  5.14  is  used  because  values  for  Re^y  were  always  greater  than  10. 


For  the  HEATEX  insert,  the  data  were  compared  to  the  two  correlations 
developed  by  Mazzone  [Ref.  6].  For  the  13.26  mm  inside  diameter  tube,  the  equation 
is  given  by: 

Num  =  Q.226Re0MPr0'46  (5.16) 

and  for  the  10.16  mm  inside  diameter  tube  by: 

Num  =  0.063 Re016Pr0M  (5.16) 

B.  ADDITIONAL  TUBE  LENGTHS 

In  addition  to  the  above  calculations  for  the  entire  tube  length,  mean  Nusselt 
numbers  for  three  additional  tube  "lengths"  could  be  determined.  This  was 
accomplished  by  sectioning  the  heated  length  of  the  tube  into  four  cells  as  shown  in 
Figure  5.3.  Each  cell  contained  one  longitudinal  thermocouple  position.  Three 
additional  tube  lengths  were  then  determined  by  adding  cell  lengths  together.  The 
first  length,  Llt  consisted  of  cell  1  alone,  L;  consisted  of  cells  1  and  2  and  Lj 
consisted  of  cells  1,  2  and  3.  The  fourth  length  was  actually  the  entire  heated  length 
as  previously  discussed. 

The  coolant  temperature  rise  along  the  heated  length  of  the  tube  was  assumed 
to  be  linear.  The  correction  to  the  wall  heat  flux  to  account  for  variation  in  electrical 
resistance  with  temperatures  was  assumed  to  be  negligible.  Therefore,  the 
temperature  rise  along  each  cell  could  be  determined  by: 
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SECTIONED  TUBE 


*”HL 

T/C  POS  1  T/C  POS  2  T/C  POS  3  T/C  POS  * 

l_ _ 

INLET 

CELL  1 

«lf  2 

' 

CELL  3 

I 

CELL,  4 

OUTLET 

r.  : - i 

Figure  5.3  Additional  Tube  Lengths. 

ar_ .  Epi.  (5.17) 

4 

Average  inside  wall  temperatures  were  based  on  the  thermocouples  within  that  tube 
"length."  Therefore,  the  temperature  rise,  heat  flux  and  average  inside  wall 
temperatures  were  all  known  for  each  tube  length.  Applying  the  same  calculation 
method  as  before,  the  mean  Nusselt  number  for  these  three  additional  tube  "lengths" 
could  be  calculated. 

This  adaption  to  the  program  was  not  done  until  after  the  smooth  tube 
experiments  had  been  completed;  however,  this  extra  data  were  taken  for  the 
remaining  two  tubes. 


54 


C.  LOCAL  NUSSELT  NUMBERS 


The  program  was  also  adapted  to  determine  local  inside  heat  transfer 
coefficients  at  the  four  axial  thermocouple  positions.  Corrections  discussed  earlier 
to  determine  the  temperatures  at  the  inlet  and  outlet  of  the  heated  length  were  first 
accounted  for.  Exact  thermocouple  longitudinal  positions  were  entered  into  the 
program  as  L/D  ratios  so  that  a  local  x"  could  also  be  calculated.  Again,  the 
temperature  rise  along  the  heated  length  was  assumed  to  be  linear;  therefore,  the 
bulk  coolant  temperature  at  the  particular  thermocouple  position  could  be 
determined  by  the  simple  ratio  method  discussed  earlier  (equation  5.10). 

The  local  heat  transfer  coefficient,  h^,  was  then  determined  by: 


T  -T 

m  c 


(5.13) 


where  T*.  is  the  average  inside  wall  temperature  based  on  the  three  thermocouples 
at  that  position  and  Tc  is  the  bulk  coolant  temperature  at  that  position.  The  local 
Nusselt  number  was  then  given  by: 


M<,  .  (5.19) 

The  local  Nusselt  number  for  all  four  thermocouple  positions  was  then  determined. 

For  the  smooth  tube  experiments,  local  Nusselt  numbers  were  calculated  only 
for  position  1.  The  program  was  later  adapted  to  include  local  calculations  at 
positions  2,  3,  and  4  for  the  Cu/Ni  finned  tube  and  Korodense  tube. 
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VI.  RESULTS  AND  DISCUSSION 


A.  GENERAL 

The  experimental  heat  transfer  results  are  initially  presented  in  terms  of  the 
mean  Nusselt  and  Reynolds  numbers  for  all  insert  conditions.  From  these  graphs,  it 
can  be  seen  how  the  Nusselt  number  is  inlet  temperature  dependent,  indicating  the 
need  for  a  Prandtl  number  (raised  to  an  exponent)  in  the  final  correlation.  However, 
although  not  specifically  shown,  the  data  shows  no  dependence  on  the  heat  flux 
setting.  In  addition,  the  Nu  vs.  Re  graphs  for  the  no  insert  condition  show  when  a 
*  transition  to  turbulent  flow  occurs. 

The  heat  transfer  results  are  then  compared  to  existing  correlations  based  upon 
the  insert  condition.  For  the  no  insert  condition,  the  Nusselt  number  is  plotted 
against  x*  (dimensionless  x  coordinate).  Note  that  for  a  given  tube,  L/D  is  constant 
and  x‘  varies  due  to  Re  and  Pr  changes.  The  data  are  then  compared  to  the  Hausen 
[Ref.  10]  correlation  (equation  2.10)  and  the  Shah  and  London  [Ref.  9]  correlation 
(equation  2.14).  Recall  that  the  Hausen  correlation  is  for  a  uniform  wall  temperature 
boundary  condition  and  was  used  by  Mazzone  [Ref.  6]  in  the  refrigerant 
condensation  experiments.  Since  the  present  data  were  taken  under  a  constant  heat 
flux  boundary  condition,  comparison  with  the  Shah  and  London  correlation  was  also 
required.  The  twisted  tape  insert  data  are  compared  to  the  correlation  developed  by 


56 


Hong  and  Bergles  [Ref.  15]  (equation  2.20)  which  was  also  used  by  Mazzone.  The 
HEATEX  insert  data  are  compared  to  the  two  correlations  developed  from  the 
limited  data  received  from  CAL-GAVIN  by  Mazzone  (equations  2.27  and  2.28). 

Correlations  are  developed  for  each  tube  and  insert  condition  based  upon  a 
least  squares  linear  regression.  An  r-squared  value,  a  measure  of  the  accuracy  of  the 
least  squares  curve  fitting  to  the  data,  is  used  to  determine  the  uncertainty  of  the 
correlation.  The  correlations  and  the  corresponding  r-squared  values  are  summarized 
in  Table  6.1.  Note  that  the  0.2  m/s  and  0.3  m/s  data  were  eliminated  prior  to 
developing  all  correlations  due  to  the  high  degree  of  uncertainty  in  these  data  points. 
TABLE  6.1  SUMMARY  OF  CORRELATIONS 

SUMMARY  OF  CORRELATIONS 


NO  INSERT 

TWISTED 

TAPE 

HEATEX 

SMOOTH 

Nu  ■ 

Nu  • 

Nu  - 

4.13  7Xstar*-0. 234 

.5l2Pr.35(R0S/yr.5e8 

.104Pr*46Re~04l 

TUBE 

r-sq  -  71.9% 

r-sq  -  99.5% 

r-sq  -  98.9% 

EXT 

Nu  ■ 

Nu  - 

Nu  * 

FINNED 

1.84Xstar*-0.341 

148Pr‘35(R0S/y)'.687 

.126Pr'46Re'.857 

TUBE 

r-sq  -  93.0% 

r-sq  •  96.7% 

r-sq  *  95  9 

KORODENSE 

Nu  • 

Nu  • 

Nu  ■ 

2O47Xstar*-O320 

509Pr'.35(Res/y)*498 

282Pr“40Re*  499 

TUBE 

r-sq  •  89.2% 

r-sq  ■  98.9% 

r-sq  *  94.3% 
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B.  SMOOTH  TUBE 


1.  Smooth  Tube  -  No  Insert 

The  first  tube  tested  was  the  copper  smooth  tube  with  no  insert  Heat 
transfer  data  is  presented  in  Figure  6.1  in  terms  of  the  Nusselt  number  as  a  function 
of  Reynolds  number  for  each  nominal  inlet  temperature.  The  first  thing  to  notice  is 
the  sudden  increase  in  Nu  at  a  Re  number  of  about  2300.  This  is  indicative  of  a 
change  from  laminar  to  transition  flow  where  the  heat  transfer  is  greatly  increased 
due  to  unproved  mixing.  This  contradicts  Mazzone’s  [Ref.  6]  assumption  of  a 
transition  to  turbulence  occurring  at  Re  >  4000.  The  data  in  the  transition  region 
was  primarily  for  the  +20°C  and  +10°C  nominal  inlet  temperatures.  Figure  6.1  also 
shows  distinct  groupings  of  data  at  given  inlet  temperatures,  indicating  a  Prandtl 
number  effect  These  groupings  are  shown  more  clearly  in  Figure  6.2  which  represent 
data  under  laminar  flow  conditions  only  (Re  £  2300)  and  Figure  6.3  for  data  in  the 
transition  region  only  (Re  >  2300). 

To  develop  a  useful  heat  transfer  correlation,  the  mean  Nusselt  number 
is  plotted  as  a  function  of  x*  in  Figure  6.4.  Recall  that  x*  can  be  determined  by: 


x'  = 


UP 

ReDPr 


(2.9) 


For  a  given  tube,  L/D  is  a  fixed  quantity.  Therefore,  for  a  given  Prandtl  number 
(essentially  fixed  inlet  temperature),  x*  is  dependent  only  upon  the  Reynolds  number. 
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Figure  6.4  also  shows  a  comparison  of  the  experimental  data  with  the 
Hausen  correlation  (equation  2.10)  for  a  uniform  wall  temperature  (used  by 
Mazzone)  and  the  Shah  and  London  correlation  (equation  2.14)  for  a  constant  heat 
flux.  The  Hausen  correlation  appears  to  significantly  underestimate  the  heat  transfer 
coefficient  found  in  the  present  experimental  data,  i.e.,  the  present  experimental  data 
indicate  a  lower  inside  thermal  resistance.  This  is  consistent  with  Mazzone’s  [Ref.  6] 
results  in  that  the  Hausen  correlation  was  predicting  larger  inside  resistances  than 
the  overall  measured  resistance.  However,  it  must  always  be  remembered  that  the 
Hausen  correlation  is  for  a  uniform  wall  temperature  boundary  condition  and  this 
would  indeed  lead  to  lower  heat* transfer  coefficients. 

The  data  also  fall  slightly  above  (but  pairallel  to)  that  predicted  by  the 
Shah  and  London  correlation.  Due  to  the  constant  heat  flux  boundary  condition  used 
by  Shah  and  London,  one  may  expect  better  agreement  (as  seen).  However,  the  fact 
that  the  data  fall  slightly  above  can  be  attributed  to  secondary  flow  induced  by  the 
apparatus  inlet  conditions  (90°  bend  due  to  the  thermocouple  well).  Also,  as  before, 
a  departure  from  laminar  flow  can  be  seen  in  Figure  6.4.  Since  the  transition  data 
forms  two  distinct  lines  (+20°C  and  +1Q°C),  it  indicates  that  x*  is  not  a  good 
correlating-  parameter  for  transitional  flow. 

After  removal  of  the  transitional  flow  data,  a  correlation  was  developed 
based  upon  the  form  of  the  Shah  and  London  correlation.  The  resulting  present 
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correlation  for  the  mean  Nusselt  number  for  laminar  flow  with  no  insert  is  given  by: 

Num  *  4.Yyi(x')~0:a*  (6.1) 

Equation  6.1  and  the  Hausen  correlation,  equation  2.10,  are  both  shown 
in  Figure  6.5  with  the  present  laminar  flow  data.  The  general  form  of  Figure  6.5  is 
explained  by  the  fact  that  at  low  x*  (high  Re)  the  heat  transfer  is  improved  due  to 
slightly  better  mixing.  As  x*  increases  (coolant  velocity  decreases),  the  heat  transfer 
decreases. 

Figure  6.5  also  shows  a  high  degree  of  scatter  in  the  data.  Initially,  it  was 
thought  that  this  scatter  was  due  to  the  flow  not  being  hydrodynamically  fully 
developed.  However,  it  is  seen  to  be  more  pronounced  at  large  x*.  Recall  that  for  a 
given  nominal  inlet  temperature,  large  x*  represents  a  low  Reynolds  number.  For  a 
low  Reynolds  number,  the  hydrodynamic  entry  length  (given  by  X/D  =  0.05  Re)  is 
the  shortest  Also,  since  we  are  dealing  with  relatively  large  Prandtl  numbers  (SO¬ 
HO),  the  assumption  that  the  velocity  profile  is  fully  developed  throughout  the 
thermal  entry  region  should  be  fairly  valid.  Therefore,  it  is  felt  that  this  scatter  is  not 
due  to  the  flow  being  hydrodynamically  undeveloped. 

This  scatter  is  believed  to  be  due  to  the  effects  of  the  different  inlet 
temperatures,  indicating  that  x'  is  not  an  ideal  correlating  parameter.  That  is  to  say, 
the  Prandtl  number  exponent  needs  to  be  adjusted  to  account  for  the  different  inlet 
temperatures.  It  should  also  be  noted  that  due  to  the  high  degree  of  detail  in  Figure 
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6.5,  the  scatter  appears  more  pronounced.  However,  in  Figure  6.4  with  larger  scale 
increments,  the  scatter  seems  minimal. 

2.  Smooth  Tube  -  Twisted  Tape 

Results  for  the  smooth  tube  with  a  twisted  tape  insert  are  shown  in 
Figures  6.6  and  6.7.  A  definite  Prandtl  number  dependence  can  be  seen  in  Figure  6.6 
where  the  mean  Nusselt  number  is  presented  as  a  function  of  Reynolds  number. 
Following  the  guidance  of  Hong  and  Bergles  [Ref.  15]  the  dimensionless  parameters 
Re,  and  y  were  combined  as  one  parameter  Re^y  and  the  mean  Nusselt  number  and 
Prandtl  number  were  combined  as  Nu/Pr0Js.  The  final  correlation  shown  in  Figure 
6.7  is: 

Num  =  0.512  Pr'XiReJy)0-**  (6-2) 

This  correlation  is  in  close  agreement  with  the  Hong  and  Bergles  correlation 
(equation  2.20).  The  slight  differences  can  be  attributed  to  differences  in  the 
experimental  apparatus  (including  the  insert  itself)  and  the  possible  existence  of 
secondary  flows.  Figure  6.7  shows  little  inlet  temperature  dependence,  indicating  that 
the  Prandtl  number  exponent  of  0.35  is  probably  fairly  accurate.  Furthermore,  there 
is  little  experimental  data  scatter,  giving  increased  confidence  in  this  correlation. 

3.  Smooth  Tube  •  HEATEX  Insert 

HEATEX  insert  results  for  a  smooth  tube  are  shown  in  Figures  6.S  and 
6.9.  Again,  on  a  Nu  vs.  Re  basis,  it  can  be  seen  from  Figure  6.8  that  there  is  a 
definite  dependence  on  the  inlet  temperature.  Using  the  same  form  as  the 
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correlation  of  Oliver  and  Aldington  [Ref.  20],  the  resulting  correlation  for  the  mean 
Nusselt  number  is: 

Num  *  0 .lMPr^Re06*1  (6.3) 

This  is  compared  in  Figure  6.9  to  the  correlation  developed  and  used  by  Mazzone 
(equation  2.23).  It  is  clearly  evident  that  Mazzone’s  correlation  underpredicts  the 
inside  thermal  resistance  and  therefore  gives  a  larger  inside  heat  transfer  coefficient. 
Again,  Figure  6.9  shows  little  dependence  on  inlet  temperature  indicating  that  the 
Prandtl  number  exponent  is  fairly  accurate.  Comparisons  of  the  enhancement  in  heat 
transfer  gained  by  the  use  of  HEATEX  and  twisted  tape  insert,  over  the  no  insert 
condition  are  given  section  F  of  this  chapter. 
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SMOOTH  TUBE  -  NO  INSERT 


Nu  vs  Re 


Reynolds 

.  +20C  ♦  +10C  A  OC  .  -IOC 

SMN1RENU 


Figure  6.1  Smooth  tube  -  No  insert  (NuB  vs.  Re)  results. 
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SMOOTH  TUBE  -  NO  INSERT 


LAMINAR  FLOW 


Reynolds 

•  +20C  ♦  +  10C  A  OC  .  -IOC 

Nu  vs  Re 
SMN  INURE  LAM 


Figure  6,2  Smooth  tube  -  No  insert  (Nu.  vs.  Re)  laminar  flow  results. 
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Figure  6.3  Smooth  tube  -  No  insert  (NuB  vs.  Re)  transitional  flow  results. 


SMOOTH  TUBE  -  NO  INSERT 


LAMINAR  FLOW 
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Figure  6.5  Smooth  tube  -  No  insert  laminar  flow  correlation  results. 


SMOOTH  TUBE  -  TWISTED  TAPE 


Nu  vs  Re 


„  +10C  ♦  OC  4  -IOC 
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Figure  6.6  Smooth  tube  •  Twisted  tape  (Num  vs.  Re)  results. 
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SMOOTH  TUBE  -  TWISTED  TAPE 


(Nii/Pr^SS)  vs  (Res/y) 


Res/y 


«  +10C  ♦  OC  A  -IOC 

Res  *  (mDi)/(Asjj) ,  y  =*  H/Di 
SMTTNUPR 


Figure  6.7  Smooth  tube  -  Twisted  tape  correlation  results. 
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Nu/PrA.46 


SMOOTH  TUBE  -  HEATEX 


(Nu/Pr\46)  vs  Re 


Reynolds 


.  +2oc  ;  +ioc  *  oc  «  -loc 

SMHXNUPR 

Figure  6.9  Smooth  tube  -  HEATEX  correlation  results. 
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C.  CU/NI  EXTERNALLY  FINNED  TUBE 


1.  Cu/Ni  Finned  Tube  •  No  Insert 

The  second  tube  tested  was  the  Cu/Ni  externally  finned  tube.  Recall  that 
this  tube  has  a  smaller  inside  smooth  diameter  than  the  smooth  tube  (10.15  mm  vs. 
13.26  mm).  The  no  insert  results  for  the  mean  Nusselt  number  as  a  function  of 
Reynolds  number  are  shown  in  Figure  6.10  and  in  more  detail  in  Figures  6.11  and 
6.12  for  laminar  and  transitional  flows,  respectively.  It  can  be  seen  that  a  slightly 
higher  Reynolds  number  (Re  2700)  for  the  departure  from  laminar  flow  exists  for 
this  tube.  Again,  the  dependence  on  inlet  temperature  is  clearly  evident 

The  data  are  compared  to  the  Hausen  correlation  (equation  2.10)  and  the 
Shah  and  London  correlation  (equation  2.14)  in  Figure  6.13.  As  expected,  the 
Hausen  correlation  again  underpredicts  the  inside  heat  transfer  coefficient  However, 
the  laminar  flow  data  is  in  closer  agreement  with  the  Shah  and  London  correlation 
than  they  were  for  the  smooth  tube.  This  can  be  attributed  to  the  different  diameter 
affecting  the  thermal  boundary  layer  development  The  resulting  correlation  for  the 
laminar  flow  mean  Nusselt  number  is  shown  in  Figure  6. 14  and  is  given  by: 

Num  =  1.84  (jc*)-0-341  (6.4) 

The  same  conclusions  about  data  scatter  can  be  drawn  as  for  the  smooth  tube 
results,  although  in  general  the  scatter  seems  to  be  less. 
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2.  Cu/Ni  Finned  Tube  -  Twisted  Tape  Insert 

Figures  6.15  and  6.16  show  the  results  for  the  Cu/Ni  finned  tube  using  a 
twisted  tape  to  enhance  tubeside  heat  transfer.  The  inlet  temperature  dependence 
in  Figure  6.15  seems  to  be  reduced,  but  still  evident  The  form  of  the  correlation  was 
again  that  of  Hong  and  Bergles’  correlation.  It  was  developed  in  the  same  manner 
as  for  the  smooth  tube  and  is  given  by: 

Num  *  0. 148/V 025  (Rejy)0-6*7  (6.5) 

The  data  and  resulting  correlation  are  compared  to  the  Hong  and  Bergles  correlation 
(equation  2.20)  in  Figure  6.16  where  the  latter  is  seen  to  overpredict  the  mean 
Nusselt  number.  This  difference  can  be  attributed  to  several  factors.  A  small 
diameter  tube  with  a  larger  diameter  smooth  entrance  length  could  cause  severe 
secondary  flows.  Also,  because  of  the  smaller  diameter  tube  a  different  size  twisted 
tape  (twist  ratio,  y=3)  was  used,  which  had  a  slightly  looser  fit  than  the  tape  used 
with  the  smooth  tube.  Additionally,  Prandtl  number  effects  need  to  be  investigated 
more  closely,  as  Figure  6.16  shows  a  definite  trend  with  inlet  temperature. 

3.  Cu/Ni  Finned  Tube  -  HEATEX  Insert 

Data  for  the  Cu/Ni  finned  tube  using  the  HEATEX  insert  are  presented 
in  Figures  6.17  and  6.18.  Figure  6.17  (Nu  vs.  Re)  again  shows  a  dependence  of  the 
data  on  the  coolant  inlet  temperature.  Since  the  HEATEX  element  used  for  this 
tube  was  a  smaller  diameter  (10.16  mm)  with  a  different  loop  density  than  that  used 
for  the  smooth  tube,  a  new  correlation  needed  to  be  developed.  Again,  using  the 
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form  of  the  correlation  proposed  by  Oliver  and  Aldington,  the  final  proposed  correlation 
is  given  by: 

Num  -  QA26Pr°MRe9a7  (6-6) 

The  data  and  resulting  correlation  are  compared  to  the  correlation  developed 
by  Mazzone  (equation  2.24)  for  this  insert  in  Figure  6. 18.  It  can  be  seen  that  Mazzone’s 
correlation  is  adequate  for  predicting  the  inside  heat  transfer  coefficients  at  low  Reynolds 
number  (Re  <  1400),  where  both  correlations  merge.  However,  at  higher  Reynolds 
numbers,  Mazzone’s  correlation  overpredicts  the  present  experimental  data.  The 
differences  in  the  two  correlations  are  primarily  due  to  the  limited  amount  of  data 
available  to  Mazzone  when  developing  his  correlation.  Furthermore,  the  data  he  used 
were  taken  on  a  different  experimental  facility  (by  the  manufacturers).  As  seen  before, 
the  data  shown  in  Figure  6. 18  still  indicate  a  dependence  on  inlet  temperature  and  hence 
the  value  of  0.46  for  the  Prandtl  number  exponent  needs  to  be  reinvestigated  and 
optimized. 
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Figure  6.10  Cu/Ni  externally  finned  tube  -  No  insert  (Nu.  vs.  Re)  results. 
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Mean  Nusselt 


Cu/Ni  FINNED  TUBE  -  NO  INSERT 

LAMINAR  FLOW 
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Nu  vs  Re 
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Figure  6.11  Cu/Ni  externally  finned  tube  -  No  insert  (NuB  vs.  Re)  laminar  flow 
results. 


76 


Mean  Nusselt 


Cu/Ni  FINNED  TUBE  -  NO  INSERT 


TRANSITIONAL  FLOW 
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Figure  6.12  Cu/Ni  externally  finned  tube  -  No  insert  (Nu„  vs.  Re)  transitional  flow 
results. 
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Cu/Ni  FINNED  TUBE  -  NO  INSERT 


Nu  vs  Xstar 


Xstar 

a  -IOC  ♦  OC  A  +10C  .  +15C 
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Figure  6.13  Cu/Ni  externally  finned  tube  -  No  insert  (NuB  vs.  x’)  results. 
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Figure  6. 14  Cu/Ni  externally  finned  tube  -  No  insert  laminar  flow  correlation  results. 
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Figure  6.15  Cu/Ni  externally  finned  tube  -  Twisted  tape  (Num  vs.- Re)  results. 
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Figure  6.16  Cu/Ni  externally  finned  tube  -  Twisted  tape  correlation  results. 
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Figure  6.18  Cu/Ni  externally  finned  tube  -  HEATEX  correlation  results. 
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D.  KORODENSE  TUBE 


1.  Korodense  Tube  •  No  Insert 

The  final  tube  tested  was  the  Korodense  tube.  This  tube  was  first  tested 
with  no  insert  The  results  for  the  mean  Nusselt  number  as  a  function  of  Reynolds 
number  are  shown  in  Figure  6.19  and  in  more  detail  for  laminar  and  transitional  flow 
in  Figures  6.20  and  6.21,  respectively.  As  was  seen  in  the  smooth  tube  case,  a 
transition  to  turbulent  flow  occurred  at  a  Reynolds  number  of  approximately  2300. 

Figure  6.22  compares  the  data  to  Hausen’s  correlation  (equation  2.10)  and 
the  Shah  and  London  correlation  (equation  2.14).  As  before,  the  transitional  flow 
data  are  seen  as  two  distinct  lines  fqr  data  taken  at  +20°C  and  +10°C  nominal  inlet 
temperatures.  Under  laminar  flow  conditions,  the  data  are  seen  to  fall  very  closely 
to  Shah  and  Londons’  correlation,  but  significantly  above  Hausen’s  correlation.  This 
was  also  found  for  the  other  two  tubes.  This  result  was  expected  since  a  mildly 
corrugated  wall  is  not  good  in  enhancing  laminar  flow  and  the  tube  behaves 
essentially  like  a  smooth  tube.  Significant  enhancement  for  Korodense  tubes  is  not 
expected  until  ReD  >  104  [Ref.  21]. 

The  resulting  correlation  for  laminar  flow  data  is  shown  on  Figure  6.23 
and  given  by: 


Num  *  2.05(*-)-°J“  (6.7) 

The  Hausen  correlation  (equation  2.10)  used  by  Mazzone  again  under  predicts  the 
present  measured  inside  heat  transfer  coefficient.  It  was  expected  that  the  mean 
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Nusselt  number  would  be  the  same  (or  maybe  slightly  higher)  as  for  the  smooth  tube 
case  for  the  reasons  already  mentioned.  However,  the  Korodense  data  fall  slightly 
below  the  smooth  tube  data  and  closer  to  Shah  and  Londons’  smooth  tube 
correlation  (equation  2.14).  No  explanation  can  be  given  for  this  result  Data  scatter 
explanations  given  for  the  other  tubes  also  apply  for  the  Korodense  tube. 

2.  Korodense  Tube  -  Twisted  Tape  Insert 

The  results  for  the  Korodense  tube  with  a  twisted  tape  insert  are  shown 
in  Figures  6.24  and  6.25.  The  same  twisted  tape  that  was  used  for  the  smooth  tube 
was  used  for  the  Korodense  tube.  Note  that  the  Korodense  tube  nominal  inside 
diameter  (13.35  mm)  is  slightly  larger  than  the  smooth  tube  inside  diameter  (13.26 
mm).  This  small  difference  in  inside  diameters  (together  with  corrugation  in  the 
Korodense  tube)  led  to  a  looser  fit  for  the  twisted  tape  insert 

The  mean  Nusselt  number  as  a  function  of  Reynolds  number  is  shown  in 
Figure  6.24  which  again  clearly  indicates  the  inlet  temperature  dependence.  The  data 
were  again  correlated  using  the  form  of  Hong  and  Bergles’  correlation.  The  present 
data  were  compared  to  their  correlation  (equation  2.20)  in  Figure  6.25.  The 
correlation  developed  for  the  present  data  is: 

Num  =  0.569Pr035(ReJy)0m  (6:8) 

This  correlations  falls  below  the  Hong  and  Bergles  correlation  (equation  6.20)  as  well 
as  .for  the  smooth  tube.  This  result  is  thought  to  be  primarily  due  to  the  looser  fit  of 
the  twisted  tape  inside  the  Korodense  tube.  The  corrugations  also  leave  areas  of  no 
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contact  between  the  tape  and  the  inside  wall  of  the  tube.  Stagnant  flow  or  even  air 
bubbles  may  form  in  these  gaps,  especially  at  low  Reynolds  numbers,  causing  an 
interruption  of  the  swirled  flow.  The  thermal  boundary  layer,  therefore,  is  not  as 
efficiently  mixed  as  in  the  smooth  tube.  For  these  reasons,  lower  inside  heat  transfer 
coefficients  than  those  assumed  by  Mazzone  should  be  expected. 

3.  Korodense  Tube  -  HEATEX  Insert 

The  last  tube-insert  condition  tested  was  the  Korodense  tube  with  the 
HEATEX  insert  The  same  HEATEX  insert  used  for  the  smooth  tube  was  used 
here.  Figure  6.26  shows  the  mean  Nusselt  number  as  a  function  of  Reynolds  number 
and  again  shows  an  inlet  temperature  dependence.  The  resulting  correlation  for  the 
present  data  is  given  by: 

Num  =  0.28 2Pr0MRe0A”  (69) 

The  data  and  this  correlation  are  compared  with  the  correlation  Mazzone  used 
(equation  2.23)  in  Figure  6.27.  It  is  evident  that  Mazzone’s  chosen  correlation 
overpredicts  the  inside  heat  transfer  coefficient  significantly.  This  again  is  thought 
to  be  primarily  due  to  the  incomplete  mixing  of  the  thermal  boundary  layer  caused 
by  non-contact  regions  existing  between  the  tube  wall  and  HEATEX  petals. 
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Figure  6.20  Korodense  tube  -  No  insert  (Nu,,,  vs.  Re)  laminar  flow  results. 
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Figure  6.21  Korodense  tube  -  No  insert  (NuB  vs.  Re)  transitional  flow  results. 
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Figure  6.23  Korodense  tube  -  No  insert  laminar  flow  correlation  results. 


91 


0 


500 


1000 


Res/y 

a  +20C  *  +10C  A  OC  .  -IOC 

Res  •  (mDi)/(AsM) ,  y  »  H/Di 
KDTTNUPR 


Figure  6.25  Korodense  tube  -  Twisted  tape  correlation  results. 
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Figure  6.27  Korodense  tube  -  HEATEX  correlation  results. 
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E.  INSERT  COMPARISON  BETWEEN  TUBES 

Figures  6.28,  6.29  and  6.30  each  compare  the  three  correlations  for  each  tube 
for  no  insert,  twisted  tape,  and  HEATEX  inserts,  respectively.  It  can  be  seen  that 
with  no  insert  (Figure  6.28),  the  Korodense  correlation  is  about  10%  below  the 
smooth  tube  data.  Although  one  may  expect  no  increase  in  the  heat  transfer  with  a 
roped  tube  under  laminar  flow  conditions,  one  certainly  would  not  expect  to  find 
such  a  decrease  with  a  similar  diameter  tube.  This  difference,  therefore,  can  only  be 
attributed  to  experimental  uncertainty,  probably  in  the  wall  temperature 
measurement  The  finned  tube  correlation  is  lower  than  both  the  Korodense  (=5%) 
and  smooth  tube  (=15%)  correlations;  this  tube  had  a  smaller  inside  diameter  and 
it  may  be  that  different  induced  flow  patterns  are  set  up  causing  a  reduction  in  heat 
transfer.  Experimental  uncertainty  must  also  not  be  ruled  out  as  to  the  cause  of  the 
difference. 

Figure  6.29  shows  the  same  three  tube  correlations,  but  for  the  twisted  tape 
insert  The  smooth  tube  again  has  the  highest  performance.  The  finned  tube  is  =20% 
below  at  a  high  Re  number  and  nearly  40%  lower  at  a  low  Re  number.  This  could 
be  due  to  a  different  twist  ratio  used  for  the  finned  tube  due  to  the  smaller  inside 
diameter.  The  Korodense  tube  is  similar  to  the  finned  tube,  giving  slightly  lower  heat 
transfer  at  high  Re  number  and  slightly  higher  heat  transfer  at  low  Re.  The  reason 
for  this  tube  being  significantly  lower  than  the  smooth  tube  is  though  to  be  due  to 
the  tape  not  contacting  the  inner  surface  of  the  tube  as  well  as  for  the  smooth  tube. 
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thereby  leaving  regions  of  stagnant  flow.  These  regions  would  tend  to  inhibit  mixing 
and  lead  to  a  lower  heat  transfer  performance. 

For  the  HEATEX  comparison  (Figure  6.30),  similar  conclusions  can  be  made, 
although  the  finned  tube  seems  to  have  better  relative  performance  indicating  that 
the  small  diameter  HEATEX  insert  used  in  the  small  diameter  tube  is  more  effective 
than  the  twisted  tape.  The  difference  between  the  smooth  and  Korodense  tubes  are 
of  similar  magnitude  to  that  shown  in  Figure  6.29  (about  40%  at  high  Re). 
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Figure  6.28  No  insert  comparison  under  laminar  flow  for  all  three  tubes. 
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Figure  6.29  Twisted  tape  insert  comparison  under  laminar  flow  for  all  three  tubes. 
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re  6.30  HEATEX  insert  comparisons  under  laminar  flow  for  all  three  tubes. 
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F.  LAMINAR  FLOW  ENHANCEMENT  RATIO 


The  comparisons  made  in  Figures  6.28  through  6.30  are  made  more 
complicated  by  the  fact  that  the  uncertainty  associated  with  the  wall  temperature 
measurement  is  unknown  and  most  likely  different  for  each  tube.  For  this  reason,  it 
is  not  known  how  much  of  the  divergence  seen  in  these  figures  is  due  to  differences 
in  thermocouple  attachment  (giving  uncertainty  to  the  wall  temperature 
measurement)  and  how  much  is  an  actual  flow  induced  phenomenon.  A  comparison 
of  the  relative  effects  of  the  three  inserts  for  a  given  tube  is  more  accurate,  however, 
due  co  tne  fact  that  any  uncertainty  for  a  given  tube  wall  thermocouple  attachment 
combination  would  be  the  same. 

To  determine  the  degree  of  tubeside  heat  transfer  enhancement  achieved  by 
the  use  of  the  twisted  tape  and  HEATEX  inserts  under  laminar  flow,  an 
enhancement  ratio  of  the  mean  Nusselt  number  for  the  insert  condition  (Nu^,,)  to 
that  for  the  no  insert  condition  (NuMjlBert)  was  calculated.  The  enhancement  ratio  is 
therefore  defined  as: 


Nu. 


Nu 


(6.10) 


no  insert 


Figures  6.31  through  6.33  present  mean  Nusselt  number  as  a  function  of  Reynolds 
number  under  laminar  flow  conditions  for  each  tube  under  the  three  insert 
conditions.  For  the  purposes  of  this  study,  enhancement  ratios  were  calculated  at  the 
two  extreme  Reynolds  numbers  of  200  and  2300  with  a  constant  Prandtl  number  of 
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70  using  the  correlations  developed  in  the  previous  sections.  Table  6.2  summarizes  the 
enhancement  ratios  for  each  tube  and  insert  condition. 

It  is  evident  from  these  figures  that  the  greatest  enhancement  is  achieved  using  the 
HEATEX  insert.  An  enhancement  factor  of  over  seven  was  achieved  for  the  smooth  tube 
using  the  HEATEX  insert  at  Re  =  2300.  This  reduced  to  about  2.6  at  low  Re. 
Significantly  smaller  enhancements  were  achieved  for  the  twisted  tape  inset  for  the 
smooth  tube,  indicating  that  for  laminar  flow,  the  stripping  and  mixing  effect  of  the 
HEATEX  element  is  considerably  more  effective  than  the  simple  swirl  flow  effect 
provided  by  the  twisted  tape.  It  should  be  mentioned  here  that  HEATEX  was  developed 
specifically  with  laminar  flow  enhancement  in  mind.  In  turbulent  flow  (Re  >  10*),  the 
enhancements  obtained  are  much  smaller  (around  1.5  to  2)  and  both  twisted  tape  and 
HEATEX  elements  give  similar  enhancement  values. 

TABLE  6.2  SUMMARY  OF  LAMINAR  FLOW  ENHANCEMENT  RATIOS 
(defined  by  equation  6.10) 


SMOOTH  TUBE 

FINNED  TUBE 

KORODENSE  TUBE 

TWISTED 

TAPE 

HEATEX 

TWISTED 

TAPE 

HEATEX 

TWISTED 

TAPE 

HEATEX 

Low  Re 
(Re  >200) 

1.79 

2.66 

1.49 

2.86 

1.92 

2.58 

High  Re 
(Re  >2300) 

4.04 

7.19 

3.46 

6.71 

2.93 

4.25 
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For  the  finned  and  Korodense  tubes,  the  HEATEX  element  also  provides  the 
highest  enhancements.  At  high  Re,  the  enhancements  provided  by  the  twisted  tape 
and  HEATEX  inserts  are  lower  than  those  obtained  with  the  smooth  tube  for  the 
reasons  mentioned  previously.  However,  at  low  Re,  the  enhancements  for  the  three 
tubes  are  about  the  same  for  a  given  insert,  indicating  the  "stagnant"  regions 
mentioned  earlier  have  less  of  an  effect  on  the  total  heat  transfer. 
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Figure  6.31  Insert  condition  comparison  for  laminar  flow  in  a  smooth  tube. 
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Figure  6.32  Insert  condition  comparison  for  laminar  flow  in  a  Cu/Ni  externally 
finned  tube. 
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FigUr'  6  33  condition  comparison 
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G.  REPROCESSED  REFRIGERANT  CONDENSATION  DATA 

Figure  6.34  shows  the  condensation  data  of  Mazzone  [Ref.  6]  with  the  outside  heat 
transfer  coefficient  recalculated  using  the  present  correlations.  The  data  presented  are  for 
the  top  copper  smooth  tube.  Similar  recalculations  can  also  be  carried  out  for  the  other 
two  tubes  (Cu/Ni  finned  and  Korodense  tubes)  as  well  as  for  condensation  bundle  data. 
The  data  of  Mazzone  are  reprocessed  using  the  Modified  Wilson  '  .*■  echnique  (see 
Mazzone  [Ref.  6]  for  details).  This  technique  essentially  lets  the  leading  coefficient  of 
the  correlation  ’float’  to  try  to  smooth  ou  in  mcertainty  in  the  correlations  being  used. 
As  previously  stated,  the  chosen  correlation  for  the  no  insert  condition  (even  when 
reprocessing  using  the  Wilson  Plot  Technique)  gave  values  of  inside  thermal  resistances 
greater  than  the  measured  overall  thermal  resistance,  thereby  yielding  negative  values  of 
the  outside  heat  transfer  coefficient,  h°.  This  explains  why  Figure  6.34  has  no  data  for 
Mazzone’s  no  insert  condition. 

It  can  be  seen  that  by  using  the  present  developed  correlations  for  the  no  insert 
condition,  the  calculated  values  of  h,  are  significantly  improved  (i.e.  positive  values). 
In  addition,  for  the  other  two  inserts,  the  new  correlations  increase  h„  by  about  10%  for 
the  twisted  tape  insert  and  nearly  20%  for  the  HEATEX  insert  The  newly  calculated 
values  of  h„  when  using  these  inserts  are  now  very  close  and  just  above  the  Nusselt 
prediction  as  may  be  expected  for  low  vapor  velocity  film  condensation.  Mazzone 
commented  that  his  data  were  below  Nusselt  by  about  10% 
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and  thought  that  this  was  a  consequence  of  the  length  of  his  condenser  tubes  violating  the 
uniform  wall  temperature  assumption  of  the  Nusselt  theory.  The  fact  that  these  data  are 
so  close  gives  confidence  in  the  new  correlations,  as  the  insert  itself  should  have  no 
influence  on  the  magnitude  of  tv  By  significantly  reducing  the  inside  resistance,  the 
insert  simply  increase  the  accuracy  of  lv  Of  course,  the  no  insert  condition  should  also 
give  similar  values  of  h,.  Figure  6.34  shows  that  they  are  too  high,  indicative  of  the  fact 
that  the  correlation  for  the  no  insert  condition,  although  much  improved,  is  still  not  ideal. 
If  one  recalls  the  scatter  obtained  in  the  present  data  for  the  no  insert  case,  together  with 
the  fact  that  the  inside  resistance  is  still  dominating,  then  the  results  shown  in  Figure 
6.34  a re  not  too  surprising  and  still  very  encouraging.  It  is  recommended,  therefore,  that 
all  future  refrigerant  condensation  experiments  carried  out  on  the  bundle  facility  should 
be  conducted  with  the  inserts  using  the  newly  formulated  correlations. 

One  final  note  is  the  fact  that  the  present  correlations  were  formulated 
using  a  constant  wall  heat  flux,  whereas  for  the  condensation  experiments,  one 
should  ideally  use  correlations  developed  under  uniform  wall  temperature  conditions. 
One  may  expect  that  the  inside  heat  transfer  coefficient  would  increase  slightly  under 
constant  wall  heat  flux  conditions,  thereby  decreasing  the  outside  heat  transfer 
coefficient.  However,  from  Figure  6.34  it  is  clear  that  the  use  of  the  constant  wall 
heat  flux  correlations  seems  to  increase  h,.  Thus  it  appears  that  with  the  inserts,  the 
inside  resistance  is  so  reduced  that  the  boundary  condition  imposed  has  little 
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Outside  heat  transfer  coeff./(W/m2K) 


influence.  Future  investigators  should  look  into  this  boundary  value  problem  in 
more  detail. 


Smooth  Tube  Condensation  Data 


Coolant  Velocity/(m/s) 


Figure  6.34  Reprocessed  refrigerant  condensation  data  using  new  correlations. 
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VII.  CONCLUSIONS 


The  following  conclusions  are  based  upon  results  obtained  in  this  study. 


A.  GENERAL 


1.  A  test  apparatus  was  successfully  designed  and  built  to  accurately  determine 
inside  heat  transfer  coefficients  for  the  same  tubes  and  insert  conditions  used 
in  the  condensation  experiments. 

2.  Correlations  to  accurately  determine  inside  heat  transfer  coefficients  under 
laminar  flow  conditions  both  with  and  without  inserts  were  successfully 
developed. 

3.  Reprocessing  of  Mazzone’s  [Ref.  6]  refrigerant  condensation  data  using  these 
new  correlations  has  yielded  more  accurate  values  of  the  outside  heat  transfer 
coefficient  (Note:  the  name  of  the  program  used  for  the  reprocessing  is 
DRPCON7.) 


B.  SPECIFIC 


1.  A  departure  from  laminar  flow  into  a  transitional  flow  exists  at  Reynolds 
numbers  of  approximately  2300  for  the  smooth  and  Korodense  tubes  and  2700 
for  the  Cu/Ni  finned  tube.  This  contradicts  Mazzone’s  [Ref.  6]  assumption  of 
laminar  flow  up  to  a  Reynolds  number  of  approximately  4000. 

2.  Laminar  flow  results  for  the  no  insert  condition  were  in  close  agreement  with 
the  Shah  and  London  [Ref.  9]  correlation  for  a  constant  wall  heat  flux 
boundary  condition  for  all  three  tubes  tested. 

3.  The  Hausen  [Ref.  10]  correlation  used  by  Mazzone  [Ref.  6]  significantly 
overpredicted  the  inside  thermal  resistance  and  led  to  nonsensical  results  for 
bo- 
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4.  The  smooth  tube/twisted  tape  data  were  in  close  agreement  with  the  Hong 
and  Bergles  [Ref.  13]  correlation.  Since  this  correlation  was  used  by  Mazzone 
[Ref.  6],  his  results  for  outside  heat  transfer  coefficient  should  be  fairly 
accurate  for  this  tube. 

5.  For  the  Cu/Ni  finned  and  Korodense  tubes  with  the  twisted  tape  insert,  the 
Hong  and  Bergles  correlation  predicted  higher  inside  heat  transfer  coefficients 
than  the  experimental  data  over  the  entire  Reynolds  number  range  tested. 
This  could  have  been  due  to  a  different  twist  ratio  (for  the  Cu/Ni  finned  tube) 
and  stagnant  regions  (for  the  Korodense  tube). 

6.  Experimental  data  for  all  three  tubes  using  the  HEATEX  insert  yielded 
•  results  that  were  well  below  the  correlations  developed  and  used  by  Mazzone 

[Ref.  6]  over  the  entire  Reynolds  number  range  tested. 

7.  HEATEX  inserts  provided  the  best  inside  heat  transfer  enhancement  for  all 
three  tube.  Enhancement  ratios  of  over  a  factor  of  seven  were  found. 
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VIII.  RECOMMENDATIONS 


Based  on  the  results  of  this  study,  the  following  recommendations  are  made: 


1.  Further  refinement  of  the  correlations  to  better  account  for  different  coolant 
inlet  temperatures  should  be  made  by  optimizing  the  Prandtl  number 
exponent 

2.  These  correlations  should  be  used  to  reprocess  all  refrigerant  condensation 
data  and  should  be  used  in  future  refrigerant  condensation  experiments 
conducted  at  NPS.  Note  that  since  these  correlations  were  developed  by 
omitting  coolan.  velocities  of  0.2  and  0.3  m/s,  then  the  condensation  data 
should  do  the  same. 

3.  HEATEX  inserts  provide  the  best  enhancement  of  tubeside  heat  transfer. 
Due  to  the  data  reduction  technique  used,  HEATEX  inserts  should  therefore 
be  used  in  all  tubeside  refrigerant  condensation  experiments  to  increase  the 
accuracy  of  the  outside  heat  transfer  coefficient 

4.  Transitional  flow  correlations  can  also  be  developed  for  the  no  insert 
condition  since  condensation  experiments  typically  move  into  this  flow  region. 
However,  since  the  transition  region  is  not  generally  well  developed,  any 
developed  correlations  for  this  region  will  be  useful  only  for  this  particular 
apparatus. 
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APPENDIX  A.  FLOW  METER  CALIBRATION 


The  flowmeter  used  to  measure  flow  through  the  instrumented  tube  was 
initially  calibrated  at  0°C  and  24°C  by  Mazzone  [Ref.  6].  However  an  additional 
calibration  temperature  of  -15°C  was  needed  for  this  study.  The  flowmeter  was 
calibrated  using  the  same  procedure  as  Mazzone  and  is  summarized  as  follows. 

The  outlet  of  the  flowmeter  was  disconnected  from  the  test  section  and 
connected  to  an  empty  55  gallon  drum.  The  drum  was  placed  on  a  scale  with  a 
resolution  of  0.5  lbs.  The  coolant  pump  was  started  with  the  flow  through  a  bypass. 
The  flowmeter  was  initially  opened  to  maximum  flow  and  then  throttled  back  to 
10%.  The  weight  of  fluid  in  the  drum  was  recorded  and  the  timer  started.  The  weight 
was  then  recorded  Ll  the  end  of  the  calibration  time.  This  procedure  was  repeated 
in  10%  increments  up  to  100%  flow.  The  calibration  run  was  then  repeated  from  the 
maximum  flow  back  to  10%  flow  in  10%  increments.  The  mass  flow  rate  was  then 
calculated  in  kg/s  based  on  an  average  of  the  two  readings.  The  data  for  the 
flowmeter  calibrations  as  a  function  of  nominal  flow  rate  are  shown  in  Table  A.  1. 

A  least  squares  linear  regression  was  performed  for  the  flowmeter  at  each  of 
the  three  temperatures  for  which  the  calibration  runs  were  carried  out  The  linear 
regression  equation  takes  the  following  form: 
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rit  *  a+bN  (^*1) 

where  a  is  the  y-intercept  in  kg/s,  b  is  the  slope  in  kg/s/%  and  N  is  the  flowmeter 
setting  in  %.  The  results  are  shown  in  Table  A.2. 


TABLE  A.1.  MASS  FLOW  RATE  (KG/S)  AT  -15,  0  AND  24°C 


]  FLOW  METER  A 

NOMINAL  FLOW 

MASS  FLOW  RATE  (kg/s)  1 

RATE 

(%) 

-15°C 

o°c 

24°C 

10 

1 

0.0172 

20 

1 

0.0403 

30 

0.0483 

Si 

0.0621 

40 

0.0679 

0.0794 

0.0851 

50 

0.0866 

0.1003 

0.1120 

60 

0.1078 

0.1179 

0.1339 

70 

0.1258 

0.1388 

0.1556 

80 

0.1453 

0.1563 

0.1769 

90 

0.1613 

0.1846 

0.2028 

100 

0.1824 

0.2019 

0.2151 

TABLE  A-2.  FLOWMETER  CALIBRATION  REGRESSION  RESULTS 


FLOW  METER  A 

-15*C 

o°c 

24°C 

'  Slope 

0.0019 

0.00206 

0.00225 

Y-Intercept 

-0.0098 

-0.00414 

-0.00329 

In  order  to  obtain  flowrates  at  any  intermediate  temperatures,  a  simple 


interpolation  procedure  was  incorporated  into  the  data  reduction  program 
"DRPSING."  The  results  of  the  calibrations  for  this  flowmeter  are  shown  in  Figure 
A.l. 
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FLOWMETER 

CALIBRATION  @  -15,  0,  24C 


Figure  A.l  Flowmeter  Calibration. 
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APPENDIX  B.  HEAT  LOSS  TO  AMBIENT  EXPERIMENT 


A  series  of  heat  loss  experiments  were  conducted  to  account  for  the  amount 
of  heat  lost  to  ambient  To  accomplish  this,  the  coolant  was  drained  from  the  test 
tube  and  the  inlet  and  outlet  capped  off  to  prevent  any  air  flow  through  the  tube. 
Heater  power  was  set  at  a  desired  level  and  wall  temperatures  monitored  by  the 
DAU.  Average  wall  temperatures  were  determined  every  minute  for  the  first  15 
minutes  and  then  every  four  minutes  until  an  equilibrium  temperature  was  reached 
(after  approximately  60  minutes).  At  that  point,  any  additional  heat  added  could  be 
considered  lost  to  ambient  This  procedure  was  repeated  at  five  nominal  power 
settings,  1W,  3W,  5W,  7W,  and  9W.  All  data  were  taken  at  an  ambient  temperature 
of  22.0°C. 

Figure  B.l  is  a  plot  of  average  wall  temperature  versus  time  for  the  1W  and 
7W  power  settings.  It  can  be  seen  that  an  equilibrium  average  wall  temperature  was 
reached  after  approximately  60  minutes  for  these  two  power  settings  (as  it  was  for 
the  other  three  power  settings.)  This  graph  shows  that  if  the  average  wall 
temperature  in  the  experiments  reached  46°C  or  74°  C,  then  the  heat  lost  to  the 
atmosphere  was  1W  or  7W,  respectively. 

Table  B.l  shows  the  data  for  the  equilibrium  average  wall  temperature  as  a 
function  of  the  actual  power  setting  (qUp€).  The  data  are  shown  graphically  in  Figure 
B.2.  A  least  squares  linear  regression  was  performed  and  the  resulting  equation  was: 
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_  (7V-44.15)°C  (B.l) 

Hbm  4.39  wrc 

Therefore,  for  a  given  average  wall  temperature,  the  amount  of  heat  lost  to 
ambient  could  be  determined.  Equation  B.l  was  incorporated  into  the  data  reduction 
program  "DRPSING." 


TABLE  B.1  EQUILIBRIUM  AVERAGE  WALL  TEMPERATURE  RESULTS 


NOMINAL  POWER 

ACTUAL  POWER  (W) 

EQUILIBRIUM  T.  (°C) 

1W 

0.9412 

46.47 

3W 

3.122 

59.23 

5W 

5.126 

67.90 

7W 

6.790 

74.34 

9W 

9.560 

84.85 
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HEAT  LOSS  TO  AMBIENT 

Temperature  vs  Time 
1W  and  7W 


Temperature  (C) 


0  20  40  60  80  100 


Time  (min) 

Tamb  •  22.0 

Figure  B.l  Heat  loss  to  ambient;  Average  wall  temperature  versus  time  for  1W 
and  7W. 


118 


HEAT  LOSS  TO  AMBIENT 

Temperature  vs  Heat  Input 


Temperature  (C) 


Tamb  •  22:0 

Figure  B.2  Heat  loss  to  ambient;  Average  wall  temperature  versus  heat  input. 
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APPENDIX  C.  FRICTIONAL/HEAT  LEAKAGE  EXPERIMENT 


As  coolant  flows  through  the  thermocouple  wells/mixing  box  and  along  the 
tube,  a  small  temperature  rise  results  from  frictional  dissipation.  When  the  coolant 
inlet  temperature  drops  below  ambient,  an  additional  temperature  rise  can  be 
attributed  to  heat  leakage  into  the  coolant  from  ambient  Since  there  was  no  means 
available  to  determine  these  terms  separately,  they  were  considered  as  one  heat 
leakage  term  (q^J. 

To  determine  this  heat  leakage  term,  the  power  to  the  heating  elements  was 
secured  and  coolant  passed  through  the  tube  being  tested  at  a  given  flowmeter 
setting.  Any  temperature  rise  detected  was  therefore  due  to  frictional  effects  or  heat 
leakage  from  ambient  This  temperature  rise  was  determined  at  flowmeter  settings 
of  15%  and  100%  for  each  tube  and  insert  condition  at  all  four  nominal  inlet 
temperatures. 

Data  showed  that  flowmeter  setting,  tube  and  insert  condition  had  no 
substantial  effect  on  the  temperature  rise  and  that  only  the  coolant  inlet  temperature 
made  any  difference.  It  cah  therefore  be  concluded  that  the  heat' leakage  term  was 
primarily  due  to  heat  addition  from  ambient  The  results  are  shown  in  Table  C.l. 
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TABLE  C.1  HEAT  LEAKAGE  RESULTS 


INLET  TEMPERATURE 

TEMPERATURE  RISE 

+20°C 

o.o°c 

+  10°C 

o.orc 

0°C 

0.08°C 

-10°C 

0.18°C 

These  corrections  were  incorporated  into  the  data  reduction  program 
"DRPSING'’  by  directly  subtracting  this  temperature  rise  from  the  measured  outlet 
temperature. 
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APPENDIX  D.  POWER  CALIBRATION 


Actual  heating  element  power  (and  therefore  heat  flux)  was  based  upon  the 
applied  voltage  and  total  heating  element  resistance.  To  measure  voltage  and 
current,  the  Data  Acquisition  Unit  made  use  of  a  voltage  and  current  sensing  unit 
However,  a  power  factor  needed  to  be  applied  to  these  sensed  values  to  determine 
actual  values  required  for  the  power/heat  flux  calculation. 

Voltage  and  current  meters  were  placed  in  the  circuit  downstream  of  the  power 
source  and  prior  to  the  first  heating  element  as  shown  in  Figure  D.l.  The  voltage 
rheostat  was  varied  over  its  entire  range  and  voltage  and  current  sensed  by  the 
meters  and  in  the  DAU  were  recorded  (Table  D.l).  For  each  setting,  actual  power 
was  calculated  by: 

P  =  VI  (D.l) 

and  total  heating  element  resistance  by: 

R  =  ~  .  JL  (D.2) 

IV2 

In  addition,  voltage  (PFV)  and  current  (PF,)  power  factors  were  determined  by: 
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PFV  *  Vacwal  and  PF,  =  Iacwal  (D.3) 

9  ^CH2S 

where  the  subscripts  "CH25"  and  "CH29"  refer  to  the  DAU  sensing  unit  channel 
assignments  and  "actual"  to  the  meter  readings.  All  the  above  readings  and 
calculations  were  done  for  five  rheostat  settings  as  shown  in  Table  D.l. 

An  average  heater  resistance  (R),  current  (PF,)  and  voltage  (PFV)  power  factors 
were  then  calculated.  Heater  power  could  then  be  calculated  by  one  of  the  following 
forms  of  equation: 

P  .  (PF,m.Waa>  'an,  <D  4> 

or 

p  >  <D-5) 

The  current  power  factor  varied  too  much  to  be  utilized  in  a  simple  calculation. 
Therefore,  it  was  decided  to  make  use  of  equation  (D.5).  The  final  equation  to 
determine  the  average  heating  element  power  (q^  in  Watts)  for  a  given  rheostat 
setting  was: 

V  =[59'66  Vamf'*9-6  <D-6) 

This  equation  was  incorporated  into  the  data  reduction  program  "DRPSING." 
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TABLE  D.1  POWER  CALIBRATION  RESULTS 


POWER 

(W) 

V 

(VOLTS) 

I 

(AMPS) 

R 

(OHMS) 

VCBn 

(VOLTS) 

less 

(AMPS) 

PFV 

PF, 

3131 

30.7 

0.61 

5032 

03135 

0.0106 

59.78 

5730 

76.14 

61.4 

124 

4931 

1.032 

0.0430 

59.49 

28.80 

30330 

1223 

248 

4939 

2054 

0.1050 

59.63 

23.60 

75036 

1924 

3.90 

4930 

3226 

0.1746 

59.64 

2230 

113.84 

234.0 

4.76 

4931 

3.915 

02170 

59.77 

21.90 

POWER  CALIBRATION 
ELECTRICAL  SCHEMATIC 


Figure  D.l  Electrical  schematic  of  power  calibration  experiment 


APPENDIX  E.  TUBE  AND  INSERT  SPECIFICATIONS 


TABLE  E.1  TUBE  CHARACTERISTICS/SPECIFICATIONS 


PHYSICAL 

CHARACTERISTICS 


Inside  Diameter  (D;) 


Outside  Diameter  (D0) 


T/C  Diameter  (Dr/C)  t/c 


Wall  Thickness 


Length 


Overall  Length  (L) 


Heated  Length  (LH) 


Heated  Finned  Length 

(Lf) 


Smooth  Section  Length 

(Ls) 


T/C  Location  (1) 


Pos  1 


Pos  2 


Pos  3 


Pos  4 


T/C  X/D  for  Lh  (2) 


Pos  1 


Pos  2 


Pos  3 


Pos  4 


Inside  Surface  Area  for  LH 


A, 


T/C  Attachment  Method 


Cu 

SMOOTH 


13.26mm 


15.88mm 


14.36mm 


1.31mm 


152.4cm 


140.6cm 


15.2cm 


55.9cm 


96.5cm 


137.2cm 


6.77 


37.41 


68.06 


98.71 


0.0586m2 


Silver 

Solder 


Cu/Nl 

FINNED 


10.16mm 


14.00mm 


12.62mm 


1.22mm 


153.7cm 


143.0cm 


122.3cm 


15.7cm 


17.8cm 


57.2cm 


96.5cm 


135.9cm 


2.05 


40.80 


79.55 


118.27 


0.0340m2  (3) 


Soft  Solder 


Cu/Ni 

KORODENSE 


13.35mm 


15.85mm 


14.53mm 


1.25  mm 


152.7cm 


142.4cm 


15.6cm 


56.2cm 


96.8cm 


134.5cm 


7.91 


38.35 


68.79 


96.99 


0.0600m2 


Tach  Weld / 
Liquid  Metal 


measured  trom  tube  inlet; 
finned  heated  length. 


r«i 


measured  trom  start  ot  heated  length; 
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TABLE  E.2  INSERT  CHARACTERISTICS/SPECIFICATIONS 


I  TWISTED  TAPE 

NARROW 

WIDE 

thickness,  S 

0.559mm 

0.559mm 

width,  w 

10.16mm 

13.26mm 

180°  pitch,  H 

3.05cm 

5.30cm 

twist  ratio,  y 

3 

4 

material 

brass 

brass 

HEATEX 

SMALL 

LARGE 

diameter 

10.16mm 

13.26mm 

material 

stainless  steel 

stainless  steel 
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APPENDIX  F.  UNCERTAINTY  ANALYSIS  FOR  THE  INSIDE 
HEAT  TRANSFER  COEFFICIENT 

Data  file  SMNI03  (smooth  tube,  no  insert,  +10°C/10kW/m2),was  chosen  to 
conduct  an  uncertainty  analysis  of  the  experimental  determination  of  the  inside  mean 
heat  transfer  coefficient,  hj.  In  order  to  assess  the  effects  of  coolant  velocity  on  the 
uncertainty  in  the  measurement  of  hj,  this  uncertainty  analysis  was  performed  at 
coolant  velocities  of  0.4  and  1.4  m/s  (for  the  same  data  set).  The  calculations  below 
are  shown  for  the  0.4  m/s  coolant  velocity.  Table  F.l  summarizes  the  results  of  the 
uncertainty  analysis  for  both  0.4  m/s  and  1.4  m/s  coolant  velocities. 

The  uncertainty  analysis  utilizes  the  procedure  suggested  by.  Kline  and 
McClintock  [Ref.  24].  This  procedure  states  that  if 

R=R  . •*) 

then  the  uncertainty  in  R,  namely,  SR  is  given  by  the  equation: 


where  x,  is  the  measured  variable  and  Sx„  is  the  uncertainty  in  the  measured  variable. 

The  uncertainty  in  the  coolant  temperature  measurements  is  due  to  uncertainty 
in  the  voltage  measured  by  the  thermocouples.  For  the  thermocouples  used,  this 
value  was  assumed  to  be  4  microvolts  or  0. 1°C.  The  uncertainty  in  A T  is  found  from 
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the  uncertainty  in  the  coolant  inlet  and  outlet  temperatures,  Ta  and  and  can  be 
calculated  as  follows: 


T  m  T(14)+T(15) 
-  2 


Taking  the  partial  derivatives  we  obtain: 

3T«  _  1 
37(14)  37(15)  2 

The  uncertainty  in  the  measurement  of  7^  is  thus: 


37. 


dT 


l  37(14) 


- ST 


dT 


{  37(15) 


-57 


<u 


Based  on  the  thermocouple  uncertainty,  ST  *  O.TC,  we  obtain: 


ST  =  0.071  °C 

CO 


The  coolant  inlet  temperature  is  calculated  from  a  single  thermocouple.  The 
uncertainty  in  this  measurement  is  calculated  as: 


37V 

37(13) 


1 


and 


then  the  uncertainty  in  AT  is  calculated  from: 


Mr  *  [(STJ^fSTjf 


or, 


SAT  *  0.1225°C 

For  a  measured  AT  of  3.7°C,  the  uncertainty  is  therefore: 


_  =  0.0332  =  3.32 % 
AT 
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Coolant  velocity  is  calculated  from: 


m 


For  a  coolant  velocity  0.4  m/s,  m  was  calculated  to  be  0.05954  kg/s.  The  uncertainty 
in  the  reading  of  the  flowmeter  is  defined  as  the  scale  interpolation.  This  value 
corresponds  to  one-half  the  value  of  the  smallest  marked  increment  In  this  case,  the 
scale  interpolation  factor,  &&  is  0.5%,  which  corresponds  to  0.00101  kg/s  resolution. 
In  addition,  a  time-wise  jitter,  of  0.5%  was  noted  during  the  runs.  Therefore,  the 
uncertainty  in  m  is  given  by: 


8m 

-  -  3 

» 

m 

Lb*  i 

l*  *)J 

For  a  the  measured  coolant  mass  flow  rate: 


—  *  0.0239  =  239% 
m 


The  uncertainty  in  the  cross-sectional  area  of  the  tubes  was  estimated  from: 


H 

Ac 


' 8D '2 


2 


where  5D;  is  given  as  0.1  mm  based  on  tolerances  supplied  from  the  manufacturer. 
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temperature  rise  where: 


*± 

4 


(8rit\2  /SA7V 

N*M. 


or 


H  *  0.0409  =  4.09% 

4 


The*  uncertainty  in  the  heat  flux  is  then: 


where  the  uncertainty  in  As  is  given  by: 


'  sd' 


\D‘, 


2 

+ 


OlS 


Here,  the  uncertainty  in  the  length  of  the  tube,  SL,  is  assumed  to  be  5  mm  while  8D„ 
is  assumed  to  be  0.1  mm  (again  from  given  manufactures  tolerances).  Substitution 
of  the  appropriate  numbers  yields: 


0.0147  *  1.47% 
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Finally,  since  h;  is  given  by: 


1  LMTD 


the  uncertainty  in  the  LMTD  must  be  calculated.  The  LMTD  is  given  by: 


LMTD  = 


T  -T  . 
In  .  JL  — 

t^-t 


In  order  to  calculate  the  uncertainty  in  the  measurement  of  the  LMTD,  the  following 
equations  were  used: 


dLMTD  =. _ 1 

3AT  r  t 


T  -7 
In  *  ° 

T  -T 


dLMTD 

dT. 


In  LlL-  (T.-Ta)(T„-TJ 

1  m  1  co 


dLMTD 


T  -T 

In  ■■  — q-  (T  -T  ) 

m  co* 
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dLMTD 

*T* 


A  T 


In 


T  -T 

m  a 

T  -T 

m  cm 


(7>7V) 


The  uncertainty  in  the  LMTD  is  therefore  given  by: 


8 LMTD  = 


1  dLMTD  hat) 

2 

+ 

dLMTD  aT  ' 

d&T 

dT  . 

_  V  > 

1  “•  / 

dLMTD 

ar 


2 

+ 


dLMTD 

dT 

a 


\ 

ST 

a 

/ 


2 


Before  the  uncertainty  in  LMTD  can  be  calculated,  the  uncertainty  in  the 
measured  average  inside  wall  temperature  must  be  determined.  The  average  inside 
wall  was  determined  from  a  total  of  twelve  thermocouples,  three  at  each  of  the  four 
longitudinal  positions.  The  uncertainty  in  a  single  wall  temperature  measurement  is 


due  to  the  uncertainty  in  the  voltage  measurement  of  the  thermocouple  (5T  =  0.1°C) 
and  the  uncertainty  due  to  the  influence  caused  by  the  various  attachment  techniques 
used.  Obviously,  the  variations  seen  in  the  wall  thermocouple  readings  were  due  to 
this  latter  uncertainty  and  some  account  has  to  be  made  of  this  as  this  is  by  far  the 
dominant  uncertainty.  If  the  average  temperature  at  a  given  longitudinal  position  is 
assumed  to  be  the  required  temperature,  then  the  standard  deviation  of  the  three 
thermocouple  readings  at  that  position  could  be  assumed  to  be  this  latter  uncertainty. 
Therefore,  the  uncertainty  at  a  given  longitudinal  position  can  be  represented  by: 


i 

7 


134 


135 


Considering  position  1,  the  average  inside  wall  temperature  is: 


r  .  r(io)>r(ii)-fT(i2) 

f*  3 

For  the  three  thermocouple  readings  of 

T(10)  «  25.99*C,  T(ll)  =  25.74*C  and  T(12)  =  26.01*C, 
the  average  inside  wall  temperature  for  position  1  is 

Tpo*  -  25-91*  C 

The  standard  deviation  at  this  location  is  given  by: 

S.D.  «  -  (ECr.-T^Inf*  -  0.123* C 

Therefore, 

87* 

— =£  »  4.7  x  10-5  -  0.47* 

T~l 

Repeating  this  for  the  other  three  locations  gives: 

57^  -  0.497° C 

*  0.634°  C 

57^  -  0.145°C 
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Therefore,  the  uncertainty  in  the  measurement  of  average  inside  wall  temperatures 


is: 


S7\  =|(0.25)(0.123)f  +[(0.25)(0.497)f  +[(0.25)(0.634)f  +[(0.25)(0. 145)f =  0.207 

After  substituting  the  appropriate  values  for  dLMTD/dAT,  dLMTD/dT^, 
dLMTD/dT,,,,  and  dLMTD/dT^  the  uncertainty  in  the  measurement  of  LMTD  is: 


8LMTD  =  1.122°C 


and 


SLMTD 

LMTD 


1.122 

24.13 


0.0465  =  4.65% 


The  uncertainty  in  the  overall  heat  transfer  coefficient,  U&  is  then  given  by: 


or 


/  SLMTD  \ 
[  LMTD  J 


i 

7 


=  [(0.0417)2+(0.0463)2]T 
h. 


Sh. 

* 


0.0625 


6.25% 
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__L  =  0.0625  *  6.25% 

Table  F.l  summarizes  the  results  of  the  uncertainty  analysis  for  the  smooth  copper 
tube  with  no  insert  for  coolant  velocities  of  0.4  and  1.4  m/s,  respectively  at  a  nominal 
inlet  coolant  temperature  of  +10°C. 

The  uncertainty  analysis  neglects  any  uncertainty  in  the  physical  properties  of 
the  coolant  For  the  most  part  this  is  probably  a  reasonable  assumption.  The 
exception  to  this  may  be  the  viscosity  due  to  its  sensitivity  to  temperature.  Hence,  the 
uncertainty  reported  here  is  likely  to  be  conservative. 

It  is  curious  that  the  relative  magnitudes  of  the  uncertainty  for  the  low  and  high 
coolant  velocities  are  nearly  the  same.  Two  different  mechanisms  dominate  in  this 
calculated  uncertainty  at  low  and  high  coolant  velocities.  In  terms  of  the  mass  flow 
rate,  the  time-wise  jitter  and  scale  interpolation  tenns  are  the  same  for  both  coolant 
velocities.  However,  the  actual  mass  flow  rate  is  less  at  the  lower  flow  rates.  Hence, 
the  uncertainty  in  the  mass  flow  rate  calculation  is  more  significant  at  the  lower 
coolant  velocities.  At  higher  coolant  flow  rates,  the  uncertainty  in  the  calculation  of 
LMTD  becomes  the  dominant  term.  This  is  because  the  coolant  temperature  rise  is 
smaller  at  high  coolant  velocities  but  the  uncertainty  in  the  temperature 
measurements  remains  the  same. 


TABLE  F.l  SUMMARY  OF  RESULTS  FOR  SAMPLE  UNCERTAINTY  ANALYSIS 


Variable 

Ve  =  0.4  m/s 

Vc  =  1.4  m/s 

We 

Td 

0.98% 

0.87% 

5  T 

CO 

0.51% 

0.56% 

T 

CO 

SAT 

A  T 

3.32% 

10.67% 

5m 

2.39% 

0.69% 

m 

&4, 

4 

1.07% 

1.07% 

SV 

K 

2.62% 

1.27% 

SRec 

W 

2.32% 

1.47% 

Sq 

4 

4.09% 

10.69% 
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0.57% 

0.71% 

r* 

* 

4.17% 

10.72% 

q'J 

SLMTD 

LMTD 

4.65% 

16.32% 

*1 

*1 


6.25% 


19.53% 


APPENDIX  G.  PROGRAM  LISTING 


im : 

ta. 

ms  ' 

1930 

1835 

1048 

1045 

106S 

1000 

1085 

1000 

1005 

1125 

1130 

1135 

1140 

1145 

1150 

1155 

1160 

1170 

1180 

1215 

1220 

122S 

1230 

1235 

1240 

1245 

1250 

1255 

1260 

1261 

1263 

1264 

1265 

1266 

1267 

1268 
1260 

1270 

1271 

1273 

1274 

1275 

1276 

1277 

1283 

1284 

1285 

1286 
1287 
1200 
1300 
1310 
1315 
1325 
1345 
1355 


PukloSE  fSu^progran  eollacts  and  proctsm  »maia  tnjtrunantad  tuba  da 

CREATED :  APRIL  1.  1002  <5.  rtEMORV ) 

SEEP 

PRINTER  IS  1 

PRINT  USIN6  "4X,  "SELECT  OPTION"* 

PRINT  US  INS  "SX,"0  T  AKINS  OATA  OR  R6PR0CESINS  PREVIOUS  DATA" 

PRINT  USINS  "SX,"1  PURSE  FILES'" 

PRINTER  IS  701 
INPUT  leal  1 

IF  IeaIl-0  THEN  CALL  Nam 
IF  Ieall-1  THEN  CALL  Purga 
ENO 

SU0  Ha in 
CON  /Ce/  Cl  7) 

CON  /Fid/  Ift  . 

Km  /Nua/  lm,T»at,Odpl,Hnu»,Kf .Rhof .Hfg.Huf ,Qo, I  tuba 
CON  /Will/  0oa<4>ipia(4),Kna<4),Ia«t,0root(*).LdI31 
COH  /Mtl2/  0alta.I*dt,Naata,Had.Cta<4),Aloaa<3> 

Oin  Enf  1 17  ),Tpt  4  > ,  T (  17  )  ,Ho(  4  ),Qdp(  4  >,Uo<  4  l.Pel  4  I.Xitarl  4  >,Nu*l  4  )  ,Nu»»<  4  ) 
Oin  X**l(4),X»*2(  4 ),Xxx31 4 ),Xxx4( 4  >,X*xS<  4 ) ,X**6(  4  > 

Oin  Nuwl  4  ),Nuag<  4 )  .Lddl  4  )  ,Nu(  4  ) ,  X*<  4 > 

OATA  1.0.5.172,0.226,0.063,5.172 
RE AO  Cial * ) 

OATA  0 . 10086001 . 25727.04360, -767345 . 8205. 7802SS0S. 81 
OATA  -0247486580 , 6 . 07688E 1 1 . -2 . 66 102E1 3 , 3 . 04078E l 4 
READ  C< • » 

OATA  0 .015875 ,0.81 4000 , 0 . 015850 ,0.0150,0.0 
OATA  0.013250,0.010160,0.013350,0.01346.0.0 
OATA  386. 0,42. 08. 42. 08, 21. 0.0-0 
OATA  0.0005588.3 

REAO  Ooal  ♦ I.Olat • >,Knat • l.Oalta, Hod  1  Hod-H/Oi 

IF  I  tuba-0  THEN 

OATA  6.765,37.41.68.06,08.71 

REAO  Ldl •  )  '  LEN6TH/QIAHETER  RATIO  FOR  4  T/C  L0C3  ISHOOTH) 

L- 1.406  !  Haatad  langth  of  snooth  tuba 

ENO  IF 

IF  l tuba- 1  THEN 

OATA  2.05,40.80,70.55,118.20  . 

REAO  Ldl*  I  1  FIN  LEN6TH/QIANETER  RATIO  FOR  4  T/C  LOCS  (CU/NI  FIN) 
L-1.223  1  Haatad  fiAnad  langtn  of  Cu/Ni  f»n  tuba 

ENO  IF 

IF  I  tuba-2  THEN 

OATA  7.008,38.35.63.702,36.087  . 

ftEftO  Ld(.)  .  L£N6TH/0IAnETER  RATIO  FOR  4  T/C  LOCS  IKOROOENSE) 

L-1.422  1  Haatad  langth  of  Korodensa  tube 

ENO  IF 

Jsat-0 

Okaeept-1 

BEEP 

INPUT  -ENTER  aonth,  DATE  and  TINE  rnn:OO:HM.nn:SS>".Otg0 
OUTPUT  703;-TO".Otg0 
8EEP 
I hard- 1 

INPUT  "WANT  A  HAROCOPV  PRINTOUT  < 1 -OEF-/ES , 0-N >  .Ihard 

SEEP 

Iin-1 

IF  I  hard* l  THEN  PRINTER  IS  701 

8EEP 
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1360 
1 365 
1370 
1375 
1360 


INPUT  “GIVE  A  NAME  FOR  THE  NEW  OATA  FILE*. File* 

CREATE  BOAT  File*. 20 
ASSISN  ®File  TO  File* 

SEEP 

INPUT  "ENTER  TU8E  CODE  (0-SMOOTH  CU.l-FINNED  CU/NI ,2-KOROOENSE, 3-TIT 


FIN  )" , I  tube 

1383  BEEP 

1384  INPUT 


"0-NO  INSERT,  1-TUISTEO  TAPE  (WIDE),  2-HEATEX  (L>. 


3-HEATEX  (S) 


,  4-TUISTEO  TAPE  < NARROW >'. Insert 


1385 

13B6 

1388 

1389 

1390 

1391 

1392 

1393 
1395 
1400 
1405 
1430 
1435 
1595 
1600 
1605 
1610 
ISIS 
1620 
1625 
1630 
1635 


IF  Insert-1  THEN  Hod-4 
IF  Insert-4  THEN  Hod-3 
BEEP 

INPUT  "NOMINAL  SUMP  TEMP.  0—10  C,  1-0  C,  2—10  C,  3—20  C'.Itenp 
SEEP 

Egret -54 

INPUT  "ENTER  ES  CONCENTRATION  <  UT  PERCENT )  .OEFAULT  -  S4X",£grat 
BEEP 

ENTER  709;OtgS 
OUTPUT  0File:Otg* 

OUTPUT  SFila; Itube, Egrat,0dl,Qd2,0d3,0d4,0dS 
PRINT  USING  "10X.*"FILE  NAME:  *" , 12A* ;FilaS 

PRINT 
lout-1 
BEEP 

INPUT  "WANT  TO  CREATE  AN  OUTPUT  FILE?  < 1-OEF-YES.0-N0 >" , lout 
IF  lout— 1  THEN 
•  BEEP 

INPUT  "ENTER  A  NAME  FOR  OUTPUT  FILE’.FoutS 
CREATE  80AT  FoutS.S 
ASSIGN  0Fout  TO  Pout* 

ENO  IF 


1640  i 

1646 

1647 

1649 

1650 

1651 
1653 
1655 
1660 
1661 

1665 

1666 

1667 

1668 

1670 

1671 

1672 

1673 

1674 

1675 

1676 

1677 
1679 
1679 
1690 
1631 
1682 
1693 
1634 


IF  Itube-0  THEN  Otc-. 01436  1  OIAMETER  TO  T/C  8URIAL  (SMOOTH) 

IF  Itube-l  THEN  Otc-. 0126  1  OIAMETER  TO  T/C  BURIAL  (CU/NI  FIN) 

IF  I  tube-2  THEN  Otc-. 01453  1  OIAMETER  TO  T/C  BURIAL  (KORODENSE) 

Oo-Ooa< Itube ) 

Ool-Ooa( Itube  > 

Qi-Oia< I  tube ) 

Kn-Kwa( Itube ) 

Ax-PI *Ox "2/4  i  Cross-sectional  area 

IF  Itube-2  THEN  A*-. 00014  ■  8A5EQ  ON  WOLVERINE  DESIGN  OATA 


Ao-PI *Oo«L 


Ai-Pr*Oi*L 

IF  Itube-2  THEN  Ao-.05»L 
IF  Itube-2  THEN  Ai-Ao/1.186 
Perin-PI»Oi 

Anarea-Pl*( Do ‘2-Oi*2 > / 4 . 


BASED  ON  UOLVERINE  DESIGN  OATA 
BASED  ON  UOLVERINE  OESIGN  OATA 


Ldd(  l  )-L/<  4.  »Oi  ) 

Ldd( 2  )-L/( 2. »Oi  ) 

Ldd<  3  )«3.  »L/<  4.  »Oi  ) 

Ldd(  *  )-L/Oi 
IF  Itube-0  THEN 

Ulenl-,0627  '  UNHEATEO  INLET  ENO  OF  COPPER  SMOOTH  TUBE  ONLY 

Ulen2* .0572  •  UNHEATED  OUTLET  ENO  OF  COPPER  SMOOTH  TUBE  ONLY 


ENO  IF 

IF  Itube-l  THEN 

Ulenl- .045  '  UNHEATED  INLET  ENO  OF  CU/NI  FIN  TUBE  ONLY 

Ulen2» .06  1  UNHEATED  OUTLET  ENO  OF  CU/NI  FIN  TUBE  ONLV 

Ocuo-. 015875  '  UNHEATED  OUTSIDE  OIAMETER  OF  CU/NI  FIN  TUBE 
Ocui-.0132S8  1  UNHEATED  INSIDE  OIAMETER  OF  CU/NI  FIN  TUBE 
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1685 

1686 

1687 

1688 
1689 
1698 

1691 

1692 

1693 

1694 
1790 

1701 

1702 

1703 

/  9  )? 

1704 

1705 

1706 

1707 

1708 

1709 

1710 

1711 

1712 

1713 

1714 

1715 

1716 

1717 

1718 

1719 

1720 

1721 

1722 

1723 

1724 

1725 

1726 

1727 

1728 

1729 

1730 

1731 

1733 

1734 

1735 

1737 

1738 

1739 

1740 

1741 

1742 

1744 

1745 

1746 

1747 

1748 

1749 
1790 
1751 


lcui-.112  i  HEATEQ  INLET  SMOOTH  ENO  OF  CU/NI  FIN  TueE 
Lcuo-.097  i  HEATEO  OUTLET  SMOOTH  ENO  OF  CU/NI  FIN  TUBE 
Per im»PI»Qcui 

Anarea-PI»<  Gcuo“2“Ocut*2  >/4. 

ENO  IF 

IF  I  tube-2  THEN 

Ulenl-.0S0  '  UNHEATEO  INLET  ENO  OF  CU/NI  KORQOENSE  TUBE  ONLY 
Ulen2-.0S5  '  UNHEATEO  OUTLET  ENO  OF  CU/NI  KORQOENSE  TUBE  ONLY 
ENO  IF 

Hr-50  !  HEATER  RESISTANCE 

! 

PRINTER  IS  1 
SEEP 

INPUT  "WHAT  IS  REQUIREO  NOMINAL  COOLANT  VELOCITY  IN  n/s  (MAX.  -  1.4  n 

' , Vegst 

OUTPUT  709; ”AR  AF13  AL13  VRS" 

OUTPUT  709;* AS  SA* 

Vsum-0 

FOR  J-l  TO  S 
ENTER  709 ;E 
Vsum-Vsum+E 
NEXT  J 

Enf < 13 )“UauM/S 
T( 13>-FNTvav<E*f< 13) ) 

Tavgl«T< 13 ) 

IF  Tavgl<-15  OR  Tavgl>24  THEN 

PRINT  'COOLANT  TEMP.  NOT  IN  RANGE  >15  TO  24  C’ 

GOTO  1703 
ENO  IF 

Rhoegl-FNRhoegf  Tavgl , Egret ) 

Hdot l-«hoegl »Ax«Uegat 
PcalS-tMdotl*.  0098)/.  0019 
Pca0-<  Mdotlf .004 14  )/ . 00206 
Pca24-<  Mdo 1 1  + . 00329 ) / . 00225 
IF  Tavgl <0  THEN 

Pc(0)-PcalS*<  <Tavgl  +  l5)»(Pca0-PcalS)  >/ IS 
ELSE 

Pc<  0  )-Pca0+<  Tavgl»(  Pca24-Pca0  )  )/24 
ENO  IF 

PRINT  "SET  FLOUMETER  READINGS  CORRESPONDING  TO  - 
PRINT  USING  -SX.--Z  OF  METER  A  -  ” ,000.0" ;Pc( 0  ) 

PRINT 

PRINT  "HIT  CONTINUE  UH£N  READY" 

PAUSE 
Repeat:  i 

'  PRINTER  IS  701 

Ido-1 

ON  KEY  <3,15  RECOVER  1734 

PRINT  USING  ”4X , - "SELECT  OPTION 

PRINT  USING  ‘6X  ,  “  "0-TAKE  DATA  *“ 

PRINT  USING  -SX,-* 1-SET  Heater  Tape  Heat  Flux  (DEFAULT)  -” 

PRINT  USING  'SX," "2-CHECK  FLOWMETER  Z  SETTING  •** 

PRINT  USIN6  ' 4X  ,  ” "NOTE ;  KEY  0  -  ESCAPE'” 

1  Ido  •  desired  option 

BEEP 

INPUT  Ido 

1  Set  default  value  for  input 

IF  I  do >2  THEN  Ido*2 
1  Take  data  option 

IF  Ido-0  THEN  1816 
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17S2 

1754  • 

1755 

1756 

1758 

1759 

1761  ' 

1762 

1764 

1765 

1766 

1767 

1768 

1769 

1770 

1771 

1772 

1773 

1784  ! 

1785 
1789 

1793 

1794 

1795 

1796 

1797 

1798 

1799 

1800 
1801 

1803 

1804 

1805 
(  12> 

1806 

1807 
) 

1808 

1809 

) 

1810 
1811 

) 

1812 

1913 

1814 

1315  ' 

1816 

1818 

1819 

1820 
1821 
1822 

1823 

1824 

1825 

1826 

1827 

1828 
1829 


IF  Ido-2  THEN  1703 

Loop  to  set  h«at  flux  and  check  temperatures 
IF  Ido-1  THEN 

INPUT  "ENTER  OESIREO  HEAT  FLUX’ .Gdpuan 

PRINT  USING  -4X,*-  Qdp  (wanted)  Qdp  (actual)  Tm  Qelt’" 
PRINT  USING  “4X , "  (U/m’2>  (U/m’2>  (C>  <C> 

Read  thermocouple  voltages 
OUTPUT  709; -AR  API  AL15  UR5‘ 

FOR  1-1  TO  IS 

OUTPUT  709; ’AS  SA' 

Usum-0 

FOR  J-l  TO  20 
ENTER  709 ;E 
Usum-Usum+E 
NEXT  : 

Emf( 1  >-Vsu m/20 
T< I  )-FNTvsv(Emf< I  ) ) 

NEXT  I 

Compute  average  temperature  of  coolant  outlet 
Tout-<  T(  14  )+T(  IS  )  )/2 
Oelt-Tout-T< 13 ) 

OUTPUT  709; -AR  AF2S  AL2S  URS’ 

OUTPUT  709; ‘AS  SA’ 

ENTER  709; Amp 

OUTPUT  709; -AR  AF29  AL29  URS* 

OUTPUT  709; -AS  SA’ 

ENTER  709;Uolt 
Otape-( 60«Uol t )'2/H- 
Qdptape-Otape/Ao 
PRINT 

PRINT  USING  ’ 4X ,S<  HZ . 3DE , 3X  )* . Qdowan , Qdp tape , T( 15 ) ,0ei t , Q tape 
PRINT 

PRINT  USING  ’4X, "TEMPS  AT  POSITION  1  * ’ , 3( ODO .DO , 3X ) * ; T( 10 ) , T(  1 1 ) , T 
PRINT 

PRINT  USING  *4X , * ’TEMPS  AT  POSITION  2”  ,3< 000 .DO ,3X >* ; T( 7  ) , T( 8 > , Tc 9 
PRINT 

PRINT  USING  -4X, "TEMPS  AT  POSITION  3"  ,3(  ODD  .00 , 3X  )’ ;  T(  4  > ,  T(  S  > ,  Tv  6 
PRINT 

PRINT  USING  ’4X,  "TEMPS  AT  POSITION  4", 3'  00 .00 . 3X  )  ’  .  T<  1  ) .  T(  2  I ,  T<  2 

WAIT  2 
GOTO  1762 
END  IF 

TAKE  OATA  IF  Im-0  LOOP 

Twm-0 

Tul-0 

PRINTER  IS  701 

OUTPUT  709. ’AR  AF 1  AL1S  URS’ 

FOR  l-l  TO  15 

OUTPUT  709; ’AS  SA’ 

Usum-0 

FOR  I-t  TO  20 
ENTER  709, E 
Usum-Usum+E 
NEXT  J 

£mf< I  )«Usum/20 
NEXT  I 
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1840  I 

1845  '  OAT A  ANALYSIS 
13S0  ' 

I860  FOR  1-1  TO  IS 

1865  T< I  )«FNTvsv( Enf  < I ) ) 

1870  NEXT  I 

1871  FOR  1-1  TO  12 

1872  Twi-Tui+T < I  ) 

1873  NEXT  I 

1874  Tuiav-Tui/12. 

1876  Tout-<  T(  14  >+T<  IS  )  )/2 

1880  Oelt-Tout-T< 13) 

1881  Hlc-0.  !  Heat  leakage  correction 

1883  IF  I  tube-0  THEN 

1884  IF  Itenp-0  THEN  Hlc-.18  «  HEAT  LEAKAGE  CORRECTION  AT  -10  C 

1885  IF  Itenp-1  THEN  Hlc«.08  •  HEAT  LEAKAGE  CORRECTION  AT  0  C 

1886  IF  Itenp-2  THEN  Hlc-.07  i  HEAT  LEAKA6E  CORRECTION  AT  *10 

1887  IF  Itenp-3  THEN  Hlc-0.  !  NO  HEAT  LEAKA6E  CORRECTION  AT  + 

1888  Oelt-Oelt-Hlc 

1890  END  IF 

1891  Tout-Tout-Hlc 

1893  OUTPUT  709; *AS  SA“ 

1894  ENTER  709;Anp 

1895  OUTPUT  709: ’AR  AF29  AL29  VRS’ 

1896  OUTPUT  709; "AS  SA* 

1897  ENTER  709;Uolt 

1898  Qtape-<60*wolt )*2'Hr 

1923  Qdptape-Qtape/Ao 

1925  Jsat«Jset*l  — 

1926  IF  Qkaccpt-0  THEN  Jset-Jset-l 

1930  PRINT 

1935  PRINT  USING  ’10X,”Oata  set  nunber  -  ”,DO*;Jset 

1940  PRINT 

1941  IF  I  tube-0  THEN  PRINT  USING  * 10X , * *Snooth  tube’*’ 

1942  IF  I  tube-2  THEN  PRINT  USING  *  10X  , ’’Korodense  tube”* 

1943  IF  Itube-l  THEN  PRINT  USING  ’ 10X ,* ’Copper /Nickel  fin**’ 

1945  IF  Itenp-0  THEN  PRINT  USING  * 10X, ’’Nominal  sunp  temp.  *  -10  C”’ 

1946  IF  Itenp-1  THEN  PRINT  USING  ' 10X , ’’Nominal  sunp  tenp.  -  0  C”’ 

1948  IF  Itenp-2  THEN  PRINT  USING  ’ 10X , ’’Noninal  sunp  tenp.  *  *10  C*” 

1949  IF  Itenp-3  THEN  PRINT  USING  ’ 10X, ’’Noninal  sunp  tenp.  -  *20  C’** 

1951  IF  Insert-0  THEN.  PRINT  USIN6  ’10X,”No  insert  used”* 

1952  IF  Insert-1  THEN  PRINT  USING  *10X,”Uide  twisted  tape  insert  used”’ 

1953  IF  Insert-4  THEN  PRINT  USING  ’ 10X ’Narrow  twisted  tape  insert  used" 

1954  IF  Insert-2  THEN  PRINT  USING  *10X,  "Large  HEATEX  insert  used*** 

1955  IF  Insert-3  THEN  PRINT  USING  ’10X,"Small  HEATEX  insert  used”* 

195G  PRINT  USING  “  10X ,’ "Noninal  set  heat  flux  -  ”  , M  2 . 30E" . Qdp t ape 

1960  PRINT  USING  ' 10X ,’ "Noninal  set  coolant  velocity  -  ' " ,0.0" . Uegst 

1970  PRINT 

1990  Tavg«T< 13  )+0el t • . 5 

1996  Rhoeg-FNRhoeg< Tavg.Egrat > 

2000  Nuueg-FNNueg< Tavg.Egrat  ) 

2005  Hueg-Nuueg*Rhoeg 

2010  Cpeg-FNCoeg( Tavg.Egrat  ) 

2015  Keg-FNKegt Tavg.Egrat  ) 

2020  Preg-Cpeg-Hueg/Keg 

2025  Mdot -FNFncal ( T<l3),Pc(0)> 

2030  Ueg-Hdot/(  Rhoeg»A>.  ) 

2035  Reeg-Ueg»Oi/Nuueg 

2040  fles-4»Mdot / < PI •Hueg* ( 0 i -4 »0el t a > > 

2045  Qdot-Mdot»Cpeg»Oelt 
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r-j  o 


2048 1 
2048 

2050 

2051 

2052 

2053 
2058 

2057 

2058 

2059 

2060 
2061 
2062 

2063 

2064 

2065 

2066 

2067 

2068 

2069 

2070 

2071 

2072 

2073 

2074 

2075 

2076  » 

2077 

2078 

2079 

2080 
2081 
2082 

2083 

2084 

2085 

2086 

2087 

2088 

2089 

2090 

2091 

2092 

2093 


ITERATE  TO  FINO  INSIOE  WALL  TEMPERATURES 
Twl-Tyiav 
Tn2"Tyiav 
Tnl-<Tmi+Tm2>/2. 

Cti»Qdot»LOG<  Otc/Oi  >/<  2.*PI»L»Km  > 

T*2"Tuiav-Ct i 

IF  ABS< Tn2-Tnl >>.001  THEN  SOTO  2051 
FOR  1-1  TO  12 
T ( I  >-T<  I  >-Ct  l 
Tufl“Tw«+T<I) 

NEXT  I 
Tyw-Tyn/ 12 

Nuuegy-FNNueg<  Tum.Egrat  ) 

Nuegu-Nuuegu *Rhoeg 
Ty«l-<  T(  12  >+Tt  1 1  )>T <  10  )  >/3 
Tyn2“(  T<  9  >+T<  8  >+T<  7  )  )/3 
Tun3-<  T(6)  +  T(S>+T(4)>/3 
Ty*4-<  T(  3  )+T<  2  )+T<  l  )  )/3 

Tynl2-(T(  12  )+T<  1 1  HT<  10  )+T(  9  )+T(  8  >+T(  7  )  >/S 

Tyrtl23-<  T(  12  >+T<  11  )+T<  10  >+T<  9  )+T<  8  >+T<  7  >+T(  6  )+T(  5  >+T<  4  >  )/9 

Tln-<  Tout-T<  13  )  )/L06<  <  Tu«-T(  13  >  >/<  Tun-Tout  ) ) 

Qloss-< Tyn-44. IS  1/4 . 39  '  HEAT  LOSS  TO  ATMOSPHERE 

IF  Twm<45  THEN  Qloss-0 

Qtape-Qtape-Qloas 

Qdpi-Qdot/Ai 

Alfai-Qdpl/Tlfi 

FIN  CORRECTION  TO  ACCOUNT  FOR  CONDUCTION  TO  ENDS  OF  TUBE 

Ff  1»<  Kn-Perin-Anarea  )*  .5  _ 

Alfail-Alfai 

Alfai2*Alfail 

Ala«l"<  <  A1  fai2*Perim)/<Km»Anarea  )  >*  .5 
Ff3-Ffl«Alfai2'  .5»<Tyni-T(  13  >  > 

Ff4-Ff l»Alfai2' .S»<  Tun-Tout  J 
Func6»Ff 3*FNTanh<  Alanl •Ulenl  ) 

Func7»Ff4»FNTanb<  Alanl •Ulen2  > 

Func8-Al fai2*Ai •( Tyn-Tavg ) 

Func*,,FuncS+Func7-»'Func8-Qdot 

Ofunc6-< .5*Func6/Al faiZ )+2 . •FfS’Uleni /( 1 . +FNCosh<  2 . -Alan t  »Ulenl  >  > 
0func7-< .5»Func7/Al fai2 )+2, »Ff4»Ulen2/( 1 . +FNCosh< 2 . »Alanl *Ulen2 ) ) 
Ofunc8-Ai •< Tun-Tavg ) 

0funcx»0func6+Qfunc7*0func8 
A1  f at .  “A l  f  a i-~Func.</0 funcx 
IF  A8S(Alfail-Alfai2)’.0S  THEN  SOTO  2079 
Alfai2*Alfail 


2095  Cfeg-<Mueg/Muegu >* . 14 

2096  '  FIND  TEMPERATURE  INCREASE  OUE  TO  UNHEATEO  LENGTHS 

2097  Oelt in»FuncS/<Mdot »Cpeg > 

2098  Oe 1 t out -Func7/ ( Mdot *Cpeg  ) 

2099  Thlm-T(  13  HOeltin 

2100  Thlout-Tout-Oel tout 

2101  Oelhl»<  Thlout -Thl in  > 

2105  IF  I  tube- 1  THEN 

2106  Heat in-Lcui»Oelhl/( L+Lcui ♦Lcuo >  1  HEATED  SMOOTH  INLET  SECTION  FOR  CO 
/N I  FIN  TU8E 

2107  Hea t ou t -Lcuo-Qelhi / ( L*Lcu i *Lcuo >  1  HEATED  SMOOTH  OUTLET  SECTION  FOR 
CU/NI  FIN  TUBE 

2109  Thl in«Thl m+Heat in 

2110  Thlout-Thlout-Heatout 

2111  Oeihl-< Thlout -Thl in > 

2114  ENO  IF 
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2115 

Qhl»Mdot  »Cpeg«Oelhl 

2116 

Qdphl*Qhl/Ai 

2117 

Qdpxx“Func8/Ai 

2113 

Tcol«Thlm+<  Oelhl/4 ) 

2119 

Tco2“Thl in+< Oelhl/2  ) 

2120 

Tco3*Thl in+< 3*0elhl/4) 

2121 

Ttcl»T< 13  )>(  Ld(  0  ) *0i*0elhl/L  ) 

2122 

Ttc2-T<  13  )+<  Ld<  1  >»0i»0elhl/L > 

2123 

Ttc3-T< 13>+<Ld<2)»Qi»0elhl/L> 

2124 

Ttc4-T(  13 )-*•< Ld<  3  )*0i«0alhl/L> 

2126 

Tl«all“<  Thlout-Thlin )/LQ6< ( Tun-Thlin  )/( Tum-Thlout ) )  * 

ALL  d.2.3. 

2127 

Tlnl«(Tcol-Thlin )/L0G< (Tu«l-Thlin)/(Tu«l-Tcol>)  1 

SECTION  l 

2128 

Tl«2-<  Tco2-Tcol  >/L0G<  <  Tw«i2-Teol  >/(  Twn2-Tco2  ) )  ' 

SECTION  2 

2129 

T1ki3-<  Tco3-Tco2  )/L0G(  <  Tuei3-Tco2  >/<  Tun3-Tco3  ) )  1 

SECTION  3 

2130 

Tl«4»<  Thlout-Tco3 )/L06<  <  Tw«4-Tco3  )/( Tur»4-Thlout ) )  1 

SECTION  4 

2131 

Tlnl2-(Tco2-Thlin)/L06<  <Tunl2-Thlm)/(Tunl2-Tco2  >  >  < 

SECTION  1, 

2132 

Tlnl23-(Tco3-Thlin)/l.0S(  <T«Bl23-ThUn)/(Tunl23-Tco3  ) ) 1 

SECTION  1, 

2133 

Alpall“Qdphl/Tlnall 

2134 

Alpl-Qdphl/Tlnl 

2135 

Alp2“Qdphl/Tl«2 

2136 

Alp3*Qdphl/Tl«3 

2137 

Alp4«Qdphl/Tl«4 

2138 

AUocl-Qdphl/<  Twnl-Ttcl  ) 

2139 

A1 loc2”0dphl/<  Tun2-Ttc2  ) 

2140 

Alloe3«Qdphl/(  Tw«3-Ttc3 ) 

2141 

A1 loc4“Qdphl/(  Twm4-Ttc4  ) 

2143 

Alpl2-Qdphl/T'«12 

2144 

Alpl23-Qdphl/Tlnl23 

214S 

Kegl»FNKeg<  (  Tcot*Thl in )/2 .Egrat ) 

* 

2146 

Keg2«FNKog< ( Tco2+Tcol )/2 .Egrat > 

2147 

Keg3“FNKeg<  CTco3+Tco2 )/2 .Egrat ) 

2148 

Keg4»FNKeg(  <  Thlout+Tco3 )/2 , Egrat ) 

2149 

Kegl2»FNKeg<  <  Tco2+Thlin )/2, Egrat ) 

2150 

Kegl23*FNKeg<  (Tco3+Thlin  >/ 2 .Egrat ) 

2151 

Nuna 1 1 • A 1 pa 1 1 »Q i /K  eg 

2152 

Nul"Alpl»Qi/Kegl 

21S3 

Nu2*Alp2»Qi/Keg2 

2154 

Nu3"Alp3#Qi/Keg3 

2155 

Nu4»Alp4»Qi/Keg4 

2156 

Nulocl*AUocl*Qi/Kegl 

2157 

Nulac2«AUoc2»0i/Keg2 

2158 

Nuloc3*Alloc3»0i/Keg3 

2159 

Nuloc  4»A 1  Ioc4*0i/Keg4 

2161 

Nu  1 2#*A  l  p  1 2  *0 1  /  K  eg  1 2 

2162 

Nul23"Alp 123*0 i/Kegl 23 

2163  ' 

LOCAL  NUSSELT  FROM  SHAH/LONOQN  PAGE  128,  EON  240-242 

2164  i 

CALCULATE  XSTAR  FIRST  FOR  EACH  THERMOCOUPLE  LOCATION 

2165  ' 

.  OEFINE  COOLANT  peclet  number 

2166 

Peeg"Preg»Reeg 

2167 

Xs( 1  )-Ldd( 1 )/Peeg  '  FOR  MEAN  UALUES  OF  h 

2168 

Xs<  2  )*Ldd(  2  l/Peeg 

2169 

Xs(  3  )*Ldd(  3  )/Peeg 

2170 

Xa(  4  >»Ldd(  4  )/Peeg 

2171 

FOR  I -1  TO  4 

2172 

Xstar( I)-Ld( 1-1 >/Peeg  '  FOR  LOCAL  VALUES  AT  T/C  POSITIONS 

2173 

Nua( I )-1.302»Xstar< I )'( -1/3J-.S 

2175  ' 

OEVELOP ING  HYORO  ANO  THERMAL  80UN0ARY  LAYER  ( S/L  Pi 48, 

EON.  258) 

2176 

Xxl«<  l  +  <  <  220/PI  )»Xstar<  I  )  >*<  -10/9  )  >*  <  3/5  ) 

2177 

X*2-< 1  ♦(  Preg/ . 0207 >*(2/3 ) )* .5 

2178 

Xx3-PI/< 1 15 . 2  *X  at  ar( I >  > 
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Xx4-<  t+(Xx3/(Xx2»Xxl  >  )'< S/3) >*(3/10) 

XxS-< l+< (220/PI  )«Xstar( I  )  >*(  -10/9 ) )“( 3/10 ) 

Noaa( I  )-(X.,4»S.3S4.XaS  )-l 
IF  Xstar(  I  K. 00005  THEN 
Nux( I  >«1.302»Xstar( I >“<-l/3)-l 
ENO  IF 

IF  Xatar< I >>.0015  THEN 

Nux<  I  >-4.364-t-8.SB»<  l000»Xatar<  I  )  )*(  -  .506  >«EXP(  -41»Xstar<  I  )  ) 

ENO  IF 
NEXT  I 

IF  Inaert-0  THEN 

PRINT  USIN6  "3X , ' "  LOCAL  NUSSELT  VALUES'  *  * 

PRINT  USING  "3X , *  *  Exp.  Nux(TLM)  Exp.  Nux(BULK)  Pred.  Nux< 1  ) 


Prad.  Nux(2>  Xstar(tc)" 


2206 

2207 
atari  1  ) 

2208 
atar<  2  ) 

2209 
star<  3 ) 

2210 
Star! 4  ) 
2211 

2215 

2216 


PRINT 

PRINT  USING  *3X.“P0S  1  “  .  S<  MZ  .  30E.  SX  >' .  Nul  .Nulocl  ,Nux  (  1  )  ,Nuss(  1  ) ,  X 

PRINT  USIN6  '3X,”P0S  2  ' * ,5( MZ .30E.SX >' ;Nu2,Nuloc2 ,Nux( 2 ) ,Nuaa( 2  i  ,X 

PRINT  USING  *3X.”P0S  3  '  *  ,S<  MZ .  30E.SX  >' ;  Nu3,Nuloc3  ,NUx<  3  )  ,Nuss<  3  )  ,X 

PRINT  USING  "3X,"P0S  4  “ ' ,S( MZ .3DE ,SX )' ;Nu4,Nuloc4,Nijx( 4 ) ,Nusa<  4 ) , X 

PRINT 

ELSE 

PRINT  USING  '3X,  “POSITION  *  -  Exp.  Nux  ( TLM  >  Exp.  Nux  (SULK) 


Xatar  at  tc  poa. *  *  * 


2217 

2218 
star( l  ) 

2219 
atar( 2  ) 

2220 
3tar<  3  ) 
2221 
star( 4  ) 
2222 
2223 
2230 

2232 

2233 


PRINT 

PRINT  USING  '3X , " “POSITION  1  (LOCAL)  -  ”  ,3(MZ.30£,SX r ;Nul .Nulocl ,X 

PRINT  USING  '3X,  "POSITION  2  (LOCAL)  -  ”  ,3(MZ.3DE,5X >' ;Nu2,NiHoc2,K 

PRINT  USING  -3X,  "POSITION  3  (LOCAL)  -  ”  , 3( MZ . 30E , 5X  ) ' ; Nu3 , Nuloc3 , X 

PRINT  USING  '3X. "POSITION  4  (LOCAL)  -  "  , 3( MZ  .3DE  ,  SX  ) ' . Nu4 , Nuloc* . X 

ENO  IF 
PRINT 

FOR  1-1  TO  4 
IF  Inaert-0  THEN 

8bl-.0668»Reeg-Preg/Ldd< I  ) 

8b2-l  +  .  04»( Reeg*Preg/Ldd( I  )  . 6665 

Nueg( I  >-3.66  +  <  8b 1 /8b2  ) 

ENO  IF 

IF  Inaert-l  OR  Inaert-4  THEN 

Nueg<  I  )-Cia(  Insert  )•(  l +5 . 484E-3»Preg" . 7»(  Res/Hcd  >  '1.25  )  *  .5 
•  Nuo th-. 383 »Pr eg  ' . 35 •( fles/Hod )* .622 
ENO  IF 

IF  Insert-2  THEN 

Nueg( I  )-Cia( Insert  )•( Reeg  1 . 65  )»Preg" . 46 
ENO  IF 

IF  Insert-3  THEN 

Nueg< I  )-C ia( Insert  ) •( Reeg" . 76 ) »Pr eg" . 46 
ENO  IF 
SEEP 

CALCULATED  VALUE  OF  INSIOE  HT.  TRANSFER  COEFF.  FROM  CORRELATIONS 
Hi«Nueg< 4  )»Keg/Oi 

MEAN  NUSSELT  NUMBER  FROM  SHAH/LONOON  PAGE  128.  EQN  245-246 
IF  Xs(  I  >■'  .03  THEN 
Nun<  I  >-l .  953»Xs(  I  )"(  -1/3  ) 


22S6 

F.SE 

2258 

Num< I  )— 4 . 3S4+( ,0722/Xs< I > > 

2259 

END  IF 

2262 

NEXT  I 

22E3 

PRINT  USING  "3X , “ "POSITION 

t  -  1 

Xstar 

at  section”*” 

2254 

PRINT 

2256 

PRINT  USING  ”3X, ""POSITION 

1  (MEAN)  - 

,Xs(  1 ) 
2257 

PRINT  USIN6  ”3X,  "POSITION 

12  (MEAN)  - 

),Xa<2) 

2258 

PRINT  USING  *3X, "POSITION 

123  (MEAN)  ■ 

3  >  ,Xs(  3  > 
2269 

PRINT  USING  ”3X , ” “POSITION 

1234  (MEAN) 

(  4  )  ,Xs(  4  > 
2270 

PRINT 

2271 

PRINT 

2272 

PRINT  USIN6  " 10X , “ ‘Mass 

flow  rate  < 

Exp .  mean  Nu  Pr»d.  mean  Nu 


■  ,3<  MZ.30E.5X  >” ; Nul23 ,  Num( 


2273 1 

2274 

MZ.30E” iHod 

2275 


PRINI  UMNb  iWx,  insioe  tuoe  uia.  irw  - 

IF  Insert-1  THEN  PRINT  USIN6  *10X,"180  QEG  OVER  Oia.  < HCO ) 


nt.ouc  inou 

2275  PRINT  USING  '  10X , * ’ Inlat  temperature  <C>  -  .MZ.30E” ; T(  13 > 

2276  IP  Inaert-I  OR  lnsert-4  THEN 

2277  PRINT  USING  ”10X,"OELTA  Tape  Thickness  -  " .MZ.3DE” ;0elta 

2278  END  IF 

2279  PRINT  USING  ”10X,"OELT  temp  Oif.  <C>  -  .MZ.30E” ;0elt 

2280  PRINT  USING  ”10X,”H.  T.  to  coolant  (U)  -  *  *  ,flZ .  3DE" ;  Qdot 

2281  PRINT  USING  ”10X,"H.  T.  from  heater  tape  (U)  -  * • ,nZ .30E" . Qtape 

2282  PRINT  USING  ” 10X , * ”Prandt 1  number  -  "  ,MZ . 3DE” . Preg  - 

2283  PRINT  USING  ” 10X , ” 'Reynolds  number  -  *  * ,MZ . 3DE” ■ Reeg 

2284  IF  Insert-l  OR  Insert-4  THEN-PRINT  USING  ’  10X  ,  " -Reynolds  number  H 

18  S  -  "  .MZ.30E’ :Res 

2285  PRINT  USING  ”10X,"Hi  from  correl.  (U/m'2.K)  -  "  ”  ,MZ ,30E” . Hi 

2286  PRINT  USING  "l0X,"”Hi  from  heat  bal .  <U/m“2.K>  »  ” " , HZ . 30E “ ; A 1 f a i 

2 

2287  IF  Insert-0  THEN  PRINT  USING  ”l0X,"Nu  for  no  insert  (Hausen) 

■ , MZ.30E” ;Nueg<  4  ) 

2288  IF  Insert-1  OR  'nsert-4  THEN  PRINT  USING  ”10X,”’Nu  for  twisted  ta 
pe  ( Bergles2 >“ ” , MZ . 30E” ; Nuo th 

2289  IF  Insert-l  0°  Insert-4  THEN  PRINT  USING  -10X,"Nu  for  twisted  ta 
pe  ( 8ergles3  ) ”" ,MZ . 30E” ; Nueg( 4 > 

2290  IF  Insert-2  THEN  PRINT  USING  ”10X,””Nu  for  large  Heatex  (Macconel 
” " ,MZ . 30E” . Nueg<  4  ) 

2291  IF  Insert-3  THEN  PRINT  USING  ”10X,'“Ni.  for  small  Heate.*  >  Haooone  ) 
”” .MZ.30E” .Nueg(4 > 

2292  IF  Insert-0  THEN  PRINT  USING  “10X,””Nu  from  Shah/London  *  ““• 
flZ  .  30E”  •  Num(  4  ) 

2293  PRINT  USING  "  10X  ,  * ” X s t ar  for  Heated  length  -  ' * . HZ . 2DE ' . X si  4 > 

229<t  PRINT  USING  "  10X  , "  "Experimental  Nusselt  no.-  ” ' , MZ . 306 " • Numa 1 1 

2295  PRINT 


2296 

PRINT 

USING 

”3X. 

” "Ual 1 

t/c 

12 , 1 1 , 10  ( pos . 

1 ) 

(Cl¬ 

■ "  ,3<  000.00, 3X 

)” .  T<  12 
2297 

),T(  11  ),T(  10) 
PRINT 

USING 

”  3X  , 

""yall 

t/c 

9.8.7 

( pos . 

2  ) 

io  * 

”’ .3(000. DO, 3X 

>”  .  T(  9  >,T<  8  >,r<7> 

2298  PRINT 

USING 

”  3X  , 

'■Wall 

t/c 

5,5.4 

<  DOS  . 

3  ) 

(Cl¬ 

3(000. 00, 2X 

)”  ,  T(S  ) 

2299 

,T(S  ),T(  4  ) 

PRINT 

USING 

"  3X  , 

,  ”  "Ual 1 

t/c 

3.2.1 

<  DOS  . 

4  ) 

io- 

" "  ,  3(000. 00, 3X 

)■  .  T(  3  ),T(  2  >,T<  1  > 

2300  PRINT  USING  ”"3X  , "  “  Av  .  inside  wall  temp  (C)  -  ” ”  ,  HZ  .  20E ”  ■  Twm 

2301  PRINT 
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POS  123 


Nu( 1  )»Nul 

Nu<  2  )-Nul2 

Nu( 3  >“Nul23 

Mu<  4  )"Nunal  1 

IF  Insert*®  THEN 

FOR  1-1  TO  4 

XxxK  I  >-LQG(Nu<  I  )  ) 

Xx.<2<  I  i-L0G<  Xs(  I  ) ) 

X.<x3(  I  )*LQG<  Nueg<  I )  > 

Xxx4< I )"LOG<  Nun! 1  ) ) 

NEXT  I 

PRINT  USIN6  '3X ,  * ‘MEAN  VALUES  AT 
POS  1234- ■’ 


2302 

2303 

2304 

2305 
2307 

2306 

2310 

2311 

2312 

2313 

2314 

2315 
2 

2317 
xxl<2), 

2318 
xx3<  2 ) , 

2319 
xx4<  2  > , 

2320 
xx2<  2 ) , 

2321 

2323 

2324 

2325 

2326 

2327 

2328 

2329 

2  POS  123  POS  1234*** 

2331  PRINT  USING  *6X,**EXP.  L06« Nu/Pr* .46 > 
Xxxl(2  >,XxxK3).Xxxl(4) 

2332  PRINT  USING  *SX.**C0R.  LOG(Nu/Pr* .46 ) 
Xxx3(2),Xxx3<3>,Xxx3(4) 

2334  PRINT  USING  "6X , *  *LQG<  Re  > 

Xxx2<  2  ).Xxx2( 3  ).Xxx2( 4) 

2335 

2336 

2337 

2338 

2339 

2340 

2341 

2342 


PRINT  USIN6  *6X,*‘EXP.  LQG(Nu) 

Xxxl( 3  )  ,Xxxl( 4  ) 

PRINT  USING  *6X,**C0R.  LOG(Nu)  < H > 
Xxx3<  3  ) , Xxx3( 4  ) 

PRINT  USING  "SX.-’COR.  LOG<Nu>  <S/L> 
Xxx4( 3  > ,Xxx4( 4  ) 

PRINT  USING  *6X , * “LQ6( Xstar  ) 
Xxx2<3),Xxx2(4) 

ENO  IF 

IF  Insert»2  OR  Insert-3  THEN 
FOR  1-1  TO  4 

XxxK I )-L0G<Nu< I )/Preg* .46  ) 

Xxx2< I >"LQG< Reeg ) 

Xxx3<  I  >»LOG<Nueg<  1  )/Preg*.46  ) 

NEXT  I 

PRINT  USING  “3X , * ‘MEAN  VALUES  AT 


POS  1  POS  1 

‘,4(HZ.30E,4X)";Xxa1( 1 > , X 
’  ,4<fiZ.30E,4X  )';Xax3( 1  ) , X 
\4<MZ.3DE,4X  > * ; Xxx4( 1 ),X 
“  ,4(«Z.30E,4X  )* ; Xx x2 ( 1  ),X 


POS  1  POS  1 

■  *  ,4(nZ  .3DE.4X  >* ; Xxx 1!  l  > , 
•"  ,4(MZ.3QE,4X  >*  ,Xax3<  l  ), 
.4(MZ.30E,4X  )*  ;  Xx.<.2<  1  !. 


ENO  IF 

IF  Insert"!  OR  Insert«4  THEN 
FOR  1-1  TO  4 

X  •. .<  1  <  I  )-LOG<Nu(  I  >/Preg " .  35  ) 

X x ;< 2 (  I  >«LOG(  Res/Hod  ) 

Xaa3( I  )"L0G( Nuoth/Preg" .35 ! 

Xxx4<  I )«LOG<  Nu< I  )  ) 

Xxx5(  I >“LQG(  Preg*<Res/Hod>'l.78  ) 


2343 

Xxx6<  1  )*L0G(Nueg(  I  ) ) 

2344 

NEXT  I 

234S 

PRINT  USING  *3X , * 

"MEAN 

values  at 

POS  I 

POS 

2 

POS  123  POS  1234* 

'  - 

2346 

PRINT  USING  *SX,* 

*EXP . 

LQG<  Nu/Pr' . 35  > 

*  " * ,4<  HZ . COE . 4 x 

'  *  ■  X  •  •  K  1  ' 

Xxx 1 ( 2  )  ,Xxxl( 3  > , XxxK 4 > 

2347 

PRINT  USING  "6X , * 

'COR. 

LOG! Nu/Pr" .35 ) 

•  **,4<nz.30E,4X 

.  Xaa3<  1  ? 

Xxx3( 2 

> ,  X  <  x3<  3  ),Xxx3<  4  ) 

2348 

PRINT  USING  *6X , * 

'LOG!  Res/y > 

’■ ,4(MZ.30E,4X 

1*  ;  X  x  x  2 !  1  > 

Xxx2<  2 

)  ,Xxx2<  3  ).Xxx.2<  4  ) 

2349 

PRINT 

2350 

PRINT  USING  "GX , " 

*EXP  . 

LOG( Nu  )  * 

*“ ,4<MZ.:0E.4X  : 

I  ■■  ■  X  A ..  4  <  1  j 

X  x  x  4  (  2 

)  ,  Xx  .4(  3  )  ,  X/x4<  4  ) 

2351 

PRINT  USING  *6X , " 

*COR . 

L0G<  Nu  ) 

■*,4<riZ.30E.4x  ; 

• "  •  x  x  •  6  (  1  1 

Xxx6(  2  > ,Xxx6<  3  )  ,Xx*6<  4  ) 

23S2  PRINT  USING  "6X , " *LOG( Pr (Res/y  )* 1 . 78  ) 
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'*  ,4<  MZ.30E  ,4X  )* .  Xa  .5<  l  ), 


,4<NZ.30E,4X  ) * . Nu 


XxxS(  2  >,X*xS<  3  >,X.<xS(  4  ) 

2353  ENO  IF 

2354  IF  Inaart-0  ANO  Reeg>230O  THEN 

2355  Nutu».027»Reeg* .8*Preg‘ .3333»Cfeg 

2356  PRINT  USING  "SX , “ "TUR8ULENT  S/T  CORRELATION  -  "  * 
tu 

2357  Xxx7-L0G(Nutu ) 

2358  PRINT  USING  “SX , "  “TURB .  COR.  LOG(Nu)  -  *“ ,4(MZ.30E,4X  )■ ,Xxx7 

2359  ENO  IF 

2360  Okaccpt-1 

2361  8EEP 

2362  INPUT  'OK  TO  ACCEPT  THIS  SET  < 1-OEFAULT-YES,  0-NO >7* .Okaccpt 

2363  IF  Okaccpt-1  THEN  OUTPUT  9Fi le ; Pc< ♦ > , E«f ( ♦ > , Tp< • > 

2364  IF  Im-0  THEN 

2365  Okrpt-1 

2366  BEEP 

2367  INPUT  'WILL  THERE  BE  ANOTHER  RUN  < 1-VES-OEFAULT ,0-NO  )* .Ofcrpt 

2368  IF  Okrpt-1  THEN  1701 

2369  ELSE 

2370  IF  JseUNsets  THEN  1701 

2371  ENO  IF 

2372  ASSIGN  SFile  TO  « 

2373  IF  lout-1  THEN  ASSIGN  0FUe  TO  • 

2374  SUBENO 

2375  OEF  FNGrad(T) 

2376  Grad— 3-.  8778S7E-S-2»4 . 7 1 428S7E-8»T 

2377  RETURN  Grad 

2378  FNENO 

2379  OEF  FNKcu(T) 

2380  '  OFHC  COPPER  250  TO  300  K 

2381  Tk-T+273. 15 

2382  K-434- . 1 12*Tk 

2383  RETURN  K 

2384  FNENO 

2385  OEF  FNNueg<  Tc , Egr  ) 

2366  '  RANGE  OF  VALIDITY:  -20  TO  20  OEG  C 

2387  Tk-Tc+273 . IS 

2388  Nul-7. 1 196507E-3-Tk«( 7.4863347E-5-Tk •< 2 .6294 943E-7-Tk *3.083332 9E- 10 > > 

2390  Nu2-4 . 9237638E-3-Tk •<  4. 92 13912E-5-Tk •< 1 .6437S346-?-TkM . 333333  IE- 10 > ) 

2395  Nu3-8 . 6586293E-3-Tk ♦< 8 . 8837302E~5-Tk  »( 3. 0495032E-7-Tk  »3 . 49999966-10  >  ) 

2400  A2-< Nu3-2»Nu2+Nul )/200 

240S  Al»<  Nu2-Nul-940»A2  )/  10 

2410  A0-Nul-42»Al-l764»A2 

2415  Nu-A0+Egr»< Al*Egr»A2 > 

2420  RETURN  Nu 

242S  FNENO 

2430  OEF  FNCoeql Tc , Egr  ) 

2435  1  RANGE  OF  VALIOITY  0  TO  20  OEG  C 

2440  Tt -Tc*273 . 15 

2445  Cpl-l .6701550E+3+Tk»6.3 

2450  Cp2»l .  4748i2SE+3+’Tk  »S  .  25 

2455  Cp3*9.5300500E+2  +  Tk.»7.3 

2460  A2-< Cp3-2»Cp2*Cpl 1/200 

2465  Al«< Cp2-Cpl-900»A2  1/10 

2470  A0-Cpl-40.AI-1600*A2 

247S  Cp«A0+Ei3r*(  Al+Egr»A2  1 

2480  RETURN  Co 

2485  FNENO 

2490  OEF  FNRhoegt T , Egr  ) 

2495  Rol-L .0607093E+3-T.<  3.7031283E-l*T»4.0837l83E-3  > 


151 


2500  Ro2*l .0748272£*3-T»<  4.4266195E-1*T«4 . 0939706E-3  ) 

2505  Ro3-l . 088S934E+3-T* < 5 . 73S56S3E- l*T«6 . 128 140SE-3  ) 

2510  A2-<Ro3-2»Ro2+Rol >/20O 

2515  Al-<Ro2-Rol-90O*A2 >/10 

2S20  A0-Rol-40»AI-1600»A2 

2525  Ro»A0+Egr*<  Al+£gr»A2 ) 

2530  RETURN  Ro 

2535  FNENQ 

2540  OEF  FNPregl T.Egr > 

2545  Pr«FNCpeg<  T.Egr )»FNNueg( T,Egr  >»FNRoeg<  T,Egr  >/FNKeg< T.Egr > 

2550  RETURN  Pr 

2555  FNEND 

2SS0  OEF  FNKeg(Tc.Egr) 

2565  !  RANGE  OF  VALIDITY;  -20  TO  20  0E6  C 

2570  Tk-Tc+273.15 

2575  Kl-2.2824708£-l+Tk»<5.5989286E-4+Tk*3.57l4286E-7 ) 

2580  K2-2. 58466 ISE-l+Tk  ♦< 2. 397857 l£-4+Tk*7 . 1428S71E-7 > 

2585  K3«3 .2l38932E-l-Tk»<  3 .00428S7E-4-Tk»l . 42857 14E-6  ) 

2590  A2-<K3-2»K2+K1 )/200 

2595  A1-(K2-K1-900*A2)/10 

2600  A0»K1-40»A1-1600*A2 

2605  K»A0+Egr»<  Al+Egr»A2 ) 

2610  RETURN  K 

2615  FNENO 

2620  OEF  FNT«nh<  X ) 

2625  P«EXP<X> 

2630  Q-l/P 

2635  Tanh-<  P-0  )/<  P+Q  > 

2640 _ RETURN  Tanh 

264S  FNENO 

2650  OEF  FNTvav(V) 

2655  COM  /Cc/  C<7) 

2660  T«C<0> 

2665  FOR  1-1  TO  7 

2670  T-T+C< I >*V"I 

2675  NEXT  I 

2680  RETURN  T 

2685  FNENO 

2690  OEF  FN8eta<  T  ) 

2695  Roo-FNRho(T+. 1 > 

2700  Ron«FNRho< T-. I ) 

2705  8eta»-2/<  Roo+Ron )♦< Rop-Ron )/ .2 

2710  RETURN  Beta 

27  IS  FNENO 

2720  OEF  FNPaat(Tc) 

2725  '  0  TO  80  deg  F  CURVE  FIT  OF  Paat 

2730  T  f»l . 8»Tc+32 

2735  Pa*5 . 94S525+T  f •( . 15352082+Tf •< 1 . 4840963E-3+T f •9.61S0671E-6 ’ ' 

2740  Pg«Pa-l4.7 

2745  IF  Pg/0  THEN  >  ♦  •PSIG,— in  Hg 

2750  Psat-Pg 

2755  ELSE 

2760  Psat«Pg»29 .92/14.7 

2765  ENO  IF 

2770  RETURN  P sat 

277.5  FNENO 

2905  SU8  Purge 

2910  8EEP 

2915  INPUT  'ENTER  FILE  NAME  TO  8E  OELETED'  .File* 

2920  PURGE  File* 
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2925  SOTO  2910 

2950  SUBENO 

2935  SUB  Wilson 

2940  COM  /Will/  Ooa<  4>,Oia<4),Kma(4),Iact ,Qraot  <  4  )  ,Ld<  3  > 

2945  COM  /Uil2/  Oelta, Isat .Nsets ,Hod,Cia( 4 > .Aleaat 3 > 

2950  DIM  Enf( 17),T( l7),Xa( 100),Ya( l00),Pc(4),To(4) 

2955  BEEP 

2960  INPUT  'PLEASE  RE-ENTER  NAME  OF  FILE', File* 

2965  ASSIGN  9File  TO  File# 

2970  INPUT  -0-SMOOTH  COPPER  TUBE,  1-CUNI  FIN,  2-KOROOENSE,  3-FIN  TITAN', Itube 

2975  8EEP 

2976  INPUT  ‘0-NO  INSERT,  1-  TWISTED  TAPE  INSERT,  2-HITRAN1 ,  3-HITRAN2* , Insert 

2977  BEEP 

2980  INPUT  “SIUE  A  NAME  FOR  XY  FILE’.Xy* 

2985  CREATE  BOAT  Xy*,20 

2990  ASSIGN  9Xy  TO  Xy* 

2995  L-1.2192 

3000  Oo-Ooa< Itube) 

3005  Oi-Oia< Itube) 

3010  Kn-Kna< Itube ) 

3015  Ax-PI«0iA2/4  !  Cross-sectional  area 

3020  Ao-PI»0o»L 

3025  Rn-Oo»LOG<  Oo/Oi  )/(  2*Kpi  ) 

3030  ! 

3035  1  Initial  values 
3040  Tf-Tsat 

3045  Alpa-.6SS 

3046  IF  Insert-0  THEN 

3050  Ci-1.0. 

3051  ENO  IF 

3052  IF  Insert-1  THEN 

3053  Ci -5. 172 

3054  ENO  IF 

3055  IF  Insert -2  THEN 

3056  Ci -.226 

3057  ENO  IF 

3058  IF  Insert -3  THEN 

3059  Ci-. 063 

3060  ENO  IF 

3062  G-9.31 

3063  Ibeg-0 

3065  Iend-0  'CHANGE  TO  4 ,  IF  FIVE  TUBES  IN  8UNOLE 

3070  IF  Iact-10  THEN 

3075  Ibeg-Iact 

3080  lend-Iact 

3085  ENO  IF 

3090  ' 

3095  FOR  I-Ibeg  TO  lend 

3100  3a»0 

3105  Sy-0 

3110  Sas-0 

3115  5.*y-0 

3120  Jset-0 

3125  ASSIGN  UFile  TO  File* 

3130  ENTER  9Fi le ; Otg* , I  tube .Egrat ,Qdl , Dd2 ,0d3 ,0d4 ,0d5 

3135  ENTER  9F  i  le ; Pel •  ),Enf< • >,Tp( •  ) 

3140  FOR  J-0  TO  17 

3145  T< J  )-FNTvsv< Enf ( J  )  ) 

3150  NEXT  J 

3155  Tvao*<  T(  0  )+T(  1  )  +  T(  2  )  )/3 
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3160 

3161 
3165 
3170 
317S 
3130 
3185 

3190  1 

3191 
3195 

3200 

3201  f 

3202  ! 

3204  * 

3205  ! 

3206  ' 
3210  ' 
3215 
3220 
3225 
3230 
3235 
3240 
3245  ! 
3250 
3255 
3260 
3265 
3270 
3275 
3280 
3295 

3286 

3287 
3298 

3289 

3290 

3291 

3293 

3294 

3295 

3296 

3297 

3298 

3299 

3300 

3302  ' 

3303  ' 
3335 
5310 
3315 
3320 
3325  i 
3330 
3335 
3340 
3345 
33S0 
3355 
3360  ' 


Tliq-<T(3>+T<4>)/2 

Tliol-TlS) 

IF  I sat-0  THEN 
Tsat»Tl lql 
ELSE 

Tsat-Tvap 
ENO  IF 

Qfed-FNGradl  T(  I  *6  ) ) 

Tout  •  <  T<2»  1  +  10  >+T(  2*1  +  1 1  >  >'2 
Qelt-Tout-T< 1+6 ) 

Tavg-T  < 1+6 )+Oel t • .5 
IF  I-Iend  THEN 
Oalt-TU0>-T<9> 

Tavg-(T(10)-T(9))/2 
ENO  IF 

Uater/Ethylene  Glycol  Mixture  Properties 
Rhoeg-FNRhoag(  Tavg.Egrat ) 

Nueg-FNNueg( T  avg , Egr a  t ) 

Mueg-Nueg*Rhoeg 
Cpeg-FNCpeg<  Tavg.Egrat ) 

Keg-FNKeg<  Tavg.Egrat ) 

Preg-Cpeg+Mueg/Keg 

Mdot-FNFncal  ( I , T( 1+6  )  ,Pc< I > > 

Veg-Mdot/<Rhoeg+Ax > 

Roe q -Veg  +0  i ./  Nu  eg 

Res"4«ridot/(  Piueg*<  PI»0i-4«0elta  >  > 

Odot -Mdot -Cpeg-Oe 1 t 
Qdp-Odot/rto 

Lntd-0alt/L0G<  <Tsat-T< 1+6  > >/< Taat-T( 1+6  >-0elt  ) ) 

Uo-Qdp/Lntd 

IF  Insert-0  THEN 

Bbl-.0668»<  Oi/L  )«Reeg*Preg 

BbZ«l+.04»<  (Oi/L >-Reeg+Preg>".6666 

Onega-3. 66+<BbI/8b2  > 

ENO  IF 

IF  Insert-!  THEN 

Onega -< 1+S . 484E-3»Preg" . 7»< Res/Hod  )  *  1 . 25  > * .5 
END  IF 

IF  Insert-2  THEN 
Onega-CReeg" .65  )»<Preg"  .46  > 

END  IF 

IF  Insert-3  THEN 

Qnega-<  Reeg" .76  >•<  Preg* .48  > 

ENO  IF 

R-114/R113  Properties 
Hfg-FNHfg<  Tsat  > 

K  f -FNK ( Tf  ) 

Rhof-FNRho( Tf  ) 

Muf-FNMu< Tf ) 

F*<Kf3»Rhof"2fG»Hfg/<Muf*Qo»Qdo  )  >“  .33333 

Ho-P loa»F 

Tuo-Tsat-Qdp/Ho 

Tf-Tsat/3+2*Tuo/3 

Y»(  1  /Uo-fln  ) -F 

X»0o»F/(  Keg«0nega  ) 

PRINT  "OMEGA-" : Onega ; "F-* .F, "X-" .X, *Y»" .Y 
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3355 

3370 

3375 

3380 

3385 

3390 

3395 

3400 

3405 

3410 

3415 

3420 

3425 

3430 

3435 

3440 

344S 

3450 

3455 

3460 

3465 

3470 

3475 

3480 

3485 

3490 

3495 

3500 

3505 

3510 

3515 

3520 

3525 

3530 

3535 

8270 

8280 

8290 

8300 

8310 

8311 
8320 

8470 

8471 
8480 
8490 
8500 
8510 
8520 
8530 
8540 
8550 


XaUset  >-X  '  INEFFICIENT  (  MQOIFY  LATER) 

Ya<  J set  W 

Sx»Sx>X 

Sy-Sy+Y 

SxS"SxS+X»X 

Sxy»Sxy+X»Y 

Jset»Jset+l 

IF  Jset'Nsets  THEN  3135 
ASSIGN  8File  TO  . 

Slo{je*(Nset3*Sxy-Sx»Sy)/<Nsets«Sxs~Sx*2) 

Intcpt“< Sy-Slope»Sx )/Nsets 

Cic«l/Slope 

Alpac“l/ Intcpt 

Cerr- A8S( <Ci-Cic )/Cic ) 

Aerr*A8S( (  Alpac-Alpa )/Alpac ) 

IF  Cerr> .00 1  OR  Aerr,.001  THEN 
AIpa«( Alpa+Alpac  )• . 5 
Ci*<  Ci+Cic )• .5 

PRINT  "CIC»‘;Cic; “ALPA»* ;Alpa 
SOTO  3100 
ENO  IF 
SEEP 
BEEP 

PRINTER  IS  1 

PRINT  "CIC**;Cic; “  ALPA»* ; Alpa  . 

Cia< Insert  )«Cic 
Alpaa< I  )“Alpac 
PRINTER  IS  701 
FOR  J-0  TO  Nsets-l 

OUTPUT  9Xy;Xa( J),Ya( J) 

NEXT  J 

PRINTER  IS  l 
NEXT  I 

ASSIGN  0Xy  TO  • 

SU8EN0 

OEF  FNFncal< T ,Pc  ) 

rtdo t 1S«. 001 9»Pc-. 0098 
Mdo t 0«. 00206 «Pc-. 004 14 
Mdot24-.0022S»Pc-. 00329 
IF  T‘0  THEN 

Mdt-<  (  T-t-15  >•(  P1dot0-Hdot  IS  )/  15  >+Mdot 15 
ELSE 

Mdt»<  T*< Hdot24-Mdot0  )/24 )+Mdot0 
ENO  IF 
RETURN  Ndt 
FNENO 

OEF  FNCosM  X  ) 

P-EXP(X) 

Q«EXP( -X  ) 

Cosh» . S»( P+Q  > 

RETURN  Cosh 
FNENO 
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